M.Sc. Physics - I Year: Dkp11: Classical and Statistical Mechanics
M.Sc. Physics - I Year: Dkp11: Classical and Statistical Mechanics
PHYSICS – I YEAR
DKP11 : CLASSICAL AND STATISTICAL MECHANICS
SYLLABUS
1.Classical Mechanics (3rd ed.,2002) by H. Goldstein, C.Poole and J. Safko, Pearson Edition
2. Classical Mechanics - J. C. Upadhyaya- Second Edition-2005-Himalaya Publishing House
3.Classical Mechanics - G. Aruldhas-2008-PHI Learning Pvt.Ltd.
4.Classical Mechanics-A Text Book-Suresh Chandra-Narosa Publications
5.Statistical Mechanics by R. K. Pathira (2nd edition)
6.Statistical Mechanics by R.K. Pathira and P.D. Beale (3rd edition)
7.Statistical Mechanics by K.Huang
8.Statistical Mechanics by L.D.Landau and I.M.Lifshitz
I= Ldt
t1
(or) I = Ldt 0
t1
T q , q V q dt.
t2
j j j
t1
T q j , q j V q j dt = 0
t2
t1
t2
T T V
t j q j j q j j q j j dt = 0
q q q
1
T V
q j dt T q j dt 0
t2 t2
t j q j q j t j q j
1 1
t j q j q j t q j dt
1 1
1
t1
t1
dt q j
There is no coordinate variation at end points and hence
t2
q j t1
0.
t j q j q j
t j dt q j q j dt 0
1 1
t2
T V d T
t q j q j dt q j
q j dt 0
1
each δqj are independent of each other, the coefficient of every δqj should be equated to
zero. And we get
T V d T
0
j
q q j dt q j
d T
q
(T V ) 0
dt q j
j
V is not a function of q j and therefore
d (T V )
(T V ) 0
dt q j q j
d L L
0 J = 1, 2,......n.
dt q j q j
This set of equations is called Lagrange’s equations of motion and each is a second
order differential equation in terms of the time as independent variable.
1.2 DERIVATION OF LAGRANGE’S EQUATIONS FROM HAMILTON’S
PRINCIPLE:
The Lagrangian L is a function of generalized co-ordinates q j ′s and generalized
velocities 𝑞 j ′s and time t.
Manonmaniam Sundaranar University, Directorate of Distance & Continuing Education, Tirunelveli.
4
That is L L[q j (t ), q j (t ), t ]
I= L[q j (t ), q j (t ), t ]dt
t1
t2
If the Lagrangian does not depend on time t explicitly, then the variation δL can be written as
n L L
L q j q j
j 1 q q j
j
L L
t2 t2 2 t
t Ldt t j q j
q j dt t j q j q j dt 0
1 1 1
Ldt 0
t1
L L
t2 2 t
Therefore t j q j j
q dt t j q j q j dt 0
1 1
d
q j (q j )
dt
Integrating by parts the second term, we get
t2
t2
L L t2
d L
q j dt
t j q j j
q dt
j q
q
t1 t1 j dt q j
j
1 j
L d L
t2 t2
q j dt 0
t j q j j t j dt q j
q dt
1 1
(or)
d L
q j dt L q j dt 0
t2 t2
t j dt q j
t j q j
1 1
For holonomic system the generalized co-ordinates δq j are independent of each other.
Therefore the coefficients of each δq j must vanish. And we get
d L L
0 J = 1, 2,......n.
dt q j q j
This set of equations is called Lagrange’s equations of motion and each is a second
order differential equation in terms of the time as independent variable.
1.3 PRINCIPLE OF LEAST ACTION:
The time integral of twice the K.E is called the action. The principle of least action
states that
t2
2T dt 0
t1
t2
p j q j dt 0
j
t1
q j t
d q j d
q j q j t
q j t
Since q j d and q j d q j t
Any function f = f (q j , q j , t )
f f f
f = q j q j t
j q q j t
j
f
q j q j t f q j qj t f t
j q j q j t
f
f t
t
= t ....(1)
t
Here operation and time differentiation cannot be interchanged.
Proof:
t2
A = p q
t1
j
j j dt
t2
= ( L H ) dt
t1
t2
= L dt H (t
t1
2 t1 ) [since H is conserved] ....(2)
t2
A = L dt H (t 2 t1 )
t1
t2 t2
= L dt H t ....(3)
t1 t1
t2
To Solve L dt :
t1
t2
Let L dt
t1
=I so that I L
t2
L dt cannot be zero.
t1
t2 t2
L L
L dt q j q j dt
t1 t1
j
q j q j
d L
q j L d (q j ) dt
t2
=
j dt q q j dt
t1 j
L d L
[Since from Lagrange’s equation of motion ]
q j dt q j
t2 t2
d L
L dt t j dt q j j dt
q
t1 1
Putting q j q j q j t
t2
d L L
t2
L dt q j q j t dt
j dt q q j
t1 t1 j
t2
L L
= q j q j t
j q q j
j t1
L
t2 t2
L dt q j t |
t1
j q j t1
t2
= p j q j t |
j t1
= H L p j , q j t |
2 t
j t1
t2
Thus A p j , q j dt 0
t1
Here p, q
t1
j Hamiltion’s characteristic function
(i) Principle of least action in terms of arc length of the particle trajectory.
Let a system contain only one particle of mass m
2
Kinetic energy T m
1 ds
2 dt
ds element of arc traversed in time dt
1
m
2
dt = ds
2T
The principle of least action
1
t2 t2
m
2
2T dt 2T ds 0
t1 t1 2T
t2
ie 2m T 2 ds 0
1
t1
t2 1
[2m( E V )] 2 ds 0 Since T +V = E
t1
t2 1
[2m( H V )] 2 ds 0 Since E = H
t1
t2 1
The above equation represents the principle of least action in term of arc length of the
particle trajectory.
(ii) Jacobi’s form of the principle of least action:
The K.E of the system
1
T a jk q j q k
2 j ,k
2T a jk q j q k
j ,k
d
2
= ....(7)
dt
2T dt 0
t1
d
t2
2T 1 0
t1 (2T ) 2
t2
2T 2 d 0
1
t1
t2
2E V 2 d 0
1
t1
t2
2H V 2 d 0
1
t1
t2
H V 2 d 0
1
....(8)
t1
2T dt 0
t1
t2
T dt 0
t1
If T is conserved, then
dt 0
t1
(t2 – t1) = 0
(t2 – t1) = extremum.
ie, the time taken by the light ray to travel between two points is extremum.
1.4 CANONICAL TRANSFORMATION:
Transformation is one to change one set of position and momentum coordinates into
another set of position and momentum coordinates.
We assume p j and q j are the old momentum and position coordinates and P j and Q j
are new one related by
P j =P j (p j, q j,t)
Q j = Q j (p j,q j,t) .…(1)
Then if there exists a new Hamiltonian H in the new coordinates such that
H H
Pj and Q j ….(2)
Q j Pj
and H p j Q j K
…(3)
Here the position and momenta coordinates are independent.
Canonical transformations are the transformations of phase space. They are
characterized by the property that they leave the form of Hamilton’s equations of motion.
example : Cartesian to polar coordinate is an example of co-ordinate transformation.
The transformation of one set of position coordinates to new set of coordinates is
called point transformation. They are the transformations of configuration space.
1.5 HAMILTONIAN FORMALISM:
Lagrangian equations of motion are invariant in form with respect to the set of any
generalized coordinates. In the new set Q j, Lagrange’s equations will be
d L L
0
dt Q
Q
j j
ie, Lagrange’s equations are covariant with respect to point transformations. If we define Pj
as
H
P j Q j , Pj
Q j
Therefore, this transformation is extended to Hamiltonian formulation. In
Hamiltonian formulation, we admit the existence of one more independent variable called
momentum. Consequently the simultaneous transformation of the independent coordinates
and momenta q j, p j to a new set Q j, P j can be represented in the form
Q j Q j q, p, t
….(1)
Pj Pj q, p, t
t2
Pj Q j k Q, P, t dt 0 ....(2)
t1
t2
p j q j Hq, p, t dt 0 ....(3)
t1
The simultaneous validity of equations (2) and (3) does not mean that the integrands
of the two integrals are equal. We can therefore write
p j q j H Pj Q j K dt 0
t2
....(4)
t1
Equation (4) will not be affected if we add to or subtract from it a total time derivative of a
function F = F(q,p,t).Now we can write equation (4) as
is follows that
p qj j
H Pj Q j K
dF
dt
....(5)
p j q j H Pj Q j K 1 q, Q, t
dF
....(6)
dt
F1 = F1 (q, Q, t)
dF1 F F F
1 q j 1 Q j 1
dt j q j j Q j t
F F F
1 p j q j Pj 1 Q j K H 1 0 .... (7)
q
Q j t
j
Since q j and Q j are to be treated as independent variables, equation (7) can hold only
separately vanish.
if q j and Q j
F1
ie pj q, Q, t ....(8)
q j
F1
Pj q, Q, t ....(9)
Q j
F1
and K = H q, Q, t ....(10)
t
p j q j H Pj Q j K
d
dt
F2 q j , Pj , t Pj Q j
F F F
Pj Q j K 2 q j 2 Pj 2 Pj Q j Pj Q j
q j Pj t
F2 F F
Pj Q j K 2 q j 2 Pj
t q j Pj
F F F
2 p j q j 2 Q j Pj H 2 K 0 ....(12)
q P t
j j
Since q j and P j are independent variable, equation (12) can be satisfied only when
F2
pj q j , Pj , t .....(13)
q j
F2
Qj q j , Pj , t .....(14)
Pj
F2
and K H
t
q j , Pj , t .....(15)
p j q j H Pj Q j K
d
dt
F3 p j , Q j , t p j q j
F F F
Pj Q j K 3 p j 3 Q j 3 p j q j p j q j
p j Q j t
F3 F F
H p j Q j K 3 q j p j 3 Pj Q j
t p Q
j j
Manonmaniam Sundaranar University, Directorate of Distance & Continuing Education, Tirunelveli.
14
F F F
3 q j p j 3 Pj Q j H 3 K 0 ....(17)
p Q t
j j
Since p j, Q j are independent variables, equation (17) can be satisfied only when
F3
qj p j , Qj , t ....(18)
p j
F3
Pj p j , Q j , t ....(19)
Q j
F3
k H p j ,Qj ,t ....(20)
t
equation (18) gives Q j = Q j (q, p, t)
and equation (19) gives P j = P j (q, p, t)
(D) FOURTH FORM F4(p, P, t)
F4 can be obtained from F3 by replacing Qj, by Pj
F4( p j, P j, t ) = F3 (p j, Q j, t) + P j Q j
= F1 (q j, Q j, t) -p j q j - +P j Q j
F1 (q j, Q j ,t) = F4(p j, P j, t) + p j q j - P j Q j ....(21)
putting equation (23) in equation (6) we get
p j q j H Pj Q j K
d
dt
F4 p j , Pj , t p j q j Pj Q j
F F F
PjQ j K 4 p j 4 Pj 4 p j q j p j q j PjQ j PjQ j
p j Pj t
F F F
H K 4 q j p j 4 Q j P j 4
p P t
j j .… (22)
Since p j and P j are independent variables, equation (22) can be satisfied only when
F4
qj p j , Pj , t
p j
.…(23)
F4
Qj p j , Pj , t )
Pj
….(24)
F4
K H p j , Pj , t
t ….(25)
If the expression (P j dQ j- p j dq j)
(or) (p j dq j- P j dQ j)
canonical.
Proof:
p j q j H Pj Q j K
dF
dt
….(1)
must be satisfied.
(or)
p j dq j Pj dQ j dF
Q 2q e cos p
P 2q e sin p
is a canonical transformation.
Solution:
Q 2q e cos p
PdQ pdq 2q e sin p 2q e cos p dq 2q e sin p 2q e sin pdp pdq
1/ 2 1/ 2 1/ 2 1/ 2
1
sin 2 p p dq 2q sin 2 pdp
= 2
1 1
q sin 2 p pq dq q sin 2 p pq dp
= q 2 p 2
F F
dq dp
= q p
= dF
Which shows that the RHS is an exact differential of the function
1
F q sin 2 p pq
2
and hence the transformation is canonical.
2. Show that the transformation
q 2 p sin Q p 2 p cos Q
(or)
P
2
1 2
p q2 Q tan 1
q
p
is canconical
Solution:
dq (2 p)1/ 2 cos Q dQ 2 p
1/ 2
sin Q dP
1 1
P sin 2Q dQ P sin 2Q dP
Q 2 P 2
F F
= dQ dP
Q P
1
= dF where F P sin 2Q
2
3. Show that the transformation
P = q cot p
sin p
Q = log is canonical. show that the generating function is
q
F e Q 1 q 2 e 2Q 1/ 2
q sin 1 (qe Q )
sin p sin p
log dq log dp
q q p q
1 sin p 1 cos p
2 dq dp
sin p q sin p q
q q
dq
cot p dp
q
dq
pdq PdQ pdq q cot p cot p dp
q
p cot p dq q cot 2 p dp
= qp q cot pdq qp q cot p dp
q p
F F
dq dp
q p
dF exact differential
F = qp +q cot p
and hence the transformation is canonical.
sin p
Let us put Q log
q
sin p q eQ p sin 1 qeQ
cos p = 1 q 2e 2Q 1/ 2
cot p =
1 q e
2 2Q 1/ 2
qe Q
Now the generating function F q sin 1 q eQ eQ 1 q 2e1 1/ 2
1.7 HAMILTONIAN:
The quantity p j q j L is a constant of motion with the condition that L does not
H as Hamiltonian H = H (q j, p j)
= p j q j L(q j , q j )
H = p j q j L
But Lagrangian L L q j , q j , t
L L L
dL dq j dq j dt ....(3)
q j q j t
Substituting equation (3) in equation (2) we get
L L L
dH q j dp j p j dq j dq j dq j dt
q j q j t
….(4)
L L
p j and p j
q j q j
H
p j
q j
….(6)
L H
t t
Equations (6) are known as Hamilton’s canonical equations of motion and of a set of 2n first
order equations of motion.
dH H
Therefore p j q j q j p j
dt t
H
= ….(1)
t
L
=
t
H L
.…(2)
t t
L
If L is not an explicit function of time, 0
t
H dH
0 ie 0
t dt
H = constant .…(3)
Thus if L is not an explicit function of time, H is constant of motion.
For conservative systems, the Potential energy does not depend upon generalized
velocity, ie
V
0 .…(4)
q j
We know H p jq j L
L
q j L
q j
T V
q j L
q j q j
T V
q j L since 0
q j q j
1 2
H q j m q j L
q j 2
m q j L
2
= 2T – L = 2T – (T –V)
= T + V = K.E + P.E = Total energy.
H represents the total energy of the system for conservative system.
1.8 CYCLIC (OR) IGNORABLE COORDINATE:
We know that the Lagrangian L is a function of generalized coordinate q j,
generalized velocity 𝑞 j and time t. If the Lagrangian of a system does not contain a
L
particular coordinate q k, then 0. such a coordinate is referred to as an ignorable or
qk
cyclic coordinate.
(i) Generalised momentum: [conjugate (or) canonical momentum]
Consider a system of mass points acted upon by forces derived from potentials
dependent on position only.
Now Lagrange’s equations of motion are
d L L
0
dt q j q j
Suppose q j is cyclic.ie it does not occur in Lagrangian L, then for this coordinate Lagrange’s
equation reduces to
d L
0
dt q j
L
constant
q j
Generalized momentum = constant.
L T V T
xi xi xi xi
Manonmaniam Sundaranar University, Directorate of Distance & Continuing Education, Tirunelveli.
21
1
ml xi y i zi
xi 2
2 2 2
mi xi pix
= x component of the linear momentum associated with the ith particle.
Generalising the concept, a momentum associated with the coordinate q j shall be
T L
pj = generalized momentum.
q j q j
Thus p j constant
L
1
2
m r 2 r 2 2
k
r
does not occur in L and is therefore an ignorable coordinate so that corresponding
momentum
L
p mr 2 l a constant
It appears from the expression L that we can solve Kepler problem without
considering which is ignorable. But this is not so in Lagrangian formation because which
requires how varies with t. Thus we consider . However Routhian procedure which
eliminates this consideration. We want to find a function R called Routhian function such
that it does not contain generalized velocities corresponding to ignorable coordinates.
L Lq1 , q 2 ,.....q n , q 1 , q 2 .....q n , t
If coordinates q1......q k are ignorable then
k
L n
L n
L L
(or) L q q j q j t . … (1)
j 1 q
j j j k 1 q j j k 1 q
j t
k k
R L q j p j q j p j
J 1 J 1
k k
L p jq j q j p j
J 1 J 1
k
L k
L q j q j p j
J 1 q
j J 1
n
L n
L L k
J k 1 q j
q j
J k 1 q
j
q j t q j p j
t J 1
….(3)
so that R = L p
L m r 2 2
=
1
2
m r 2 r 2 2
k
r
R
mr
r
R k
mr 2 2
r r
Now we get the equation of motion as
d
mr mr2 k2 0 [from equation (5)]
dt r
k
mr mr2 2 0
r
k
mr mr 2
r2
mv 2
f (r )
r
Which is the equation of motion of a particle under central force.
1.10 DERIVATION OF HAMILTON’S CANONICAL EQUATIONS FROM
VARIATIONAL PRINCIPLE:
Hamilton’s principle is stated as
t2
I Ldt 0
t1
H p j q j L
L p j q j H
t2
p q
t2
d H q j , p j , t dt 0
j j
t1
t
But 0 since time of travel along every path is same. Also
t2
q j t2
d q j
pj
t1
dt p j
t1
dt
dt
q j t2 t2
q j
pj | p j dt
t1 t1
t2
q j
p j dt
t1
q j
vanishes at limits t1 and t2,
Now equation (2) becomes
t2
p j q H q j H p j
I d q j p j j dt 0
t1
q j p j
q j
Putting d q j
p j
d p j
t2
H H
We get I p j q j p jq j q j p j dt 0
t1 q j p j
t2
H H
p j q j q j p j dt 0
t1
p j q j
Since p j and q j are independent variables, the above equation can be satisfied only if
H H
qj 0 q j
p j p j
H H
p j 0 p j
q j q j
which are the desired canonical equations of motion.
= p L
1
= ml 2 2 ml 2 2 mgl (1 cos )
2
1 2 2
= m l mgl (1 cos ) ….(3)
2
=T+V
Then the system is conservative.
Putting equation (2) into equation (3), we get
2
1 p
H ml 2 2 mgl (1 cos )
2 ml
giving
H p
2
p ml
….(4)
H
mgl sin
Thus Hamilton’s equations of motion for this system will be
H p
2 _....(5)
ml
H
p mgl sin ….(6)
From equation (5), we have
p ml 2
and hence equation (6) becomes
l g sin
l g sin 0
l g 0
g 0
l
which represents the equation of motion of a simple pendulum with period 2 l / g
2. Compound pendulum :
1 2
T I
2
V = mglcos
2
Now the Lagrangian L = I mgl cos
2
L
Then p
= I …. (1)
Hamiltonian H = p j q j L
= p L
1
I 2 I 2 mgl cos
2
1 2
= I mgl cos
2
2
1 p
= I mgl cos
2 I
2
p
= mgl cos ….(2)
2I
=T+V
Then the system is conservative.
The Hamilton’s equations of motion for and p are
H
P
….(3)
H
p
From equation (2) we find
H
mgl sin
p
Now equation (3) becomes ….(4)
I
and p mgl sin ….(5)
mgl
which is the equation of motion of compound pendulum with
I
3. Linear Harmonic Oscillator :
The system is conservative and constraint is independent of time. Hamiltonian will
represent the total energy of the system. The Lagrangian
L T V
1 1
mx 2 kx2
2 2
L T
p mx
x x
p
x
m
giving Hamiltonian
H T V
2
1 p 1
m kx2
2 m 2
p2 1
kx2
2m 2
and
H p r
r
p r m
H p
2
p mr …. (2)
1 vA
mv j q j j
2
…. (1)
2 c
L L
pj
q j v j
q
mv j Aj
c
Hamiltonian H q jp j L
c
q
mv 2 v A L
c
q 1
2T v A T q v A
c c
T q
1
mv 2 q
2
such that H can be interpreted as the sum of kinetic and electrostatic potential energies of the
particle.
T
1
2
m x 2 y 2 z 2
The applied force on the body is its weight acting in negative z direction.ie
V
F Fz mg
z
V mgz
Lagrangian
L T V
1
2
m x 2 y 2 z 2 mgz
L T
p x mx
x x
px
giving x .
m
py pz
Similarly y and z
m m
Hamiltonian for such a system is conserved. ie
H T V
1
2
m x 2 y 2 z 2 mg z
H p y
y
p y m
H p z
z ….(2)
p z m
From equation (2) we get
p x
x 0
m
p y
y 0
m
p z
z g
m
which shows that the acceleration along z direction is the acceleration due to gravity and is
true.
UNIT 2
Poisson bracket, Special cases of Poisson bracket , Poisson theorem, Poisson bracket
and canonical transformation, Jacobi identity and its derivation, Lagrange bracket and its
properties, the relationship between Poisson and Lagrange brackets and its derivation, the
angular momenta and Poisson bracket, Liouville’s theorem and its applications; Theory of
small oscillations:
Formulation of the problem, Eigenvalue equation and the principle axis transformation,
frequencies of free vibration and normal coordinates, free vibrations of a linear triatomic
molecule
32
CLASSICAL AND STATISTICAL MECHANICS
Also
= 1 if i = j
Proof :
X ∂X ∂(Y + Z) ∂X ∂(Y + Z)
[X, Y + Z] = −
j
∂qj ∂pj ∂pj ∂qj
X ∂X ∂Y X ∂X ∂Y
∂X ∂Z ∂X ∂Z
= + − +
j
∂q j ∂p j ∂q j ∂p j j
∂p j ∂q j ∂pj ∂qj
X ∂X ∂Y ∂X ∂Y
X
∂X ∂Z ∂X ∂Z
= − + −
j
∂qj ∂pj ∂pj ∂qj j
∂qj ∂pj ∂pj ∂qj
= [X, Y ] + [X, Z]
Similarly
X ∂qi ∂qj ∂qi ∂qj
[qi , qj ]q,p = −
k
∂qk ∂pk ∂pk ∂qk
∂qj ∂qi
= =0
∂pk ∂pk
and hence
[qi , qj ]q,p = 0 = [pi , pj ]q,p
since
∂qi ∂pj
and =0
∂pk ∂qk
33
CLASSICAL AND STATISTICAL MECHANICS
Also
∂qi ∂pj
and =0
∂qk ∂pk
sothat
X
[qi , pj ]q,p = δik δjk
= δij
= 0 if i 6= j
= 1 if i = j
Excercises
We have
X ∂φ ∂ψ ∂ψ ∂φ
[φ, ψ] = −
i
∂qi ∂pi ∂qi ∂pi
∂ X ∂ ∂φ ∂ψ ∂ψ ∂φ
[φ, ψ] = −
∂t i
∂t ∂q i ∂p i ∂qi ∂pi
X ∂ ∂φ ∂ψ X ∂φ ∂ ∂ψ
= +
i
∂q i ∂t ∂p i i
∂qi ∂pi ∂t
X ∂
∂ψ ∂φ X ∂ψ ∂ ∂φ
− −
i
∂qi ∂t ∂pi i
∂qi ∂pi ∂t
X ∂
∂φ ∂ψ ∂ψ ∂ ∂φ
= −
i
∂qi ∂t ∂pi ∂qi ∂pi ∂t
X ∂φ ∂ ∂ψ
∂
∂ψ ∂φ
+ −
i
∂qi ∂pi ∂t ∂qi ∂t ∂pi
dφ dψ
= , ψ + φ,
dt dt
34
CLASSICAL AND STATISTICAL MECHANICS
2.If {pl , qi }, {ql , qi } are the Lagrange’s brackets and [pl , pj ], [ql , pj ] are the Poisson brackets,
then prove that
n
X n
X
{pl , qi }[pl pj ] + {ql , qi }[ql , pj ] = 0
i=1 i=1
We have
{pl , qi } = −{qi , pl }
= −δil
and
{ql , qi } = 0
[pl , pj ] = 0
[ql , pj ] = δlj
=0
3.If {ql , qi }, {pl , qi } are the Lagrange’s brackets and [ql , qj ], [pl , qj ] are the Poisson brackets,
then prove that
n
X n
X
{ql , qi }[ql , qj ] + {pl , qi }[pl , qj ] = δij
i=1 i=1
We know
{ql , qi } = 0
[ql , qj ] = 0
Also
{pl , qi } = −{qi , pl }
= −δil
[pl , qj ] = −[qj , pl ]
= −δjl
35
CLASSICAL AND STATISTICAL MECHANICS
n
X n
X
{ql , qi }[ql , qj ] + {pl , qi }[pl , qj ] = 0 × 0 + [−δil × −δjl ]
i=1 i=1
= δij
Hence proved.
dF ∂F
= [F, H] +
dt ∂t
∂F
[F, H] + =0
∂t
[F, H] = 0
Proof :
X ∂X ∂Y ∂X ∂Y
[X, Y ]Q,P = −
i
∂Qi ∂Pi ∂Pi ∂Qi
X ∂X ∂Y ∂qj ∂Y ∂pj
∂X ∂Y ∂qj ∂Y ∂pj
= + − +
i,j
∂Q i ∂q j ∂P i ∂p j ∂P i ∂P i ∂q j ∂Q i ∂pj ∂Qi
X ∂Y X ∂X ∂qj
∂X ∂qj X ∂Y X ∂X ∂pj ∂X ∂pj
= − + −
j
∂qj i ∂Qi ∂Pi ∂Pi ∂Qi j
∂pj i ∂Qi ∂Pi ∂Pi ∂Qi
X ∂Y ∂Y
= [X, qj ]Q,P + [X, pj ]Q,P
j
∂q j ∂p j
36
CLASSICAL AND STATISTICAL MECHANICS
Further
Similarly
X ∂Y ∂X ∂Y ∂X
[X, Y ]Q,P = − +
j
∂q j ∂p j ∂pj ∂qj
= [X, Y ]q,p
37
CLASSICAL AND STATISTICAL MECHANICS
∂pj ∂ 2 F2 ∂Qi
= = (6)
∂Pi ∂Pi ∂qj ∂qj
∂qj ∂ 2 F3 ∂Pi
=− = (7)
∂Qi ∂Qi ∂pj ∂pj
and
∂qj ∂ 2 F4 ∂Qi
= =− (8)
∂Pi ∂Pi ∂pj ∂pj
and similarly
Excercises
q
tanQ =
p
1
P = (q 2 + p2 )
2
38
CLASSICAL AND STATISTICAL MECHANICS
∂Q 1 ∂Q −q
sec2 Q = sec2 Q = 2
∂q p ∂p p
∂P ∂P
= q =p
∂q ∂p
Then
∂Q ∂P ∂Q ∂P
[Q, P ] = −
∂q ∂p ∂p ∂q
1 2 q 2
= cos Q p + cos Q q
p p2
q2
2
= cos Q 1 +
fflp2
cos2 Q 1 + tan2 Q
=
= cos2 Q × sec2 Q = 1
= (e−2q − p2 )1/2
39
CLASSICAL AND STATISTICAL MECHANICS
and
Q
P = tan−1
p
Now
Q
tanP =
p
Q2
1 + tan2 P = 1 +
p2
Q2 + p2
1 + tan2 P =
p2
Q2 + p2
sec2 P =
p2
p
Q2 + p2
sec P =
p
p
cos P = p
Q2 + p2
p
=p
e−2q − p2 + p2
p
=√
e−2q
p
= −q
e
= p eq
P = cos−1 (p eq )
Q = (e−2q − p2 )1/2
P = cos−1 (p eq )
Then
40
CLASSICAL AND STATISTICAL MECHANICS
Therefore
∂Q ∂P ∂Q ∂P
[Q, P ] = −
∂q ∂p ∂p ∂q
−e−2q
−1
=
(e−2q − p2 )1/2 (e−2q − p2 )1/2
−p −p
−
(e−2q − p2 )1/2 (e−2q − p2 )1/2
e−2q p2
= −2q −
(e − p2 ) (e−2q − p2 )
(e−2q − p2 )
= −2q =1
(e − p2 )
Furthermore [Q.Q] = 0 and [P, P ] = 0
Hence the transformation is canonical.
Ḟ = [F, H] (1)
If the Poisson bracket of F with H vanishes then F =constant of motion. This requirement
does not however require that H should be a constant of motion. Suppose such dynamical
variables are qj and pj , then
q˙j = [qj , H]
and
p˙j = [pj , H] (2)
The above equations are identical with Hamilton’s canonical equations of motion
X ∂qj ∂H ∂qj ∂H
[qj , H] = −
i
∂H ∂p i ∂pi ∂qi
Since
∂qj
=0
∂pi
41
CLASSICAL AND STATISTICAL MECHANICS
we get
∂H
[qj , H] = δij
∂pi
∂H
= (f or i = j)
∂pj
Therefore
∂H
q˙j = = [qj , H]
∂pj
Similarly
∂H
pj = −
= [pj , H]
∂qj
Equations (2) can thus be known as equations of motion inpoisson bracket form.
If Poisson bracket [pj , H] vanishes, then
p˙j = 0
pj = constant
That is the linear momentum is conserved and hence the corresponding co-ordinate is cyclic.
Thus all functions whose Poisson bracket with Hamiltonian vanish will be constants of
motion and conversely Poisson brackets of all constants of motion with H must vanish.
We have
∂[X, Y ] ∂Z ∂Z ∂[X, Y ]
[X, [Y, Z]] = −
∂qj ∂pj ∂qj ∂pj
∂ ∂X ∂Y ∂Y ∂X ∂Z ∂Z ∂ ∂X ∂Y ∂Y ∂X
= − − −
∂qj ∂qj ∂pj ∂qj ∂pj ∂pj ∂qj ∂pj ∂qj ∂pj ∂qj ∂pj
2 2 2
∂Y ∂ 2 X ∂Z
∂ X ∂Y ∂X ∂ Y ∂ Y ∂X
= + − −
∂qj 2 ∂pj ∂qj ∂qj ∂pj ∂qj 2 ∂pj ∂qj ∂qj 2 ∂pj
2
∂ 2 Y ∂X ∂ 2 X ∂Y
∂Z ∂ X ∂Y ∂Y ∂X
− + − −
∂qj ∂pj ∂qj ∂pj ∂p2j ∂qj ∂pj ∂qj ∂pj ∂p2j ∂qj
∂ 2 X ∂Y ∂Z ∂ 2 X ∂X ∂Z ∂ 2 Y ∂X ∂Z ∂ 2 Y ∂X ∂Z
= + − +
∂qj 2 ∂pj ∂pj ∂p2j ∂qj ∂qj ∂qj 2 ∂pj ∂pj ∂pj 2 ∂qj ∂qj
∂ 2X ∂ 2Y
∂Y ∂Z ∂Y ∂Z ∂X ∂Z ∂X ∂Z
− + + +
∂pj ∂qj ∂pj ∂qj ∂qj ∂qj ∂pj ∂qj ∂qj ∂pj ∂pj ffl∂qj
42
CLASSICAL AND STATISTICAL MECHANICS
Similarly
∂ 2 Z ∂X ∂Y ∂ 2 Z ∂X
2
∂ 2 X ∂Z ∂Y
∂Y ∂ X ∂Z ∂Y
[[Z, X], Y ] = + − +
∂qj 2 ∂pj ∂pj ∂p2j ∂qj ∂qj ∂qj 2 ∂pj ∂pj ∂pj 2 ∂qj ∂qj
∂ 2Z ∂ 2X
∂X ∂Y ∂X ∂Y ∂Z ∂Y ∂Z ∂Y
− + + +
∂pj ∂qj ∂pj ∂qj ∂qj ∂pj ∂pj ∂qj ∂qj ∂pj ∂pj ffl∂qj
∂ 2 Y ∂Z ∂X ∂ 2 Y
2
∂ 2 Z ∂Y ∂X
∂Z ∂X ∂ Z ∂Y ∂X
[[Y, Z], X] = + − +
∂qj 2 ∂pj ∂pj ∂p2j ∂qj ∂qj ∂qj 2 ∂pj ∂pj ∂pj 2 ∂qj ∂qj
∂ 2Y ∂ 2Z
∂Z ∂X ∂Z ∂X ∂Y ∂X ∂Y ∂X
− + + +
∂pj ∂qj ∂pj ∂qj ∂qj ∂pj ∂pj ∂qj ∂qj ∂pj ∂pj ffl∂qj
Adding all we get
[X, [Y, Z]] + [Y, [Z, X]] + [Z, [X, Y ]] = 0
taken over an arbitrary two dimensional surface S of the 2n dimensional (q,p)phase space is
invariant under canonical transformation.
Position of a point on any two dimensional surface is expressed as
qj = qj (u, v)
pj = pj (u, v)
(3)
Transforming the integral (2) in terms of (u,v), we write
∂(qj , pj )
dqj dpj = dudv (4)
∂(u, v)
43
CLASSICAL AND STATISTICAL MECHANICS
with
∂qj ∂pj
∂u ∂u
∂(qj , pj )
= (5)
∂(u, v)
∂qj ∂pj
∂v ∂v
as the Jacobian.
Further
Z Z X Z Z X
dqj dpj = dQj dPj (6)
S j S j
where (Qj , Pj ) is another set of canonical co-ordinates to which the set(qj , pj )has been trans-
formed.
Now relation (6) becomes
Z Z X Z Z X
∂(qj , pj ) ∂(Qj , Pj )
dudv = dudv
S j ∂(u, v) S j ∂(u, v)
S is arbitrary and area dudv is arbitrary. Therefore expressions on boyh the sides will be
equal only when
X ∂(qj , pj ) X ∂(Qj , Pj )
=
j
∂(u, v) j
∂(u, v)
∂qj ∂pj ∂Qj ∂Pj
X ∂u ∂u
∂u ∂u
X
=
j ∂q ∂pj j ∂Q ∂Pj
∂vj ∂v
∂vj ∂v
or
X ∂qj ∂pj ∂pj ∂qj
X ∂Qj ∂Pj ∂Pj ∂Qj
− = −
j
∂u ∂v ∂u ∂v j
∂u ∂v ∂u ∂v
44
CLASSICAL AND STATISTICAL MECHANICS
X ∂qj ∂pj
∂pj ∂qj
{u, v} = −
j
∂u ∂v ∂u ∂v
X ∂pj ∂qj ∂qj ∂pj
=− −
j
∂u ∂v ∂u ∂v
X ∂qj ∂pj ∂pj ∂qj
=− −
j
∂v ∂u ∂v ∂u
= {v, u} (8)
{qi , qj } = 0
{pi , pj } = 0
{qi , qj } = δij
X ∂qk ∂pk ∂qk ∂pk
{qi , qj } = − (9)
k
∂qi ∂qj ∂qj ∂qi
q’s and p’s are independent and hence
∂pk ∂pk
= 0 and =0
∂qj ∂qi
Then
{qi , qj } = 0
{pi , pj } = 0 (10)
Now
X ∂qk ∂pk ∂qk ∂qk
{qi , pj } = −
k
∂qi ∂pj ∂pj ∂qi
q’s and p’s are independent and hence
∂qk
=0
∂pj
45
CLASSICAL AND STATISTICAL MECHANICS
Then
X ∂qk ∂pk
{qi , pj } = (11)
k
∂qi ∂pj
We have
∂qk ∂pk
= δki and = δkj
∂qi ∂pj
Hence
X
{qi , pj } = δki δkj
k
= δij (12)
(1)
X ∂pk ∂uj X 2n
∂qk ∂ul X ∂pk ∂uj ∂qk
· = ·
k,m=1n
∂u i ∂pm
l=1
∂u l ∂qm
k,m
∂u i ∂pm ∂qm
X ∂pk ∂uj
= · δkm
k,m
∂ui ∂pm
But
∂pm
δkm =
∂pk
The first four terms on R.H.S on multiplication is
X ∂pk ∂uj X 2n
∂qk ∂ul X ∂pk ∂uj ∂pm
· = ·
k,m=1n
∂ui ∂pm l=1 ∂ul ∂qm k,m
∂ui ∂pm ∂pk
46
CLASSICAL AND STATISTICAL MECHANICS
X ∂pk ∂uj
=
k,m
∂ui ∂pk
(2)
The last four terms will be
X ∂qk ∂uj X 2n
∂pk ∂ul X ∂qk ∂uj ∂pk
· = ·
k,m=1n
∂u i ∂qm
l=1
∂u l ∂pm
k,m
∂u i ∂qm ∂pm
X ∂qk ∂uj
= · δkm
k,m
∂ui ∂qm
X ∂qk ∂uj ∂qm
= ·
k,m
∂ui ∂qm ∂qk
X ∂qk ∂uj
=
k,m
∂ui ∂qk
(3)
The second term is
X ∂pk ∂uj X 2n
∂qk ∂ul
− · =0
k,m=1n
∂u i ∂qm
l=1
∂u l ∂pm
since
2n
X ∂qk ∂ul ∂qk
= =0
l=1
∂ul ∂pm ∂pm
Similarly the third term will be zero.
Hence R.H.S of equation (1) is
X ∂pk ∂uj X ∂qk ∂uj X ∂uj ∂pk ∂uj ∂qk
+ = +
k
∂ui ∂pk k
∂ui ∂qk k
∂pk ∂ui ∂qk ∂ui
∂uj
= (qk , pk )
∂ui
∂uj
= = δij
∂ui
(4)
Therefore
2n
X
{ul , ui } [ul , uj ] = δij
i=1
47
CLASSICAL AND STATISTICAL MECHANICS
The angular momentum l can be expressed in terms of linear momentum p and radius
vector r as
l=r×p
giving
lx = ypz − zpy ; ly = zpx − xpz ; lz = xpy − ypx
Therefore
48
CLASSICAL AND STATISTICAL MECHANICS
= (xpx − ypx ) = lz
and
[lz , lx ] = ly
dρ
i.e =0
dt
d
i.e (δΓ) = 0
dt
Consider any fixed element of volume of phase space located between q1 and q1 + δq1 , ...qf
and qf + δqf , p1 and p1 + δp1 , ...pf and pf + δpf .
If ρ is the density of phase points, the number of phase points in this volume element at any
instant t is
δN = ρ.δΓ = ρ(δq1 ...δqf δp1 ...δpf )
49
CLASSICAL AND STATISTICAL MECHANICS
q´1
[ ( ) ]
q´1 +
q´1
q1
δ q1
p
D C
p1 + δ p1
δ p1
p1
A B
δ q1
q1 q1 + δ q1
q
The number of phase points located in the volume (δq1 ...δqf δp1 ...δpf ) changes as the co-
ordinates and momenta vary. The change in the number of phase points within this volume
of phase space in time dt is
d(δN ) d(ρ.δΓ) dρ
= = dt(δq1 ...δqf δp1 ...δpf )
dt dt dt
This change is due to the number of phase points entering and leaving this volume in time
dt.
Consider two faces of hypervolume normal to the q-axis with coordinates q1 and q1 + δq1 .
Number of phase points entering the first phase in time dt is
ρ and q˙1 are the density and velocity component at (q1 ...qf ; p1 ...pf )
Number of phase points entering the second phase in time dt is
∂ρ ∂ q˙1
ρ+ δq1 q˙1 + δq1 dtδq2 ...δqf δp1 ...δpf (2)
∂q1 ∂q1
Eqn(1)-Eqn(3) gives
∂ q˙1 ∂ρ
− ρ + q˙1 dtδq1 ...δqf δp1 ...δpf (4)
∂q1 δq1
50
CLASSICAL AND STATISTICAL MECHANICS
∂ p˙1 ∂ρ
− ρ + p˙1 dtδp1 ...δqf δp1 ...δpf (5)
∂p1 δp1
The net increase in the number of phase points in time dt is in this volume of phase space
is obtained by summing the net number of phase points entering the volume through all the
faces labelled by q1 ...qf and p1 ...pf . Hence
f
d X ∂ q˙j ∂ p˙j ∂ρ ∂ρ
(δN ) = − ρ + + q˙j + p˙j dtδq1 ...δqf δp1 ...δpf (6)
dt j=1
∂qj ∂pj ∂qj ∂pj
f
∂ρ X ∂ q˙j ∂ p˙j ∂ρ ∂ρ
= − ρ + + q˙j + p˙j (7)
∂t j=1
∂q j ∂p j ∂q j ∂p j
f
X ∂ q˙j ∂ p˙j
+ =0 (8)
j=1
∂qj ∂pj
51
CLASSICAL AND STATISTICAL MECHANICS
f
∂ρ X ∂ρ ∂ρ
=− q˙j + p˙j (9)
∂t q,p j=1
∂qj ∂pj
f
∂ρ X ∂ρ ∂ρ
+ q˙j + p˙j =0 (10)
∂t q,p j=1
∂qj ∂pj
This result is known as Liouville’s theorem.This equation is identical with the equation of
continuity in hydrodynamics.
If ρ is a function of q, p and t and q, p are functions of t, then
dρ ∂ρ ∂ρ dq ∂ρ dp
= + +
dt ∂t ∂q dt ∂p dt
On generalization we get
f f
dρ ∂ρ X ∂ρ dqj X ∂ρ dpj
= + + (11)
dt ∂t j
∂q j dt j
∂pj dt
dρ
=0
dt
This form is called the principle of the conservation of density in phase space.
Now
d dρ d
(δN ) = δΓ + ρ (δΓ) (12)
dt dt dt
d
(δN ) = 0
dt
52
CLASSICAL AND STATISTICAL MECHANICS
dρ d
δΓ + ρ (δΓ) = 0 (13)
dt dt
Already we have proved that
dρ
=0
dt
It follows that
d
ρ (δΓ) = 0 (14)
dt
But ρ 6= 0, we get
d
(δΓ) = 0 (15)
dt
This equation gives the principle of conservation of extension in phase space.
Excercises :
then show that (i) the transformation is canonical and (ii)the generating function of this
transformation is
F3 = −(eQ − 1)2 tan p
Solution :
(i) For the trasformation to be canonical (pdq − P dQ) must be an exact differential.
53
CLASSICAL AND STATISTICAL MECHANICS
∂ 1 ∂ 1
= [q(p − sin 2p)]dq + [q(p − sin 2p)]dp
∂q 2 ∂p 2
1
= d[q(p − sin 2p)]
2
= an exact dif f erential
(ii) We have
eQ = (1 + q 1/2 cos p)
eQ − 1 = q 1/2 cos p
eQ − 1
q 1/2 =
cos p
(eQ − 1)2
q= (1)
(cos p)2
Now
P = 2(1 + q 1/2 cos p)q 1/2 sin p
(eQ − 1)
= 2eQ sin p
cos p
= 2eQ (eQ − 1)tan p (2)
We know that
∂F3
q=−
∂p
∂F3
P =−
∂Q
Now
∂F3 (eQ − 1)2
=− = −(eQ − 1)2 sec2 p (3)
∂p cos2 p
∂F3
= −2eQ (eQ − 1)tan p (4)
∂Q
Integrating eqn(3) we get
Z
F3 = − (eQ − 1)2 sec2 p dp + constant
54
CLASSICAL AND STATISTICAL MECHANICS
Since
Equation(5) = Equation(6)
1
p= , p = P Q2
Q
is
q
F =
Q
Solution :
1
p= , p = P Q2
Q
We find
1
pdq = (2P QdQ + Q2 dP )
Q
= 2P dQ + QdP
Now
pdq − P dQ = 2P dQ + QdP − P dQ
= P dQ + QdP
∂ ∂
= (P Q)dQ + (P Q)dP
∂Q ∂P
∂F ∂F
= dQ + dP
∂Q ∂P
55
CLASSICAL AND STATISTICAL MECHANICS
= dF
q
P =
Q2
we have
q
F = Q
Q2
q
=
Q
Hence the result.
Consider a system with n degrees of freedom and whose total energy is conserved.Suppose
the system is in an equilibrium stste at the point having the co-ordinates q01 , q02 ...q0n . Thus
the potential energy satisfies the condition
∂V
=0
∂qj q=q0
Here q = q0 represents qj = q0j for all values of j. let ηj denotes a small displacement in the
corresponding co-ordinate qj from the equilibrim position. Sothat
qj = q0j + ηj
with j=1,2,...n
Now the potential energy becomes
56
CLASSICAL AND STATISTICAL MECHANICS
The departure from the equilibrium position is small and we truncate the series after the
first non-vanishing term in ηj .
The first term on the R.H.S is constant, the second term is zero and we consider upto the
third term.
The first term gives the shift the reference of potential energy to the equilibrium point.
Thus the equation (1) becomes
n n
∂ 2V
1 XX
V (q1 , q2 , ...qn ) = ηj ηk
2 j=1 k=1 ∂qj ∂qk q=q0
n n
1 XX
= Vjk ηj ηk (2)
2 j=1 k=1
Where
∂ 2V
Vjk =
∂qj ∂qk q=q0
Here Vjk = Vkj . For the stable equilibrium, the potential energy should be minimum at the
equilibrium position and hence Vjk must be positive. The kinetic energy of the system is
n n n n
1 XX 1 XX
T = cjk q˙j q˙k = cjk η˙j η˙k
2 j=1 k=1 2 j=1 k=1
cjk be the functions of the co-ordinates qj ’s and we can expand in Taylor series form about
the equilibrium position :
T has the terms of second order in η’s and we retain only the zeroth order term. Thus we
have
Here Tjk = Tkj .Using eqns (2) and (3) the Lagrangian of the system is
n n n n
1 XX 1 XX
L = Tjk η˙j η˙k − Vjk ηj ηk
2 j=1 k=1 2 j=1 k=1
57
CLASSICAL AND STATISTICAL MECHANICS
or
T η̈ + V η = 0
Here T and V represent the square matrices of order n × n.Both of these matrices are
real as well as symmetric.Each of the n equations (4) involves all the n variables and can
be simplified by transforming them into another set of n equations each of which involves
only one variable. This is possible with the help of the normal co-ordinates and normal
frequencies.Let us try the oscillatory solutions of the form
ηk = Ak eiωt (5)
here ω is the frequency of oscillations. The real part Ak cos ωt corresponds to the actual
motion. Using eqn (5) in (4) we get
n
X
(Vjk − ω 2 Tjk )Ak = 0 (6)
k=1
58
CLASSICAL AND STATISTICAL MECHANICS
(7)
This equation is known as the secular or characteristic equation.
Equation (7) may also be wrtten as
(V A − ω 2 T A) = 0
Each mode associated with an eigen frequency is known as the principal or natural
mode.Each of the normal frequencies must be real.The eigen functions corresponding to
different eigen values must be orthogonal to each other. Thus we have
n
X
Akl Akl0 = δll0
k=1
59
CLASSICAL AND STATISTICAL MECHANICS
We can choose a certain system of body axes with respect to which the off-diagonal
elements should disappear and only the diagonal elements remain in the expression for I.
Such axes are called principal axes of transformation. The corresponding moments of inertia
are called principal moments of inertia.
if we denote this form of inertia tensor by I0 and I1 , I2 , I3 stand for the principal values,
I 0 0
1
0
I = 0 I2 0 (1)
0 0 I3
Lx = I1 ωx
Ly = I2 ωy
Lz = I3 ωx (2)
That is each of angular momentum component along a principal axis is a function of cor-
responding angular velocity component only related to it through the principal moment of
inertia about that direction.
Thus in general if a rigid body is rotating about a principal axis, the angular momentum L
and ω are directed along any of the principal axes and therefore
L=Iω
60
CLASSICAL AND STATISTICAL MECHANICS
where I is the scalar, the moment of inertia about this axis.The angular momentum L and
angular velocity ω are along the principal axes and hence
L = Lxbi + Lybj + Lz k
b = I(ωxbi + ωybj + ωz k)
b
Thus
Lx = Iωx , Ly = Iωy , Lz = Iωz (3)
or
Ixx − I Ixy Ixz ωx
ωy = 0 (4)
Iyx Iyy − I Iyz
Izx Izy Izz − I ωz
For these equations to have non-trivial solutions, the determinant of the coefficients must
vanish. That is
Ixx − I Ixy Ixz
Iyx Iyy − I Iyz = 0 (5)
Izx Izy Izz − I
This is called secular equation of inertia and its solutions the secular values
or eigenvalues. We solve the determinantal equation which will be cubic in I and therefore
will furnish three values for I through I1 , I2 , I3 which are desired principal moments of inertia
:
The direction of any one principal axis is determined by substituting for I = I1 and determine
the ratios for ωx : ωy : ωz as
ωy ωz
= λ1 , = λ2
ωx ωx
Thus
ω = ωxbi + ωybj + ωz k
b
61
CLASSICAL AND STATISTICAL MECHANICS
and
q
2
ω
b = (bi + λ1bj + λ2 k)/
b 1 + λ21 + λ2
Hence we can determine the direction of ω or the direction of principal axis corresponding
to I1 .Similarly we may find the direction of the corresponding pricipal axis if we substitute
I2 or I13 .
Example:
If the symmetry axis of the body is taken as axis of rotation and the origin of body axes
lies on this then the principal axes are the symmetry axis and any two perpendicular axes
normal to the symmetry axis. In the case of a sphere, every axis through the centre is
symmetry axis and hence any three orthogoal axes through the centre are principal axes.
The co-ordinates in the solution of equations where only one single frequency is involved
in the solution are known as the normal co-ordinates.Thus the normal co-ordinates are
defined as the generalized co-ordinates wher each of them execute oscillations with a single
frequency.On transformation from the co-ordinates uj into the normal co-ordinates denotted
by ηk the lagrangian as well as equations of motion are changed.
Let the new co-ordinates ηk be related to the original co-ordinates as
X
uj = ajk ηk
k
u = aη (1)
Here u and η are column matrices of order n × 1 and a is a square matrix of order n × n.We
shall express potential and kinetic energy in terms of η.The potential energy is
1X
V = Vjk uj uk
2 j,k
1X
= uj Vjk uk
2 j,k
62
CLASSICAL AND STATISTICAL MECHANICS
1
V = uT V u
2
= diag(ω1 2 , ω2 2 ...ωn 2 )
1X
V = λl ηl 2
2 l
1X 2 2
V = ωl ηl (2)
2 l
The kinetic energy is given by
1 XX
T = Tjk u˙j u˙k
2 j k
1 XX
= u˙j Tjk u˙k
2 j k
1
= u̇T T u̇
2
1
= (aη̇)T T (aη̇)
2
1
= η̇ T aT T aη̇
2
1
= η̇ T η̇
2
Where
aT T a = 1
1X 2
= η̇l (3)
2
63
CLASSICAL AND STATISTICAL MECHANICS
Giving
n n
∂L X ∂L X
= η̇l and =− ωl 2 η̇l
∂ η̇l i=1
∂η l i=1
Now
d ∂L ∂L
− =0
dt ∂ η̇l ∂ηl
gives
n
X
(η̈l + ωl 2 η̇l ) = 0
i=1
η¨1 + ω1 2 η˙1 = 0
η¨2 + ω2 2 η˙2 = 0
η̈l + ωl 2 η̇l = 0
(4)
Thus each co-ordinate executes only one single frequency oscillation and therefore η1 , η2 etc
are termed as normal co-ordinates.
The solution of equation
η̈l + ωl 2 η̇l = 0
is
ηl = Al cos ωl t + Bl sin ωl t if ωl 2 > 0
ηl = Al t + Bl t if ωl 2 = 0
(5)
For ωl 2 > 0, all co-ordinates remain finite and the equilibrium is stable.
But for ωl 2 = 0 and ωl 2 < 0, the co-ordinates become infinite and the equilibrium is unstable.
64
CLASSICAL AND STATISTICAL MECHANICS
Since
ω = 2πν
(6)
Here A’s and B’s are 2n arbitrary constants determined by the initial conditions.
Suppose all constants except A1 and B1 are zero. Then only η1 will vary periodically
with time. This situation corresponds to a normal mode of vibration. That is the system
is vibrating in a normal mode.Therefore there will be n normal modes of vibration and n
normal frequencies ν1 , ν2 , ν3 ...νn corresponding to each normal co-ordinate η1 , η2 , η3 ...ηn Now
eqn(6) will take the form
ηl = Al cos (ωl t + δl )
1 1
T = m(q˙1 2 + q˙3 2 ) + M q˙2 2
2 2
65
CLASSICAL AND STATISTICAL MECHANICS
m 0 0 q˙
h i 1
2T = q˙1 q˙2 q˙3 0 M 0 q˙2
0 0 m q˙3
giving
m 0 0
T = (Tij ) = 0 M 0
0 0 m
1 1
V = k(q2 − q1 )2 + k(q3 − q2 )2
2 2
1
V = k(q1 2 + 2q2 2 + q3 2 − 2q1 q2 − 2q2 q3 )
2
or
k −k 0 q
h i 1
2V = q1 q2 q3 −k 2 k −k q2
0 −k k q3
66
CLASSICAL AND STATISTICAL MECHANICS
That is
k − ω 2 m
−k 0
2 −k = 0
−k 2k−ω M
k − ω 2 m
0 −k
After expansion we get
(k − ω 2 m)[(2k − ω 2 M )(k − ω 2 m) − k 2 ] − k 2 (k − ω 2 m) = 0
(k − ω 2 m)[(2k − ω 2 M )(k − ω 2 m) − 2k 2 ] = 0
giving
ω1 = 0
r
k
ω2 =
m
s
1 2
ω3 = k +
m M
Now we want to express η1 , η2 and η3 in terms of the generalized co-ordinates q1 , q2 and q3
as
X
qj = Ajk ηk
k
giving
q1 A11 A12 A13 η
1
q2 = A21 A22 A23 η2 (1)
q3 A31 A32 A33 η3
Our problem is to find the components of eigen vectors A1 , A2 and A3 . For this purpose we
apply the relation
3
X
(Vij − ω 2 Tij )Aij = 0 i = 1, 2, 3
j=1
First case :
ω = ω1 = 0 and calculation of the components A11 , A21 and A31 of A1 . We have
k −k 0 A
11
−k 2k −k A21 = 0
0 −k k A31
67
CLASSICAL AND STATISTICAL MECHANICS
giving
kA11 − kA21 = 0
−kA21 − kA31 = 0
It gives
A11 = A21 = A31 = α
Second case :
q
k
ω = ω2 = m
and calculation of the components A12 , A22 and A32 of A2 . We have
0 −k 0 A
12
kM
−k A22 = 0
−k 2k − m
0 −k 0 A32
giving
−kA22 = 0
kM
−kA12 + (2k − )A22 − kA32 = 0
m
−kA21 − kA31 = 0
It gives
A22 = 0 A12 = −A32 = β
Third case :
q
1 2
ω = ω3 = k m
+ M
and calculation of the components A13 , A23 and A33 of A3 . Then
2 2 1 2 k 2m
ω = ω3 = k + = (1 + )
m M m M
We have
− 2mk
M −km 0 A13
− kM −k A23 = 0
−k m
0 −k − 2mk
M A33
68
CLASSICAL AND STATISTICAL MECHANICS
giving
2mk
− A13 − kA23 = 0
M
kM
−kA13 − A23 − kA33 = 0
m
2mk
−kA23 − A33 = 0
M
It leads to
2mk
A13 = −kA23
M
2m 2m
A13 = −A23 ⇒ A23 = − A13
M M
Take
A13 = γ
Then
2m
A23 = − γ
M
Similarly
2mk
−kA23 − A33 = 0
M
That is
2mk
kA23 + A33 = 0
M
Which leads to
2mk 2mk 2m
A33 = −kA23 ⇒ A33 = −k(− γ)
M M M
A33 = γ
It gives
2m
A13 = A33 = γ A23 = − γ
M
Thus the matrix A is
α β γ
A = Aij = α 0 − 2mγ
M
α −β γ
AT T A = I
69
CLASSICAL AND STATISTICAL MECHANICS
α α α m 0 0 α β γ 1 0 0
2mγ =
β 0 −β 0 M 0 α 0 − M 0 1 0
γ − 2mγ
M γ 0 0 m α −β γ 0 0 1
α α α αm βm γm 1 0 0
=
β 0 −β αM 0 −2mγ 0 1 0
γ − 2mγ
M γ αm −βm γm 0 0 1
α2 (2m + M) 0 0 1 0 0
2β 2 m = 0 1 0
0 0
0 0 2γ 2 m(1 + 2m/M ) 0 0 1
It gives
1
α2 (2m + M ) = 1 ⇒ α = √
2m + M
1
2β 2 m = 1 ⇒ β = √
2m
1
2γ 2 m(1 + 2m/M ) = 1 ⇒ γ = p
2m(1 + 2m/M )
Thus the matrix A is
√ 1 √1 √ 1
2m+M 2m 2m(1+2m/M )
1 1
A = Aij =
√2m+M 0 √ (2)
2m(1+2m/M )
√ 1 √1 √ 1
2m+M 2m 2m(1+2m/M )
The normal co-ordinates can be obtained by using equation (2) in equation (1).
70
Unit III - TWO - BODY CENTRAL FORCE PROBLEM AND H - J THEORY
Two body central force problem: Reduction to the equivalent one body problem, the
equation of motion and first integrals, classification of orbits, the virial theorem, the
differential equation for the orbit, integral power law in time in the Kelper’s problem,
scattering in central force field;
H-J Theory: H-J equation and their solutions, use of H-J method for the solution of
harmonic oscillator problem, Hamilton’s principle function, Hamilton’s characteristic
function and their properties, Action angle variable for completely separable systems, the
Kelper’s problem in action angle variables
=
m1 m2 r2 m1r
m1 m 2
m1 r
= r2
m1 m 2
m1r
r2 R ….(5)
m1 m2
m2 r
Therefore r1 R and
m1 m2
m1r
r2 R ....(6)
m1 m2
2 2
1 m 2 r 1 m1 r
Hence L = m1 R m 2 R V( r )
2 m1 m 2 2 m1 m 2
2
1 2 m 2 r 2 2Rm 2 r
= m1 R
2
m 1 m 2 2
m 1 m 2
2
1 2 m1 r 2 2Rm1r
m 2 R V( r )
2
m 1 m 2 2
m1 m2
1 2 1 2
L = M R r V( r ) ....(8)
2 2
d L L
0 ….(9)
dt R R
d L L
0 ….(10)
dt r r
L L
M R and 0
R R
L V
r
r r
Hence equations (9) and (10) become
d
dt
M R 0 ....(11)
dt
d V
r
r
0 ....(12)
Equation (11) giving M R constant
(or) R constant
That is velocity with which the centre of mass moves is constant.
V
Equation (12) giving r = f (r)
r
representing equation of motion for the system under consideration. Consequently
we can ignore the first term in equation (8) and write
1
L = r 2 V(r ) ....(13)
2
which is effective in describing the motion of the components of r. But L is the
same as a single particle of mass µ moving at a distance r from a fixed centre of force
which gives rise to the potential energy V(r). Thus two body problem can be reduced to
the equivalent one body problem.
3.2 EQUATIONS OF MOTION AND FIRST INTEGRALS:
Let us describe the position of the particle in the plane poloar co-ordinates r and
1
T = m r 2 r 2 2
2
Manonmaniam Sundaranar University, Directorate of Distance & Continuing Education, Tirunelveli.
73
V = V (r)
The Lagrangian L = T –V
=
1
2
m r 2 r 2 2 V(r )
d 1 2
r 0
dt 2
1 2
r constant .…(2)
2
ie, Areal velocity = constant
The rate at which the area swept out by the radius vector is constant which is
Kepler‟s second law of planetory motion.
Suppose dA is the area swept out by the radius vector in time dt.
r rd
1
Then dA
2
1 2
= r d
2
dA 1 2 d 1 2
r r
dt 2 dt 2
Therefore from equation (2), we write
1 2
r = constant (or)
2
dA
= constant .…(3)
dt
And
d
dt
mr 2 0 .…(7)
l2 dV
m r 3
mr dr
d l2
m r
dr 2mr 2
V
d l
mr r V r
dr 2mr 2
d 1 2 d l 2 dr
(or) mr V
dt 2 dr 2mr 2 dt
d l2
V
dt 2mr 2
d 1 2 l2
mr
V ( r ) = 0
dt 2 2mr 2
1 2 l2
(or) mr V (r ) = constant
2 2mr 2
putting l mr 2 , then
1 2 l2 1 2 m2r 4 2
mr V (r ) mr V (r )
2 2mr 2 2 2mr 2
1 1
mr 2 mr 2 2 V (r )
2 2
1
m r 2 r 2 2 V (r )
2
TV E
Thus we write
1 2 l2
mr
V (r ) E
2 2mr 2 ….(8)
E is a constant of motion. This is another first integral of motion.
1
2 l 2
2
dr
r E V
2
m 2mr dt
dr
dt 1
2 l 2
2
E V
m 2mr 2
On integration we get
r
dr
t 1 .…(9)
2 l
2
r0 2
E V
m 2mr 2
where r ro at t 0.
This equation gives r as a function of time. ie r (t).
l
We have
mr 2
Thus integration of equations (9) and (10) provides us r (t) and (t). Then we can
locate the position of the particle on the path at any time t and the solution.
3.3 CLASSIFICATION OF ORBITS:
We have already derived
d l2
m r V
dr 2m r 2
dV l2
=
dr m r 3
l2
= f (r ) 3
mr
= f ′ = effective force.
f= f (r )
mr
2 2
mr 3
= f (r ) m r 2
= f (r )
m r
2
r
2
m v
= f (r )
r
= Central force + Centrifugal force.
Thus equation takes the form
m r f
It is an equation of motion for a particle subjected to actual force f (r) and pseudo
l2 dV
centrifugal force 3
. Also f
mr dr
V = effective potential energy
l2
V = f dr f r 3 dr
mr
l2
= f (r ) dr
2mr 2
Manonmaniam Sundaranar University, Directorate of Distance & Continuing Education, Tirunelveli.
77
l2
= V (r )
2mr 2
The second term is the potential energy related to the centrifugal force.
(i) Orbits under inverse square law of force :
The effective potential energy
k l2
V (r ) dr
r2 2mr 2
k l2
=
r 2mr 2
as a function of r has been plotted.
(ii) Motion with different values of k:
corresponds to energy E1
The radial motion is confined to the values r = r1= rmin and r = r2= rmax
1 l2
E1 mr 2 V (r )
2 2mr 2
1 2
E1 V (r ) mr T
2
The particle will not be at rest at these points. The motion is confined between the
areas of two circles of radii r1 and r2.
A possible shape for an attractive inverse square law of force is an ellipse with the
focus at the centre.
When r varies from r1 to r2 and back, the radius vector turns through an angle obtained
by
l
dt 0
mr 2
dr
dt 1
2 l 2 2
E v
m 2mr 2
l
2 dr
r
1
2 l 2 2
m E V
m 2mr 2
l
2 dr
r o
1
l 2 2
2m( E V ) 2
r
2m
When , the path is a closed orbit. When is not a rational fraction of 2π, the
n
path has the shape of a rosette. Such an orbital motion is processing motion.
dV d k l2
dr dr r 2mr 2
Manonmaniam Sundaranar University, Directorate of Distance & Continuing Education, Tirunelveli.
79
k l2
=
r 2 mr 3
Putting r = ro for which V(r) is minimum.
k l2 k l2
2
3
o 2
3
ro mro ro mro
l2
ro
mk
E4 = V (r ) |min
k l2
=
ro 2mro 2
mk 2 1 mk 2 1 mk 2
=
l2 2 l2 2 l2
In this case two bounds coincide
r3 = r2 = ro
l2
Orbit is a circle of radius .
mk
: r will be imaginary and therefore no physically meaningful motion is
Case (iv) : E< Vmin
possible.
(iii) Stability of orbits and conditions for closure:
(a) Stability:
The condition of stability in radial motion is given by the existence of local
minimum in V(r), the effective potential. That is we require that
2V (r )
0 at r = ro
r 2
2V (r )
given by 0
r
For any central force, the potential energy function is given by,
V (r ) krn 1 (for bounded motion)
l2
V (r ) k r n 1
2mr 2
V (r ) l2
k (n 1)r
n
r mr
3
r r r0 mr0
3
l2
k (n 1)r0
n
3
mr0
( 3 n)
l 2r0
k (n 1)
m
2V (r ) n 1 3l 2
kn( n 1) r
r 2 r r
0 4
0
mr0
l 2 (3 n ) n 1 3l 2 4
n ro ro ro
m m
l 2 4 3l 2 4
n r0 r0
m m
l2 4
(n 3)r0
m
Which is positive if n >3
Therefore any circular orbit with r = r0 under any central force can satisfy the
stability condition if n >3.
(b) Conditions for closure:
(iii) All bound orbits are closed only if for the inverse square law of force of
electrostatic attraction or gravitational type and for Hooke‟s type linear law of
force.
(iv) The condition for bound motion is that there is a bounded domain of r in
which V(r) E. The condition for stability of circular orbits is n >3 where f
(r) r n. That is for n = 1 and n = 2 only closed orbits exist.
3.4 VIRIAL THEROEM:
We consider a system of mass points with position vectors ri and applied forces Fi .
p i Fi
….(1)
Let us write a quantity
G pi . ri
dG
pi .ri p i .ri
dt ….(2)
= mi vi
2
= 2T
and p i . ri Fi . ri
G
So that 2T Fi . ri
t
The time average over a time interval will be
z
1 dG
0 dt
dt 2T Fi . ri
dG
(or) 2T Fi . ri
dt
1
G( ) G(0) ....(3)
If the motion is periodic and is its period then G() = G(0) and hence the RHS of
equation (3) is zero. If the motion is not periods but there is an upper bound to G and the
co-ordinates and velocities remain finite, then sufficiently large, RHS again be
approximated to zero. Hence
2T Fi . ri 0
1
T Fi . ri .…(4)
2
Equation (4) is known as the virial theorem and RHS is called virial of claussius.
This theorem is used in kinetic theory of gases to prove Boyle‟s law and to obtain
equation of state for imperfect gases etc.
It the forces are derivable from a potential, then the theorem becomes.
1
T iV ri
2
1 V
rˆ . r
2 r i
1 V
= r
2 r
If the potential energy is
V k r n 1
l d du
And r
m dt d
l d du d
m d d dt
l d 2u
m d 2
l 2u 2 d 2 u
m 2 d 2
We know that the equation of motion of a particle is
dV
m r m r 2
dr
l2
f (r )
mr 3
m r 2 0
0
= constant
giving a straight line through the origin.
3.6 INTEGRABLE POWER LAW IN TIME IN KEPLER’S PROBLEM:
The inverse square law of force is most important of all the central force laws. It
results the deduction of Kepler‟s laws.
1. All the planets move in an elliptical orbit with sun at one of its foci.
2. The radius vector connecting the sun and the planet sweeps at equal areas is equal
intervals of time. ie, areal velocity is constant.
3. The square of the period of revolution of any planet about the sun is proportional
to the cube of the semi major axis.
(i) Deduction of first law:
The central force varies inversely as the square of the distance. That is
k
f (r ) k is a constant
r2
The corresponding potential energy will be
k
V (r )
r
d 2u m 1
Equation of the orbit is u 2 2 f
d 2
l u u
1 k
putting f (r ) f 2 ku2 we get
u r
u cos
mk
(or) u
l2
u 2 u cos
1 mk
r l
Where u and are constants.
If we orient our co-ordinate system so that =0, then
1 mk
2 u cos
r l
mk u l 2
2 1 cos
l mk
l2
(or) r mk .…(2)
u l 2
1 mk cos
Thus is a turning point.
We show that equation (2) represents a conic section, coincident with Kepler‟s
first law of planetary motion. That is the orbits are conic sections with the centre at one of
the foci.
We define the conic section as a curve for which the distance from a fixed point to
that from a fixed line is a constant. That is
r
constant =
d
= eccentricity
p l2
From equation (3) rmin
1 mk (1 )
Now equation (5) becomes
l 2 mk 1 mk 1
2
E k
2m l 2
l2
mk 2 mk 2
= 1 2
1
2l 2 l2
mk 2
= 1 1 2
2l 2
mk 2
= 1 ( 1)
2l 2
=
mk 2 2
2l 2
1
2l 2E 2El 2
giving 2 1 2
1
mk 2 mk 2
1
2El 2
2
such that = 1 2
.…(6)
mk
now is a known one. After putting the value of and p in equation (3) we get
l2
r mk ….(7)
1/ 2
2E l 2
1 1 cos
mk 2
l2
But p =
mk
l2
Hence a =
mk 1 2 .…(8)
l2
(or) 1 - 2 =
mka
Now E
mk 2 2
2l 2
1 becomes
mk 2 l2
= From equation (8)
2l 2 mka
k
E =
2a
which shows that all ellipses with the same major axis have the same energy.
(iv) Deduction of Kepler’s second law:
We know
p
d L d
dt dt
mr 2 0
d 2
dt
r 0
dA r rd r 2 d
1 1
2 2
dA 1 2 d
r
dt 2 dt
1 d
r2
2 dt
1 2
= r = constant
2
which is Kepler‟s second law.
(v) Deduction of third law : (period of elliptic motion )
is the ratio of the total area of the ellipse to the rate at which the area is
swept out. Suppose dA is the area swept out by radius vector is time dt, then
the rate will be dA / dt.
area
Hence
dA / dt
Area of the ellipse = a b
dA 1 2 l
And r
dt 2 2m
ab 2mab
l / 2m l
b a 1 2
1
2
putting
2ma 2
1
1 2
2
Then
l
4 2 m 2 a 4
2
l2
1 2
4 2 m 2 a 4 l 2
=
l2 mka
Assumptions:
1. The heavy nucleus and the positively charged particle to be point nucleus so that their
dimensions are not taken into account.
2. The nucleus of the atom is so heavy that it is at rest during collision.
3. The mass of the positively charged particle may be taken as constant because the
velocity is very small compared with the velocity of light.
The perpendicular distance of PO from N =MN = s. This is the shortest distance from the
nucleus to the initial direction called impact parameter.
1 cos( )
1 1
r p
1/ 2
l2 2 El 2
p and 1 2
mk mk
1 mzze2
2
p l
mzze2
1
1 cos( )
r l2
mzze2
1
1 cos
r l2
1/ 2
2 El 2
1 2
mk
1/ 2
2 El 2
1 2 2
m( zze )
1/ 2
2 El 2
1 2 2 4
mz z e
1
E and mv o 2mE
2
mv o
2
l
l mv o s and mv o
s
mvo s mr 2 l
l
mvo
s
l s 2mE
2 E s 2mE 2
1
1/ 2
mz 2 z 2 e 4
1/ 2
2 Es 2
1
2
zz e
mzze2
1
1 cos
r l2
represents a hyperbola.
The angle between initial and final directions of the positively charged particle is
called angle of scattering. That is the angle between the asymptotes is called angle of
scattering. Here is the angle of scattering.
1 cos 0
1
cos
From figure 2
2 2
1
cos
2 2
1
sin
2
cos ec
2
2
2 Es
1 cot
2
1
2
2
2 zze
2 Es
2
cot 2
2 zze 2
2Es
cot
2 zze 2
2Es
2 cot 1
zze
2
This equation gives in terms of impact parameter s, energy E, the charge on the
nucleus ze and the charge on the particle z e.
The incident intensity is defined as the number of particles crossing unit area normal to
the incident beam in unit time.
The differential of solid angle dω in the plane whose azimuth lies between and
+ d is d 2 sin d
The number of particles scattered into solid angle dω per unit time 2 I () sin d
The cross section for the particles having collision parameter between s and s+ds is the
area of ring of radius s and width ds so that
Manonmaniam Sundaranar University, Directorate of Distance & Continuing Education, Tirunelveli.
93
(s, ds) 2sds
sds
sin d
2 Es
cot
2 zze 2
zze 2
gives s cot
2E 2
zze2 d
Then ds cos ec 2
2E 2 2
zze 2 zze 2 d
cot cos ec 2
()
2E 2 2E 2 2
2 sin cos d
2 2
2
1 zze 2
cos ec 4
4 2E 2
The scattering cross section must be proportional to i) cos ec 4 , ii) the square of (ze),
2
iii) the square of (ze) and iv) inversely proportional to square of kinetic energy (E).
transformation as well as a method itself and applied when Hamiltonian involves time.
H H
Pk and Q k …(1)
Qk Pk
If we require that the transformed Hamiltonian H′ =0, then equations of motion (1)
assume the form
Pk 0 and Q k 0
Thus the coordinates and momenta are constants in time and they are cyclic.
F
H H
t
F
H ( qk , pk , t ) 0 ...(3)
t
We take the generating function F as a function of the old coordinate qk , the new
F2
pk …(4)
qk
F2 F
Therefore H (qk , , t) 2 0 …(5)
qk t
F2 F2 (qk , Pk , t )
dF2 n F n F F
2 qk 2 Pk 2
dt k 1 qk k 1 P
k t
F F2
Here Pk 0, 2 H and pk
t qk
dF2 n
Therefore pk qk H L
dt k 1
S S
Put F2 S in equation (5), we get H (qk , ,t) 0 …(7)
qk t
ie Pk k …(9)
s
Qk k
Pk
Similarly one can calculate the constants k by using initial condition at t = t0, the
known initial values of qk . Thus k and k constants are known and equation (10) will
Thus we see that the Hamilton‟s principal function S is the generaor of a canonical
transformation to constant coordinates (𝛽𝑘 ) and momenta (𝑞𝑘 ). Also in solving the
Hamilton- Jacobi equation, we obtain simultaneously a solution to the mechanical
problem.
JACOBI METHOD:-
Consider a one - dimensional Harmonic oscillator. For such a system forces are
conservative. The force acting on the oscillator at a displacement q is
𝑓 = −𝑘𝑞
k =Force constant.
q
1
Potential energy V kqdq kq2
o
2
1 𝑝2
Kinetic energy 𝑇= 𝑚𝑣 2 =
2 2𝑚
Manonmaniam Sundaranar University, Directorate of Distance & Continuing Education, Tirunelveli.
97
Hamiltonian 𝐻 =𝑇+𝑉
𝑝2 1
= + 𝑘𝑞 2
2𝑚 2
𝜕𝑆
But 𝑃 = and therefore
𝜕𝑞
1 𝜕𝑆 2 1
𝐻= + 𝑘𝑞 2
2𝑚 𝜕𝑞 2
𝜕𝑆
𝐻+ =0
𝜕𝑡
1 𝜕𝑆 2 1 𝜕𝑆
+ 𝑘𝑞 2 + =0 …(1)
2𝑚 𝜕𝑞 2 𝜕𝑡
𝑆 𝑞, 𝛼, 𝑡 = 𝑊 𝑞, 𝛼 − 𝛼𝑡 …(2)
α is a constant.
𝜕𝑆 𝜕𝑊 𝜕𝑆
= 𝑎𝑛𝑑 = −𝛼
𝜕𝑞 𝜕𝑞 𝜕𝑡
1 𝜕𝑊 2 1
+ 𝑘𝑞 2 = 𝛼 …(3)
2𝑚 𝜕𝑞 2
1 𝜕𝑊 2 1
= (𝛼 − 𝑘𝑞 2 )
2𝑚 𝜕𝑞 2
𝜕𝑊 2 1
= 2𝑚 (𝛼 − 𝑘𝑞 2 )
𝜕𝑞 2
𝜕𝑊 1
= 2𝑚 (𝛼 − 𝑘𝑞 2 )
𝜕𝑞 2
On integrating we get
1
W = 2m ( kq2 ) dq c
2
c = constant of integration
Manonmaniam Sundaranar University, Directorate of Distance & Continuing Education, Tirunelveli.
98
1
Then S = 2m ( kq2 ) dq t c
2
c is an additive constant and will not affect the transformation. Because to obtain the
𝜕𝑆
new positon coordinate 𝛽 = only partial derivative of S w. r. to α is required. Hence
𝜕𝛼
c is dropped. Thus
1
S= 2𝑚 (𝛼 − 𝑘𝑞 2 ) 𝑑𝑞 − 𝛼𝑡 …(4)
2
S 2m dq
2
t
1 2
kq
2
m dq
2
t
1 2
kq
2
m dq
2
t
k 2
kq2
2 k
m dq
k 2
t
kq2
m dq
k
t
2
kq2
k
𝑚 𝑘
= sin−1 𝑞
𝑘 2𝛼
𝑘 𝑘
sin−1 = (𝑡 + 𝛽)
2𝛼 𝑚
𝑘
sin−1 𝑞 = 𝜔 (𝑡 + 𝛽)
2𝛼
2𝛼
𝑞= sin 𝜔 (𝑡+ 𝛽) …(5)
𝑘
= 𝑞0 sin 𝜔 (𝑡 + 𝛽)
𝜕𝑆 𝜕𝑊
𝑁𝑜𝑤 𝑝 = =
𝜕𝑞 𝜕𝑞
1
= 2𝑚(𝛼 − 𝑘𝑞 2 )
2
= 2𝑚 𝛼 − 𝑚2 𝜔 2 𝑞 2 …(6)
2𝛼
𝑝= 2 𝑚 𝛼 − 𝑚2 𝜔 2 sin2 𝜔 (𝑡 + 𝛽)
𝑘
= 2 𝑚 𝛼 − 2 𝑚 𝛼 sin2 𝜔 (𝑡 + 𝛽)
= 2 𝑚 𝛼 cos 𝜔 𝑡 + 𝛽 …(7)
= 𝑝0 cos 𝜔 𝑡 + 𝛽
From (6) 2𝑚 𝛼 − 𝑚2 𝜔 2 𝑞 2 = 0
2 𝑚 𝛼 = 𝑚2 𝜔2 𝑞 2
1
𝛼 = 𝑚𝜔2 𝑞 2
2
1
= 𝑘 𝑞02 …(8)
2
𝑞 = 𝑞0 sin 𝜔 (𝑡 + 𝛽) …(9)
𝜋 𝜋
𝜔𝛽= (𝑜𝑟) 𝛽 =
2 2𝜔
Thus the new constant canonical coordinate measures the initial phase angle and in the
𝜋
present initial conditions the initial phase 𝜔𝛽 = .
2
Hence 𝐻 = 𝐸 Thus the new canonical momentum (𝑃) is the total energy of the oscillator
𝑆= 𝑝𝑑𝑞 − 𝛼𝑡
𝑃= 2 𝑚 𝛼 cos 𝜔 (𝑡 + 𝛽)
𝑞 = 𝑞0 sin 𝜔 (𝑡 + 𝛽)
2∝
= sin 𝜔 𝑡 + 𝛽
𝑘
2𝛼
𝑑𝑞 = 𝜔 cos 𝜔 𝑡 + 𝛽 𝑑𝑡
𝑘
2𝛼
𝑆= 2𝑚𝛼 cos 𝜔 𝑡 + 𝛽 𝜔 cos 𝜔 𝑡 + 𝛽 𝑑𝑡 − 𝛼𝑡
𝑘
= +2𝛼 𝑐𝑜𝑠 2 𝜔 𝑡 + 𝛽 𝑑𝑡 − 𝛼𝑡
1
= 2𝛼 [𝑐𝑜𝑠 2 𝜔 𝑡 + 𝛽 − 2] 𝑑𝑡
𝑝2 1
𝐿= − 𝑘𝑞 2
2𝑚 2
1
= 𝛼𝑐𝑜𝑠 2 𝜔 𝑡 + 𝛽 − 𝑘 𝑞0 2 𝑠𝑖𝑛2 𝜔 𝑡 + 𝛽
2
1 2𝛼
= 𝛼𝑐𝑜𝑠 2 𝜔 𝑡 + 𝛽 − 𝑘 𝑠𝑖𝑛2 𝜔 𝑡 + 𝛽
2 𝑘
= 𝛼 𝑐𝑜𝑠 2 𝑤 𝑡 + 𝛽 − 𝛼 𝑠𝑖𝑛2 𝜔 𝑡 + 𝛽
= 𝛼[𝑐𝑜𝑠 2 𝜔 𝑡 + 𝛽 − 𝑠𝑖𝑛2 𝜔 𝑡 + 𝛽 ]
= 𝛼[𝑐𝑜𝑠 2 𝜔 𝑡 + 𝛽 − (1 − 𝑠𝑖𝑛2 𝜔 𝑡 + 𝛽 )]
= 𝛼[2𝑐𝑜𝑠 2 𝜔 𝑡 + 𝛽 − 1]
1
= 2𝛼[𝑐𝑜𝑠 2 𝜔 𝑡 + 𝛽 − ]
2
Thus for Harmonic oscillator we prove that the Hamilton‟s principal function is the time
integral of Lagrangian.
𝐹2 = 𝐹2 (𝑞𝑘, 𝑃𝑘, 𝑡)
F F2
Here Pk 0 , 2 H and pk
t qk
dF2 n
Therefore pk qk H L
dt k 1
(or) F2 Ldt S
Conservative systems:-
In a particle moving under central force in which H does not depend on time t explicitly
and hence H = constant = 𝛼1 = E. In such cases, it is Hamilton‟s characteristic function
W. In Hamilton‟s principal function S, an explicit dependence on time is involved.
Hamilton – Jacobi equation for Hamilton‟s principal function S (𝑞𝑘, 𝛼𝑘, 𝑡) becomes
𝜕𝑆 𝜕𝑆
𝐻 𝑞𝑘, 𝜕𝑞 + =0 …(1)
𝑘 𝜕𝑡
𝜕𝑆 𝜕𝑊 𝜕𝑆
There fore = 𝑎𝑛𝑑 = −𝛼1
𝜕𝑞 𝑘 𝜕𝑞 𝑘 𝜕𝑡
𝜕𝑊 𝜕𝑊 𝜕𝑊
(or) H [𝑞1, 𝑞2,… 𝑞𝑛 , 𝜕𝑞 , 𝜕𝑞 , … 𝜕𝑞 ] = 𝛼1 …(4)
1 2 𝑛
𝜕𝑊 𝜕𝑊
𝑝𝑘 = 𝑎𝑛𝑑 𝑄𝑘 = 𝜕𝑃 …(6)
𝜕𝑞 𝑘 𝑘
Now if the new momenta 𝑃𝑘 are all constant of motion 𝛼𝑘 , where 𝛼1 in particular is the
𝜕𝑊
constant of motion H, then 𝑄𝑘 = 𝜕𝛼
𝑘
𝐻 𝑞𝑘 , 𝑝𝑘 = 𝛼1
𝜕𝑊
Using 𝑝𝑘 = , we obtain
𝜕𝑞 𝑘
𝜕𝑊
𝐻 𝑞𝑘, = 𝛼1
𝜕𝑞 𝑘
𝜕𝑊
Also 𝐻 ′ = 𝐻 + 𝜕𝑡
𝐻 ′ = 𝐻 = 𝛼1 (= 𝐸) …(7)
It generates a canonical transformation where all new coordinates 𝑄𝑘 are cyclic because
𝐻′ = 𝛼1, depending on one of the new momenta 𝑃1 = 𝛼1, and does not contain any 𝑄𝑘 .
𝜕𝐻 ′
Now 𝑃𝑘 = − 𝜕𝑄 = 0 𝑜𝑟 𝑃𝑘 = 𝛼𝑘
𝑘
𝜕𝐻 ′
and 𝑄𝑘 = 𝜕𝛼 = 1 𝑓𝑜𝑟 𝑘 = 1 …(8)
𝑘
𝜕𝑊
𝑄1 = 𝑡 + 𝛽1 = 𝜕𝛼 𝑓𝑜𝑟 𝑘 = 1
1
𝜕𝑊
and 𝑄𝑘 = 𝛽𝑘 = 𝜕𝛼 𝑓𝑜𝑟 𝑘 ≠ 1 …(9)
𝑘
Thus out of all the new coordinates 𝑄𝑘, 𝑄1 is the only coordinate which is not a constant
of motion. Here we observe the conjugate relationship between the time as the new
coordinate and Hamiltonian as the conjugate momentum.
The function W has a physical significance similar to the Hamilton‟s principal function S.
since W (𝑞𝑘, 𝑃𝑘 ) does not involve time t explicitly, its total time derivative is
dW n W n W
q k Pk
dt k 1 qk k 1 P
k
dW n
pk q k
dt k 1
W pk q k dt pk dqk …(10)
k k
and S Ldt [ pk q k H ] dt
k
= W Hdt
Hdt t
1
Now S W 1t
Manonmaniam Sundaranar University, Directorate of Distance & Continuing Education, Tirunelveli.
105
S (qk , k , t ) W (qk , k ) 1t
1 E = total energy.
When the Hamiltonian does not involve time explicitly, one can solve a
mechanical problem by using either Hamilton‟s principal function (or) Hamilton‟s
characteristic function.
In the motion of a system with many degrees of freedom the Hamilton- Jacobi
equation is completely separable in coordinate variables. We consider conservative
system in which the Hamiltonian does not involve time explicitly. The Hamilton-Jacobi
equation is given by
W W W
H q1 , q2 ,........qn ; , ,......... 1 ….(1)
q1 q2 qn
Wk
H k qk ; ; 1 , 2 ,.... n 1 …(3)
qk
Wk
Each of equations (3) involves only one of the coordinates qk and the corresponding
qk
defined if the orbit equations for all the ( pk , qk ) pairs describe either closed orbits or
periodic functions of q k.
Wk
Jk qk ; 1, 2, ........n dqk …(8)
qk
Since q k is a variable of integration and it will be out, when integration is over. Thus J k
H H ( J1 , J 2,.......J n ) …(12)
H
k J1 , J 2 , ....... J n
J k
k k t k …(15)
k s are constants and functions of action variables only and are frequencies of the
periodic motion. Here the frequencies of the periodic motion can be obtained without
complete solution of the problem.
Suppose that the change in angle variable ω k with the completion of one cycle by q l is
l ω k, then
k
l k dql
ql
W
Using k we get
J k
2Wk
l k dql
ql J k
W
J k q dq l
l
J k p dq
l l
J l
kl
J k
Thus k k 1 k=l
=0 k≠l …(16)
But k k 1 hence
1
k k 1 (or ) k k …(18)
k
k
For a particle of mass m, moving in an inverse square force field [V (r ) ], the
r
Hamiltonian of the system is given by
1 2 p k
2
H r
p
2m r2 r
Determine the frequency by the method of action-angle variables and discuss degeneracy.
Show that the period of the orbit is given by
m
k
2E 3
Solution :
J k pk dq k
J p d
and J r pr dr
W
p 2
and therefore
2
W
J d 2 d 2 d 2 2
0
W
J r pr d r dr
r
2mk 2 2
2m1 2 dr
r r
The motion is bounded and is elliptical path for negative value of the total energy
E. Further the limits are given by rmin and rmax values of r.
2 dr
2mk 22
2 .
r2 2m E 2
r r
2mk 22
ie 2m1 2 0
r r
k 2
r2 r 2 0
1 2m1
k 1 2 2
r k2 1 2
21 21 m
k 2 2
= 1 1 1 2 2
21 mk
k 1 21 22
ie rmin 1 and
21 mk 2
k 1 21 22
rmax 1
21 mk 2
Now
rmax
2mk 2
2
Jr = 2
rmin
2m1
r
2 dr
r
2mk 2
2
rmax 2m1 2
=2 r r dr
2mk 2
2
2m1 2
rmin
r r
rmax
2m1r 2mk dr
=2
rmin 2m1r 2mkr
2 2
dr 2 22
r 2m1r 2mkr 22
2
2
2mk
Jr = 2 2
2m1
J + Jr =
2mk
2m1
2mk 2
2m1 J J r 2
22 mk 2
1
J J r 2
But H = H = 1= E
22 mk 2
Therefore H E 1
J J r 2
H 22 mk 2 42 mk 2
2
J J J J r J J r
3
H 22 mk 2 42 mk 2
and r
J r J r J J r 2 J J r 3
42 mk 2
Thus r
J 0 J r 2
The two frequencies are equal and the motion of the system is said to be
degenerated.
Manonmaniam Sundaranar University, Directorate of Distance & Continuing Education, Tirunelveli.
111
The period of the orbit:
1 J J r
3
r 42 mk 2
22 mk 2
E
J J r 2
22 mk 2
J J r 2
E
1
22 mk 2 2
( J J r )
E
3
22 mk 2
2
1
2 2 E
4 mk
1
22 mk 2 22 mk 2
2
1
=
42 mk 2 E E
1
1 22 mk 2
2
2E E
=
m
2
= k 3
2E
k
This formula agrees with Kepler‟s third law that the semi-major axis a .
2E
Cell A Cell B
4 0
3 1
2 2
1 3
0 4
ab cd abc d
bc da bcd a
cd ab cda b
da bc dab c
If we interchange any two phase points from different cells we have different
microstates but the same macrostate.
If we interchange any two phase points in the same cell , we have same microstate
and same macrostate.
The number of microstates corresponding to a given macrostate in called the
thermodynamic probability of the macrostate.
The probability that the ensemble possessing energy E is proportional to Ω(E).
𝑃 𝐸 = 𝑐 𝛺(𝐸)
c→proportionality constant
Ω(E)→thermodynamic probability
The probability of finding one mole of the gas in the container has volume V1+V2
is
𝑉1 +𝑉2 𝑁 𝑁
𝛺2 = = 1 .…(11)
𝑉1 +𝑉2
𝑁
1
= 𝑘 𝑙𝑜𝑔 𝑉1
𝑉 1 +𝑉 2
𝑉1 +𝑉2 𝑁
∆𝑆 = 𝑘 𝑙𝑜𝑔 𝑉1
𝑉1 +𝑉2 𝑁𝑘
= 𝑙𝑜𝑔 ….(12)
𝑉1
The change in entropy when the gas changes from one state with volume V1 and
temperature T1 to another state with volume V2 and temperature T2 is given by,
𝑇2 𝑉1 + 𝑉2
∆𝑆 = 𝐶𝑣 𝑙𝑜𝑔 + 𝑅 𝑙𝑜𝑔 [ ]
𝑇1 𝑉1
𝑇
For isothermal change 𝑇2 = 𝑇1 and hence 𝐶𝑣 𝑙𝑜𝑔 𝑇2 = 0
1
𝑉1 + 𝑉2
∆𝑆 = 𝑅 𝑙𝑜𝑔 [ ]
𝑉1
𝑉1 +𝑉2 𝑅
= 𝑙𝑜𝑔 .…(13)
𝑉1
Canonical ensemble:
We have 𝑑𝑥 𝑑𝑦 𝑑𝑧 = 𝑉
So the volume in phase space = 𝑉 𝑑𝑝 𝑥 𝑑𝑝 𝑦 𝑑𝑝 𝑧
𝑝2
𝜀= → 𝑝2 = 2𝑚𝜀
2𝑚
2𝑝𝑑𝑝 = 2𝑚𝑑𝜀
𝑚
𝑑𝑝 = 𝑑𝜀
𝜑
𝑚 𝑚
= 𝑑𝜀 = 𝑑𝜀
2𝑚𝜀 2𝜀
The sum of the probabilities of the all possible state must be equal to one.
𝑚
ie, 𝑃1 + 𝑃2 + ⋯ … . . +𝑃𝑖 + ⋯ … … . +𝑃𝑚 = 𝑖=1 𝑃𝑖 =1 ….(2)
This called normalization condition.
Now equation (1) becomes
𝑚
𝑢= 𝑖=1 𝑃𝑖 𝑢𝑖
Liouville‟s theorem gives information about the rate of change of phase density in
the phase space. The theorem may be stated in two parts.
The rate of change of density of phase points in the neighborhood of a moving
phase point in the Γ space is zero. This part represents the principle of
conservation of density in the phase space.
𝑑𝜌/𝑑𝑡 = 0 ….(1)
𝑑 𝑑 𝑓
𝛿𝛤 = 𝑑𝑡 𝑖 𝑑𝑞𝑖 𝑑𝑝𝑖 = 0 ….(2)
𝑑𝑡
The change in number of phase points in volume element per unit time,
This change in the number of phase points in the given hyper volume is due to the
difference between the number of phase points entering the hyper volume through
any face and the number of those leaving the opposite face per second.
Consider two faces of hyper volume with co-ordinates 𝑞1 and 𝑞1 + 𝛿𝑞1 . If 𝑞1 is
the component of velocity of phase point at 𝑞1 , 𝑞2 , … … . 𝑞𝑓 , 𝑝1 , 𝑝2 , … … … 𝑝𝑓 ,then
the number of phase points entering the first face AD per second
= 𝜌𝑞1 𝛿𝑞2 … … . . 𝛿𝑞𝑓 𝛿𝑝1 … … . 𝛿𝑝𝑓 ….(5)
As density ρ changes with change in position and momentum co-ordinates and at
the opposite face BC the co-ordinate q1 changes to 𝑞1 + 𝛿𝑞1 and the density ρ
𝜕𝜌
changes to (𝜌 + 𝜕𝑞 𝛿𝑞1 ) at the face BC. The velocity 𝑞 1 changes to (𝑞1 +
1
𝜕𝑞 1
𝛿𝑞1 ). Therefore the number of phase points leaving the opposite face BC at
𝜕𝑞 1
𝜕𝜌 𝜕𝑞
= (𝜌 + 𝜕𝑞 𝛿𝑞1 ) (𝑞1 + 𝜕𝑞 1 𝛿𝑞1 ) 𝛿𝑞2 … … 𝛿𝑞𝑓 𝛿𝑝1 , … … … 𝛿𝑝𝑓
1 1
𝜕𝑞 1 𝜕𝜌
= 𝜌𝑞1 + (𝜌 + 𝑞1 𝜕𝑞 )𝛿𝑞1 𝛿𝑞2 … … 𝛿𝑞𝑓 𝛿𝑝1 , … … … 𝛿𝑝𝑓 ....(6)
𝜕𝑞 1 1
Subtracting (6) from (5) we get the expression for change in the number of phase
points per second corresponding to q1.
𝜕𝑞 1 𝜕𝜌
= −(𝜌 + 𝑞1 𝜕𝑞 )𝛿𝑞1 𝛿𝑞2 … … 𝛿𝑞𝑓 𝛿𝑝1 , … … … 𝛿𝑝𝑓 ….(7)
𝜕𝑞 1 1
Similarly, the expression for the change into the number of phase points per
second corresponding to p1 is
𝜕𝑝 1 𝜕𝜌
= −(𝜌 + 𝑝1 𝜕𝑝 )𝛿𝑞1 𝛿𝑞2 … … 𝛿𝑞𝑓 𝛿𝑝1 , … … … 𝛿𝑝𝑓 ….(8)
𝜕𝑝 1 1
Since the change in number of phase points per second corresponding to all
position and momentum coordinates are like equation (7) and (8), then they are
summed up.
𝑑(𝛿𝑁) 𝑓 𝜕𝑞 𝑖 𝜕𝑝 𝜕𝜌 𝜕𝜌
=− 𝑖=1 𝜌 + 𝜕𝑝 𝑖 + ( 𝑞1 𝜕𝑞 + 𝑝𝑖 ) 𝛿𝑞1 … 𝛿𝑞𝑓 𝛿𝑝1 , … 𝛿𝑝𝑓 .....(9)
𝑑𝑡 𝜕𝑞 𝑖 𝑖 𝑖 𝜕𝑝 𝑖
𝜕𝜌 𝑓 𝜕𝑞 𝑖 𝜕𝑝 𝜕𝜌 𝜕𝜌
=− 𝑖=1 𝜌 + 𝜕𝑝 𝑖 + ( 𝑞1 𝜕𝑞 + 𝑝𝑖 ) ….(10)
𝜕𝑡 𝜕𝑞 𝑖 𝑖 𝑖 𝜕𝑝 𝑖
𝜕2𝐻 𝜕2𝐻
=
𝜕𝑞𝑖 𝜕𝑝𝑖 𝜕𝑝𝑖 𝜕𝑞𝑖
𝜕𝑞 𝑖 𝜕𝑝
We get = − 𝜕𝑝 𝑖
𝜕𝑞 𝑖 𝑖
𝑓 𝜕𝑞 𝑖 𝜕𝑝
𝑖=1 𝜕𝑞 𝑖 + 𝜕𝑝 𝑖 = 0 ….(11)
𝑖
𝜕𝜌 𝑓 𝜕𝜌 𝜕𝜌
=− 𝑖=1 𝑞1 𝜕𝑞 + 𝑝𝑖
𝜕𝑡 𝑞,𝑝 𝑖 𝜕𝑝 𝑖
𝜕𝜌 𝑓 𝜕𝜌 𝜕𝜌
+ 𝑖=1 𝑞1 + 𝑝𝑖 = 0 ….(12)
𝜕𝑡 𝑞,𝑝 𝜕𝑞 𝑖 𝜕𝑝 𝑖
𝑑𝜌
(𝑞 , … . , 𝑞𝑓 , 𝑝1 , … … . , 𝑝𝑓 , 𝑡) = 0
𝑑𝑡 1
𝑑𝜌
ie, 𝑑𝑡 = 0 .…(13)
𝑑 𝑑
(𝛿𝑁) = (𝜌. 𝛿𝛤)
𝑑𝑡 𝑑𝑡
𝑑𝜌 𝑑(𝛿𝛤 )
= 𝛿𝛤 + 𝜌 ....(14)
𝑑𝑡 𝑑𝑡
As each phase point represents a definite system and systems can neither be
created nor destroyed, the number of phase points δN must remain fixed.
𝑑
i.e 𝛿𝑁 = 0
𝑑𝑡
𝑑𝜌 𝑑(𝛿𝛤)
𝛿𝛤 + 𝜌 =0
𝑑𝑡 𝑑𝑡
𝑑𝜌
from equation (14) =0
𝑑𝑡
𝑑(𝛿𝛤 )
𝜌 =0
𝑑𝑡
𝑑(𝛿𝛤)
=0
𝑑𝑡
𝛿𝛤 = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡 ….(15)
STATIONARY ENSEMBLE:
4.10 MICRO CANONICAL ENSEMBLE (ISOLATED SYSTEM):
An ensemble in which each system has the same fixed energy as well as the same
number of particles is called micro canonical ensemble.
In this ensemble, density ρ, for a closed isolated thermo dynamical system is a
function of energy and we take
ρ(E) = constant between the energy shells E and E+δE of phase space.
=𝐴 𝑖 𝑔𝑖 𝑒 −𝜀 𝑖 /𝑘𝑇
𝑁
= 𝑖 𝑔𝑖 𝑒 −𝜀 𝑖 /𝑘𝑇
𝐴
𝑍= 𝑖 𝑔𝑖 𝑒 −𝜀 𝑖 /𝑘𝑇
Z is known as partition function and Z indicates how the gas molecules of an
assembly are distributed (or) partitioned among the various energy levels.
If the energy of the 𝑖 𝑡 level is 𝜀𝑖 then the weight of an individual level is unity.
ie, 𝑔𝑖 = 1
𝑍= 𝑖 𝑒 −𝜀 𝑖 /𝑘𝑇
Here the energy term may contain the rotational, vibrational and electronic
components in addition to translational component.
„Z‟ can be used for calculating the various thermodynamic properties of
ensembles.
In classical treatment the energy distribution is continuous.
The number of energy levels of the momentum interval p and 𝑝 + 𝑑𝑝 is given by,
𝑉 4𝜋𝑝2 𝑑𝑝
𝑔(𝑝)𝑑𝑝 =
3
𝑝2 = 2𝑚𝜀
2𝑝𝑑𝑝 = 2𝑚𝑑𝜀.
𝑚
𝑑𝑝 = 𝑑𝜀.
𝑝
𝑚 𝑚
= 𝑑𝜀 = 𝑑𝜀
2𝑚𝜀 2𝜀
Now the number of energy levels in the energy range 𝜀 and 𝜀 + 𝑑𝜀 is obtained as,
𝑉 𝑚
𝑔 𝜀 𝑑𝜀 = 4𝜋 2𝑚𝜀 𝑑𝜀
3 2𝜀
2𝜋𝑉 3/2
= 2𝑚 𝜀 1/2 𝑑𝜀
3
𝑍= 𝑖 𝑔𝑖 𝑒 −𝜀 𝑖 /𝑘𝑇 .
∞
= 0
𝑔(𝜀) 𝑒 −𝜀/𝑘𝑇 𝑑𝜀
∞ 2𝜋𝑉 3/2
= 0 3
2𝑚 𝜀1/2 𝑒 −𝜀/𝑘𝑇 𝑑𝜀
2𝜋𝑉 3/2 1 𝜋
= 2𝑚
3 2 𝛽3
2𝜋𝑉 3/2 1
= 2𝑚 𝜋(𝑘𝑇)3
3 2
𝑉
𝑍= (2𝜋𝑚𝐾𝑇)3/2
3
This is the translational partition function for a gas molecule.
4.11 CLASSICAL IDEAL GAS USING MICRO CANONICAL ENSEMBLE:
Consider a micro canonical ensemble of a perfect gas.
Let there be n point particles with mass m confined in a volume V with total
energy u within the energy range 𝛿𝑢.
The corresponding volume
∆𝛤 = 𝑑 𝑞1 … … … . 𝑑𝑞3𝑛 𝑑 𝑝1 … … … . 𝑑𝑝3𝑛
𝑑 𝑞1 … … … 𝑑𝑞3𝑛 = 𝑉 𝑛 .
Hence ∆𝛤 = 𝑉 𝑛 𝑑 𝑞1 … … … … 𝑑𝑞3𝑛
The momentum space integral is to be evaluated subject to the constraint of the
ensemble
𝑢 − 𝛿𝑢 ≤ 𝑢𝑟 ≤ 𝑢.
𝑛 2
𝑢𝑟 = 𝑖=1 𝑝𝑖 /2𝑚 .
1 𝑛 2
𝑢 − 𝛿𝑢 ≤ 2𝑚 𝑖=1 𝑝𝑖 ≤ 𝑢.
𝜋 𝑓/2
𝑉𝑓 (𝑅) = 𝑅 𝑓 = 𝑐𝑓 𝑅 3
𝑓/2 !
𝜋 𝑓/2
where 𝐶𝑓 = 𝑓/2 !
𝜋 3𝑛 /2 𝛿𝑢 3𝑛/2
= (2𝑚𝑢)3𝑛/2 [1 − (1 − ) ]
3𝑛/2 ! 𝑢
𝜋 3𝑛 /2 3𝑛 𝛿𝑢
= (2𝑚𝑢)3𝑛/2 [1 − 𝑒𝑥𝑝(− . )]
3𝑛/2 ! 2 𝑢
3𝑛 𝛿𝑢
For a macroscopic system 3n=1023 ; >>u.
2 𝑢
∆𝛤 = 𝑉 𝑛 𝑑𝑝1 … … … … … … . . 𝑑𝑝3𝑛
𝜋 3𝑛 /2
= 𝑉𝑛 (2𝑚𝑢)3𝑛/2
3𝑛/2 !
𝑛
𝜋 3𝑛 /2
= 𝑙𝑜𝑔 [𝑉 (2𝑚𝑢)3𝑛/2 ]
3𝑛/2 !
= 𝑛 𝑙𝑜𝑔[𝑉 𝜋 3/2 2𝑚𝑢 3/2
] − 𝑙𝑜𝑔 (3𝑛/2)!
= 𝑛 𝑙𝑜𝑔[𝑉 𝜋 3/2 2𝑚𝑢 3/2
] − (3𝑛/2) 𝑙𝑜𝑔 (3𝑛/2) + 3𝑛/2
= 𝑛 𝑙𝑜𝑔 [𝑉 𝜋 3/2 2𝑚𝑢 3/2
] − 𝑛 𝑙𝑜𝑔 3𝑛/2 3/2
+ 3𝑛/2
𝑉 𝜋 3/2 2𝑚𝑢 3/2
= 𝑛 𝑙𝑜𝑔 + 3𝑛/2
3𝑛 /2 3/2
We know that the entropy should not depend on the unit of hyper volume ∆Γ. To
make it dimensionless we divide it by 3𝑛 .
𝜍 = 𝑙𝑜𝑔[∆𝛤/3𝑛 ]
The above equation does not satisfy the additive property and hence to satisfy the
additive property we must divide by n!
This expression satisfies the additive property because instead of V and u we have
V/n and u/n.
We shall now establish the connection of statistical quantities with corresponding
thermodynamic quantities.
(a) Internal energy(U):
1 𝜕𝜍
=
𝜏 𝜕𝑢 𝑇,𝑛
𝜕 4𝜋𝑚 3/2 3 3
= [𝑛 𝑙𝑜𝑔 𝑣– 𝑛 𝑙𝑜𝑔 𝑛 + 𝑛 𝑙𝑜𝑔 + 2 𝑛 𝑙𝑜𝑔 𝑢 − 2 𝑛 𝑙𝑜𝑔 𝑛 – 𝑛𝑙𝑜𝑔 3 ] +
𝜕𝑢 3
𝜕 5
( 𝑛)
𝜕𝑢 2
𝜕 3
= ( 𝑛 𝑙𝑜𝑔 𝑢)
𝜕𝑢 2
1 3 1
= 𝑛
𝜏 2 𝑢
3 3
𝑢 = 2 𝑛𝜏 (or) 𝑢 = 2 𝑛 𝑘𝑇
𝜏 = 𝑘 × 𝑡𝑒𝑟𝑚𝑜𝑑𝑦𝑛𝑎𝑚𝑖𝑐 𝑡𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑒.
𝜏 = 𝑘𝑇
𝑝 𝜕𝜍
We have =
𝜏 𝜕𝑉 𝑛.𝑢
3
𝜕 4𝜋𝑚 2 3 3
= 𝑛 𝑙𝑜𝑔 𝑣 − 𝑛 𝑙𝑜𝑔 𝑛 + 𝑛 𝑙𝑜𝑔 + 2 𝑛 𝑙𝑜𝑔 𝑢 – 2 𝑛 𝑙𝑜𝑔 𝑛 −
𝜕𝑉 3
𝜕
= [ 𝑛 𝑙𝑜𝑔 𝑉]
𝜕𝑉
= 𝑛/𝑉
𝑃𝑉 = 𝑛 𝜏 (or) 𝑃𝑉 = 𝑛𝑘𝑇
Which is well known ideal gas equation for a perfect mono atomic gas.
The relation between thermodynamic entropy and statistical entropy is given by,
𝑆 = 𝑘𝜍
𝑉 4𝜋𝑚 3/2 3 5
= 𝑛𝑘 𝑙𝑜𝑔 (2 𝑘𝑇)3/2 + 2 𝑛𝑘
𝑛 3 3
𝑉 3/2
5
= 𝑛𝑘 𝑙𝑜𝑔 (2𝜋𝑚𝑘𝑇) + 𝑛𝑘
𝑛3 2
3/2
𝑉 2𝜋𝑚𝑘𝑇 5
= 𝑛𝑘 𝑙𝑜𝑔 + 𝑛𝑘
𝑛 2 2
𝑉 2𝜋𝑚𝑘𝑇 3/2
= 𝑛𝑘 𝑙𝑜𝑔 𝑒 5/2
𝑛 2
This is the famous Sackur - Tetrode equation for the entropy of a perfect gas.
This formula is valid for the mono atomic gas of atoms with zero total angular
momentum.
The thermal de-broglie wavelength associated with a molecule may be defined as,
𝜆 = /(2𝜋𝑚𝑘𝑇)1/2
𝜆3 = 3 /(2𝜋𝑚𝑘𝑇)3/2
1
= (2𝜋𝑚𝑘𝑇/2 )3/2
𝜆3
𝑉 1 5
Now 𝜍 = 𝑛 𝑙𝑜𝑔 +2𝑛
𝑛 𝜆3
𝑉 1 5
𝑆 = 𝑛𝑘 𝑙𝑜𝑔 + 2 𝑛𝑘
𝑛 𝜆3
Thus the entropy of a perfect gas is determined essentially by the ratio of the
volume per particle to the volume λ3 associated with de-Broglie wavelength.
(e) Chemical potential of a perfect gas:
−𝜇 𝜕𝜍
=
𝜏 𝜕𝑛 𝑢,𝑉
𝜕 𝑉 1 5
= 𝜕𝑛 𝑛 𝑙𝑜𝑔 + 𝑛 𝑢,𝑉
𝑛 𝜆3 2
𝜕 𝜕 5
= 𝜕𝑛 [𝑛 𝑙𝑜𝑔 𝑉 − 𝑛 𝑙𝑜𝑔 𝑛 − 𝑛 𝑙𝑜𝑔 𝜆3 ]𝑢,𝑣 + 𝜕𝑛 𝑛
2 𝑢,𝑣
Manonmaniam Sundaranar University, Directorate of Distance & Continuing Education, Tirunelveli.
132
5
= 𝑙𝑜𝑔 𝑉 − 1 − 𝑙𝑜𝑔𝑛 – 𝑙𝑜𝑔 𝜆3 + 2
𝑉 3
= 𝑙𝑜𝑔 +2
𝑛𝜆 3
𝜇 𝑛𝜆 3 3
= 𝑙𝑜𝑔 −2
𝜏 𝑉
𝑛 𝑝
=
𝑉 𝜏
𝜇 𝑝𝜆 3 3
= 𝑙𝑜𝑔 −2
𝜏 𝜏
𝜆3 3𝜏
µ = 𝜏 𝑙𝑜𝑔 𝑝 + 𝜏 𝑙𝑜𝑔( 𝜏 ) − .
2
= 𝜏 𝑙𝑜𝑔 𝑝 + 𝑓(𝜏)
The micro canonical ensemble describes the systems which are perfectly insulated
and have given energy.
In thermodynamics we do not know the exact value of energy as we usually deal
with systems kept in thermal contact with a heat reservoir at a given temperature.
Thus we know only its temperature i.e its average energy.
The energy varies from instant to instant but the time average is known.
On the other hand the canonical ensemble describes those systems which are not
isolated, but are in thermal contact with a heat reservoir.
In this situation the system of interest together with a heat reservoir forms a large
closed system and the system of interest is treated as a subsystem.
If the energy of the large closed system is constant, then it would represent a
microcanonical system where as the subsystem which can exchange energy with a
heat reservoir would represent canonical system.
Thus any part of sub system of an isolated system in thermal equilibrium can be
represented by a canonical ensemble.
𝐸𝑡 = 𝐸𝑟 + 𝐸𝑠
𝜕𝛤 𝐸
∆Γ== 𝛿𝐸
𝜕𝐸 𝐸
To estimate δE:
So that,
𝜔 𝐸 𝑑𝐸 = 𝜌 𝐸 ∆𝛤(𝐸) ….(1)
𝜕𝛤 𝐸
= 𝜌(𝐸) 𝑑𝐸
𝜕𝐸 𝐸
ω(𝐸 ) δE = 1
∆𝛤 = 𝐴−1 𝑒 𝐸 /𝜏
𝜍 = 𝑙𝑜𝑔 ∆𝛤
= 𝑙𝑜𝑔[𝐴−1 . 𝑒 𝑈/𝑘𝑇 ]
𝑈
= − 𝑙𝑜𝑔 𝐴 + 𝑘𝑇
𝑈
𝑙𝑜𝑔 𝐴 = −𝜍
𝑘𝑇
𝑈−𝑘𝜍𝑇
= 𝑘𝑇
𝑈−𝑆𝑇
= ….(4)
𝑘𝑇
𝐹 = 𝑈 − 𝑇𝜍 = 𝑈 − 𝑆𝑇
𝐹
𝑙𝑜𝑔 𝐴 = 𝑘𝑇
𝐴 = 𝑒 𝐹/𝑘𝑇
𝜌(𝐸) = 𝐴𝑒 −𝐸/𝜏
= 𝑒 𝐹/𝑘𝑇 𝑒 −𝐸/𝑘𝑇
= 𝑒 (𝐹−𝐸)/𝑘𝑇
𝜌 𝐸 𝑑𝛤 = 1
𝑒 (𝐹−𝐸)/𝑘𝑇 𝑑𝛤 = 1
𝑒 −𝐹/𝑘𝑇 = 𝑒 −𝐸(𝑝,𝑞)/𝑘𝑇 𝑑𝛤
𝑍= 𝑖 𝑒 −𝐸𝑖/𝑘𝑇 (quantum)
So, 𝑒 −𝐹/𝑘𝑇 = 𝑍
−𝐹/𝑘𝑇 = 𝑙𝑜𝑔 𝑍
−𝐹 = 𝑘𝑇 𝑙𝑜𝑔 𝑍
This equation represents the expression for the Helmholtz free energy in terms of
Z.
Suppose N independent identical spinless particles, we must correct the classical
partition function dimensionally and take into account the indistinguishability of
the particles so that the correct expression is,
1
𝑍 = 𝑁! 3𝑁 𝑒 −𝐸(𝑝,𝑞)/𝑘𝑇 𝑑𝛤 (classical) ….(6)
Entropy of a system:
𝜕
= (𝜏 log 𝑍)
𝜕𝜏
𝜕
= log 𝑍 + 𝜏 (log 𝑍) ….(7)
𝜕𝜏
𝑍= 𝑖 𝑒 −𝐸𝑖 /𝜏
𝜕 𝜕 −𝐸𝑖 /𝜏
(log 𝑍) = 𝜕𝜏 log 𝑖𝑒
𝜕𝜏
𝐸𝑖
𝑒 −𝐸 𝑖 /𝜏 )
𝜏2
=
𝑖 𝑒 −𝐸𝑖 /𝜏
Manonmaniam Sundaranar University, Directorate of Distance & Continuing Education, Tirunelveli.
137
1 𝐸𝑖 𝑒 −𝐸𝑖 /𝜏
=
𝜏2 𝑒 −𝐸𝑖 /𝜏
1 𝑈
= 𝐸= ….(8)
𝜏2 𝜏2
𝑈
Now statistical entropy 𝜍 = log 𝑍 + ….(9)
𝜏
𝑈
The thermo dynamic entropy 𝑆 = 𝑘𝜍 = 𝑘 log 𝑧 + 𝜏
𝑈
𝑆 = 𝑘 log 𝑍 + 𝑇 ….(10)
N
exp
2
e N / 2 1 e
N
.…(3)
1 exp N
For the Helmholtz free energy of the system, we get
A kT ln e N / 2 1 e
N
NkT ln[e1 2 ] NkT ln[1 e ]
1
N kT ln 1 e ….(4)
2
whereby
µ = A/N. .… (5)
P=0 .… (6)
S kU k lnQ N
1
S kN
2
(e 1)
k ln e N / 2 1 e
N
1
S kN
(e 1)
k ln e N / 2 k ln 1 e
N
2
S kN
2
kN
(e 1)
kN
2
k ln 1 e
N
S kN
ln 1 e
(e 1)
Then
S kN
ln 1 e .…(7)
(e 1)
Also
1 1 1
S Nk coth ln 2 sinh
2 2 2
U ln QN
ln e N / 2 1 e
N
ln e N / 2 ln 1 e
N
1 1 1
U N N coth .… (8)
2 (e 1) 2 2
dU
and C p Cv
dT
d 1
N
dT 2 (e 1)
d
dT (e kT 1)
N
kT
e
N kT
kT 2
(e 1) 2
e
2
Nk
kT e 12
Then
e
C p Cv Nk
2
e
2 .… (9)
1
Also
2
1 1
Cp Cv Nk cos ech 2
2 2
U
U kT 2
2
T V
e
kT Nk
2 2
e
2 [From equation (9)]
1
e
Nk 2T 2
1
2
k 2T 2 e 1
2
Manonmaniam Sundaranar University, Directorate of Distance & Continuing Education, Tirunelveli.
140
e
N
2
e
1 2
e 2
Now U N
.… (10)
e 1
Thus U N
𝛺 = 𝑈 − 𝑇𝑆 − µ𝑛
𝐸𝑠 + 𝐸𝑟 = 𝐸𝑡 and 𝑛𝑠 + 𝑛𝑟 = 𝑛𝑡
Now we find the probability 𝑑𝜔𝑠 𝑛𝑠 of finding the sub systems in a state in
which sub system S contains 𝑛 particles and is found in the element 𝑑𝛤 𝑛𝑠 of its
phase space.
𝑑𝛤 𝑛𝑠 indicates that the nature of the phase space of the subsystem changes with
𝑛𝑠 .
For grand canonical ensemble
𝑑𝜔 𝑛 = 𝐴 exp 𝑛µ − 𝐸 /𝜏 𝑑𝛤 𝑛
𝛺
Normalization constant 𝐴 = exp 𝜏
exp 𝛺 + 𝑛µ − 𝐸 /𝜏 𝑑𝛤 𝑛 = 1
𝐸𝑠 + 𝐸𝑟 = 𝐸𝑡 𝑎𝑛𝑑 𝑛𝑠 + 𝑛𝑟 = 𝑛𝑡 .…(1)
𝜌 𝑛 = exp 𝛺 + 𝑛µ − 𝐸 /𝜏 ….(2)
Ω→grand potential
µ→chemical potential
𝒵= 𝑒 µ𝑛 𝑍𝑛
𝑛
= 𝑒 −𝐸𝑛 ,𝑖 /𝜏 (quantum)
i.e the grand partition function 𝒵 is the sum of canonical partition functions Z i
for ensemble with different n′ s with weighing factor 𝑒 −µ𝑛/𝜏 .
From equation (3) the partition function is given by,
𝒵 = exp −𝛺/𝜏
𝛺 = −𝜏 log 𝒵 ….(4)
𝜍 = log ∆𝛤
1
= log
𝜌 𝑛, 𝐸
= − 𝛺 + 𝑛µ − 𝐸 /𝜏
= − 𝛺 + 𝑛µ − 𝑈 /𝜏 ….(5)
𝜏𝜍 = − 𝛺 + 𝑛µ − 𝑈
𝑈 − 𝜏𝜍 = (𝛺 + 𝑛µ)
= 𝛺 + 𝑛µ ….(7)
𝐺 = 𝐹 + 𝑝𝑉 ….(8)
= 𝑈 − 𝜏𝜍 + 𝑝𝑉 ….(9)
In this case 𝐺 = µ𝑛
𝛺 + 𝑛µ + 𝑝𝑉 = 𝑛 µ
𝛺 = −𝑝𝑉
𝐹 + 𝑝𝑉 = µ𝑛 ….(11)
𝑈 − 𝜏𝜍 + 𝑝𝑉 = µ𝑛
𝑈 − 𝜏𝜍 − 𝛺 = µ𝑛
𝑈 − 𝜏𝜍 − µ𝑛 = 𝛺 ….(12)
𝜕𝛺
Then 𝑝=− 𝜕𝑉 𝜏,µ
𝜕𝛺
𝜍=− 𝜕𝜏 𝑉,𝜇
𝜕𝛺
𝑛=− 𝜕𝜇 𝑉,𝜏
From the above three relations we can evaluate thermodynamic quantities for the
grand canonical ensemble.
𝒵= 𝑛 𝑒 𝜇𝑛 /𝜏 𝑍𝑛 ….(1)
𝑛
𝜇𝑛 /𝜏 𝑓
𝑛 𝑒 𝜇 /𝜏 𝑓
𝒵= 𝑛𝑒 = 𝑛 ....(3)
𝑛! 𝑛!
𝑥𝑛
From the series expansion 𝑛 𝑛! = 𝑒𝑥
𝒵 = exp 𝑒 𝜇 /𝜏 𝑓 ….(4)
𝛺 = −𝜏 log 𝒵 = −𝜏𝑒 𝜇 /𝜏 𝑓
2𝜋𝑚𝜏 3/2
𝛺 = 𝜏𝑒 𝜇 /𝜏 𝑉 ....(5)
2
𝜕𝛺
𝑛=− 𝜕𝜇 𝑉,𝜏
𝜕 2𝜋𝑚𝜏 3/2
= − 𝜕𝜇 −𝜏𝑒 𝜇 /𝜏 𝑉
2
𝑉,𝜏
𝜕 2𝜋𝑚𝜏 3/2
= 𝜕𝜇 𝜏𝑒 𝜇 /𝜏 𝑉
2
𝑉,𝜏
2𝜋𝑚𝜏 3/2 1
=𝜏 𝑉 𝑒 𝜇 /𝜏 𝜏
2
2𝜋𝑚𝜏 3/2
𝑛= 𝑉𝑒 𝜇 /𝜏 ….(6)
2
𝛺
=−𝜏
𝛺 = −𝜏𝑛 ….(7)
2𝜋𝑚𝜏 3/2
𝛺 =−𝜏 𝑉𝑒 𝜇 /𝜏 ….(8)
2
2𝜋𝑚𝜏 3/2
−𝜏𝑛 = −𝜏 𝑉𝑒 𝜇 /𝜏
2
𝜇 𝑓
− 𝜏 = log 𝑛
𝑓
𝜇 = −𝜏 log ….(9)
𝑛
𝜕𝛺
𝜍=− 𝜕𝜏 𝑉,𝜇
𝜕 2𝜋𝑚𝜏 3/2
= 𝜕𝜏 𝜏𝑒 𝜇 /𝜏 𝑉
2
2𝜋𝑚 3/2 𝜕
= 𝑉 𝜏 5/2 𝑒 𝜇 /𝜏
2 𝜕𝜏
2𝜋𝑚 3/2 5 𝜇
= 𝑉𝜏 3/2 𝑒 𝜇 /𝜏 −𝜏
2 2
2𝜋𝑚𝜏 3/2 5 𝜇
𝜍= 𝑉𝑒 𝜇 /𝜏 −𝜏 ….(10)
2 2
5 𝜇
𝜍=𝑛 −𝜏 ….(11)
2
𝜇 𝑓 𝑓
From equation (9) we have = − log → 𝜇 = −𝜏 log
𝜏 𝑛 𝑛
𝑛
2𝜋𝑚𝜏 3/2 𝑙𝑜𝑔 5 𝑓
Now equation (10) becomes 𝜍 = 𝑉𝑒 𝑓 + log
2 2 𝑛
2𝜋𝑚𝜏 3/2 𝑛 5 𝑓
= 𝑉 + log
2 𝑓 2 𝑛
5 𝑓
=𝑛 + log
2 𝑛
𝑙𝑜𝑔2𝜋𝑚𝜏23/2𝑉𝑛 ….(12)
Thermodynamic entropy S = 𝑘𝜍
5 2𝜋𝑚𝜏 3/2 𝑉
= 𝑛𝑘 + 𝑙𝑜𝑔
2 2 𝑛
2𝜋𝑚𝜏 3/2 𝑉
= 𝑛𝑘𝑙𝑜𝑔 𝑒 5/2 ….(13)
2 𝑛
This gives the famous Sackur-Tetrode equation for the entropy of a perfect gas
and this is in agreement with the equation for microcanonical ensemble and
canonical ensemble.
(c) Internal energy:
We have 𝛺 = 𝑈 − 𝜏𝜍 − 𝜇𝑛
𝑈 = 𝛺 + 𝜏𝜍 + 𝜇𝑛
5 𝜇
𝜍=𝑛 −𝜏 [from equation (11)]
2
5 𝜇
Now 𝑈 = −𝑛𝜏 +𝜏 𝑛 2
− 𝜏 + 𝜇𝑛
5
= −𝑛𝜏 + 2 𝑛𝜏 − 𝜇𝑛 + 𝜇𝑛
3
= 2 𝑛𝜏.
3
𝑈 = 2 𝑛𝑘𝑇.
which is well known relation for the internal energy of a perfect gas.
4.16 ENERGY AND DENSITY FLUCTUATIONS IN ENSEMBLES:
𝑍= 𝑖 𝑒 −𝛽𝐸𝑖 ….(1)
1
where 𝛽 = 𝑘𝑇
𝜕𝑍
= 𝑖 𝑒 −𝛽𝐸𝑖 −𝐸𝑖
𝜕𝛽
= 𝑖 𝐸𝑖 𝑒 −𝛽𝐸𝑖
𝜕𝑧
𝑖 𝐸𝑖 𝑒 −𝛽 𝐸𝑖 = − 𝜕𝛽 .
−𝛽 𝐸 𝑖
𝑖 𝐸𝑖 𝑒
Mean energy 𝐸= −𝛽 𝐸𝑖 .
𝑖𝑒
−𝛽 𝐸 𝑖
𝑖 𝐸𝑖 𝑒
= .
𝑍
1 𝜕𝑧
𝐸=− ….(2)
𝑧 𝜕𝛽
𝜕2𝑧
= 𝑖 𝐸𝑖 2 𝑒 −𝛽𝐸𝑖
𝜕𝛽 2
𝜕2𝑧
𝐸𝑖 2 𝑒 −𝛽𝐸𝑖 = 𝜕𝛽 2
2 𝑒 −𝛽 𝐸 𝑖
𝑖 𝐸𝑖
𝐸2 = −𝛽 𝐸 𝑖
𝑖𝑒
2 𝑒 −𝛽 𝐸 𝑖
𝑖 𝐸𝑖
= .
𝑍
1 𝜕2𝑧
𝐸2 = 𝑧 ….(3)
𝜕𝛽 2
𝜕𝐸 𝜕 1 𝜕𝑧
Using equation (2) = 𝜕𝛽 − 𝑧 𝜕𝛽
𝜕𝛽
𝜕 1 𝜕𝑧 1 𝜕2𝑧 1 𝜕𝑧 2
= − 𝜕𝛽 =− − 𝑧2
𝑧 𝜕𝛽 𝑧 𝜕𝛽 2 𝜕𝛽
= − 𝐸2 − 𝐸 2
2
= − 𝛿𝐸 ….(4)
𝜕𝐸 𝜕𝛽
= 𝜕𝛽 𝜕𝑇
2 1
= − 𝛿𝐸 − 𝑘𝑇 2
1 2
𝐶𝑉 = 𝛿𝐸 ….(5)
𝑘𝑇 2
2
𝛿𝐸 = 𝑘𝑇 2 𝐶𝑉
2 1/2
𝛿𝐸 = 𝑇 𝑘𝐶𝑉 .…(6)
2
∆𝐸 𝛿𝐸
=
𝐸 𝐸
𝑇 𝑘𝐶𝑉 1/2
= ….(7)
𝐸
∆𝐸 𝑇 1/2
= 𝑘𝑁𝑘
𝐸 𝑁𝑘𝑇
𝑁 1/2
= = 𝑁 −1/2 ….(8)
𝑁
𝜕𝒵 𝑛
= 𝑛,𝑖 𝑒 (𝑛𝜇 −𝐸𝑛 ,𝑖 )/𝜏
𝜕𝜇 𝜏
1
=𝜏 𝑛,𝑖 𝑛 𝑒 (𝑛𝜇 −𝐸𝑛 ,𝑖 )/𝜏
𝜕𝒵
𝑛𝑒 (𝜇𝑛 −𝐸𝑛 ,𝑖 )/𝜏 = 𝜏 𝜕𝜇
(𝑛𝜇 −𝐸 𝑛 ,𝑖 )/𝜏
𝑛 ,𝑖 𝑛𝑒
Mean concentration 𝑛= (𝑛𝜇 −𝐸 𝑛 ,𝑖 )/𝜏
𝑛 ,𝑖 𝑒
(𝑛𝜇 −𝐸 𝑛 ,𝑖 )/𝜏
𝑛 ,𝑖 𝑒
=
𝒵
𝜏 𝜕𝒵
𝑛= … . (10)
𝒵 𝜕𝜇
𝜕2𝒵 1
= 𝜏2 𝑛,𝑖 𝑛2 𝑒 (𝑛𝜇 −𝐸𝑛 ,𝑖 )/𝜏
𝜕𝜇 2
𝜕2𝒵
𝑛,𝑖 𝑛2 𝑒 (𝑛𝜇 −𝐸𝑛 ,𝑖 )/𝜏 = 𝜏 2 𝜕𝜇 2
2 𝑒 (𝑛𝜇 −𝐸 𝑛 ,𝑖 )/𝜏
𝑛 ,𝑖 𝑛
=
𝑍
𝜏2 𝜕2𝒵
𝑛2 = . … (11)
𝑧 𝜕𝜇 2
𝜕𝑛 𝜕 𝜏 𝜕𝒵
= 𝜕𝜇
𝜕𝜇 𝒵 𝜕𝜇
𝜏 𝜕2𝒵 𝜏 𝜕𝒵 2
= 𝒵 𝜕𝜇 2 − 𝒵 2 𝜕𝜇
2𝜋𝑚𝜏 3/2
For an ideal classical gas 𝑛 = 𝑒 𝜇 /𝜏 𝑉 2
𝜕𝑛 1 2𝜋𝑚𝜏 3/2
= 𝜏 𝑒 𝜇 /𝜏 𝑉
𝜕𝜇 2
𝜕𝑛 1
= 𝑛
𝜕𝜇 𝜏
𝜕𝑛
𝜏 𝜕𝜇 = 𝑛
2
𝛿𝑛 =𝑛
2
𝛿𝑛 = 𝑛
2
∆𝑛 𝛿𝑛 𝑛 1
= = =
𝑛 𝑛 𝑛 𝑛
𝑝𝑉 = 𝑛 𝑘𝑇
𝑝𝑉
𝑛 = 𝑘𝑇
∆𝑛 𝑘𝑇 1/2
=
𝑛 𝑝𝑉
Here Na , ma and Va represent the number of particles, the mass of each particle
and volume of system a.
Here Nb , mb and Vb represent the number of particles, the mass of each particle
and volume of system b.
Entropy is an extensive quantity and satisfy the additive property.
If the entropy given by equation (3) had satisfied the additive property , then by
removing partition and allowing the gas molecules to mix freely, the entropy
of the joint system would be
𝑆𝑎𝑏 = 𝑆𝑎 + 𝑆𝑏
3 3 3
= 𝑁𝑎 𝑘[𝑙𝑜𝑔 𝑉𝑎 + 𝑙𝑜𝑔 𝑚𝑎 + 𝑙𝑜𝑔𝑇 + 𝐶] + 𝑁𝑏 𝑘[𝑙𝑜𝑔 𝑉𝑏 + 𝑙𝑜𝑔 𝑚𝑏 +
2 2 2
3
𝑙𝑜𝑔𝑇 + 𝐶] ….(5)
2
If the particles of the two system are the same and for convenience we take
𝑉𝑎 = 𝑉𝑏 = 𝑉, 𝑁𝑎 = 𝑁𝑏 = 𝑁 and 𝑚𝑎 = 𝑚𝑏 = 𝑚 ,then the entropy of the
individual system be,
𝑆𝑎 = 𝑆𝑏
3 3
= 𝑁𝑘 𝑙𝑜𝑔 𝑉 + 𝑙𝑜𝑔 𝑚 + 𝑙𝑜𝑔𝑇 + 𝐶 ….(6)
2 2
𝑔𝑖𝑛𝑖
can be replaced by 𝑊= 𝑛𝑖!
𝑉 2𝜋𝑚𝐾𝑇 3/2 5
= 𝑁𝑘 𝑙𝑜𝑔 + 𝑁𝑘. ....(9)
𝑁 2 2
The entropy given by this equation satisfies the additive property since here in the
argument of logarithm we have V/N in place of V.
In equation (9) replacing N by 2N and V by 2V , the entropy of the combined
system be given by,
2𝑉 2𝜋𝑚𝐾𝑇 3/2 5
𝑆𝑎𝑏 = 2𝑁𝑘 𝑙𝑜𝑔 + 2𝑁𝑘 .
2𝑁 2 2
𝑉 2𝜋𝑚𝐾𝑇 3/2 5
= 2 𝑁𝑘 𝑙𝑜𝑔 + 2 𝑁𝑘
𝑁 2
= 2𝑆 = 𝑆𝑎 + 𝑆𝑏 ….(10)
Thus Gibb‟s paradox is resolved .The resolution of Gibb‟s paradox is an example
of the success of the quantum theory.
A pure classical state is represented by a single moving point in phase space, that
have definite value of coordinates q1 , q 2 ,...........q f and canonical momenta
p1 , p2 ,....... p f at each instant of time.
i cij j …(2)
j
The probability that a measurement of A will give an is then simply cin*cin.
But the probability that the assembly is in the state i is i . Therefore the
probability that the measurement of A will yield an is just
i cin * cin i ij cin * cin
i j
ij c in c in*
i j
^
c in c in * i * j dq
i j
cin * cin * j dq
i
= n * n dq …(3)
Therefore the probability that the measurement on A will give an is just
n n dq nn
* …(4)
Now suppose we wish to calculate the average value of A. This will be simply
A an nn
=
n
n * an n dq
n * A n n dp
Trace A .
1 𝑁 ∗
𝜌𝑚𝑛 = 𝑁 𝑖=1 𝑐𝑖𝑛 𝑐𝑖𝑚 ….(1)
∗
𝜕𝜌 𝑚𝑛 1 𝜕𝑐𝑖𝑛 𝜕𝑐 𝑖𝑚
Again =𝑁 𝑁
𝑖=1 𝜕𝑡
∗
𝑐𝑖𝑚 + 𝑐𝑖𝑚 ….(2)
𝜕𝑡 𝜕𝑡
ħ 𝜕
𝑖 2𝜋 𝜕𝑡 𝛹𝑖 = 𝐻 𝛹𝑖
ħ 𝜕
𝑖 𝑘 𝜕𝑡 𝑐𝑖𝑘 𝜙𝑘 = 𝑘 𝑐𝑖𝑘 𝐻 𝜙𝑘 ….(3)
2𝜋
ħ 𝜕
𝑖 2𝜋 𝑘 𝜕𝑡 𝑐𝑖𝑗 𝜙𝑗 ∗ 𝜙𝑘 𝑑𝑞 = 𝑘 𝑐𝑖𝑘 𝜙𝑗 ∗ 𝐻 𝜙𝑘 𝑑𝑞
ħ 𝜕
𝑖 2𝜋 𝜕𝑡 𝑐𝑖𝑗 = 𝑘 𝑐𝑖𝑘 𝐻𝑗𝑘 [𝐻𝑗𝑘 = 𝜙𝑗 ∗ 𝐻 𝜙𝑘 𝑑𝑞] ….(4)
ħ 𝜕
−𝑖 2𝜋 𝜕𝑡 𝑐𝑖𝑗 ∗ = 𝑘 𝑐𝑖𝑘 ∗ 𝐻𝑗𝑘 ∗ ….(5)
Substituting equations (4) and (5) in equation (2) and taking Hij=Hji*
ħ 𝜕𝜌 𝑚𝑛 ħ 1 𝑁 𝜕𝑐 𝑖𝑛 ∗ 𝜕𝑐 𝑖𝑚
𝑖 2𝜋 = 𝑖 2𝜋 𝑁 𝑖=1 𝑐𝑖𝑚 + 𝑐𝑖𝑛 ∗
𝜕𝑡 𝜕𝑡 𝜕𝑡
1 1
= −𝑁 𝑘 𝑐𝑖𝑘 ∗ 𝐻𝑛𝑘 ∗ 𝑐𝑖𝑚 + 𝑁 𝑘 𝑐𝑖𝑛 ∗ 𝑐𝑖𝑘 𝐻𝑚𝑘
= − 𝜌, 𝐻 𝑚𝑛
Manonmaniam Sundaranar University, Directorate of Distance & Continuing Education, Tirunelveli.
158
Thus the time development of density matrix is given by
ħ 𝜕𝜌
𝑖 2𝜋 𝜕𝑡 = − 𝜌, 𝐻
mn n *H m dq …(2)
If we consider ρ( H )as a power series in H ,then
H a0 a1 H a2 H 2 ........
as H mn Em mn
a0 a1Em a2 Em ..... mn
2
Thus the probability of observing the eigen value En in the range E and E+E is
proportional to 1 / .
(a) MICRO-CANONICAL ENSEMBLE
For a closed, isolated thermodynamic system i.e. a system with assigned values
for the independent variables E; n1, n2 ... nr ; x1, x2 ... xs, using energy eigen
functions as basic vectors, we write
mn mn n , …(7)
1
Where n for E E n E E ...(8)
= 0 otherwise,
where E is a very small range in E, and is the number of basic states n
belonging to eigen values En in the range E≤ En≤ E + E. The constant 1 /
result from the normalization
n
n 1 ...(9)
From equation (8) it is clear that for a system known to have an energy between E
and E + E, all basic states n belonging to an eigen value En in this range have
ρ constant e H ...(10)
1
Where = constant = . From equation (1) we can write
kT
1
or constant =
n
n * e H m dq
1
= where Z n * e H m dq …(12)
Z n
H
e
and …(13)
Z
where Z (, x1,… xs, n1…nr ) is called the partition function.
p,q
The classical partition function is e H dp dq …(14)
Comparing equations (12) and (14), the integration over in classical mechanics
is replaced by a summation over quantum states.
Z may be written in alternative ways
Z e H trace e H
n nn
1
n m e Em dq
*
...(15)
Z
1 2
H
Because e 1 H ( H )
2!
and H m E m m
Now we obtain
m dq n meEm dq
* H *
n e
1
= mn e En
Z
= mn n
Where n is the probability that a system, chosen at random from canonical
ensemble, will be found in the energy state En
The partition function is written as
Z n e H n dq
*
n eE n dq
2
e En ....(16)
n
mn constant e -n n e H m dq
*
…(2)
The total probability of finding a system in one or another state N will be taken as
normalized to unity. Thus we can write,
N ,n
nn 1 …(3)
1
= …(4)
Z
where Z is called the „grand partition function‟. Therefore
e H e n
…(5)
Z
From equation (5), we get
en
Z
* H
nn n e n dq
e n
Z(n )
Z
= (n) say. …(6)
Combining equations (5) and (6), we get
(n) e H ( n )
…(7)
Z ( n)
If the basis vectors are a set of energy eigen functions, then, the probability of
observing the state n (N, x) is
e n e En
n N nn …(8)
e n e En
(d) CONDITION FOR STATISTICL EQUILIBRIUM:
When a system is in equilibrium, its corresponding ensemble must be
stationary.i.e
In the energy representation, the basic function n are the eigen functions of the
Hamiltonian H .So, the matrices H and are diagonal. Thus
mn n mn …(2)
In the representation, the density operator may be written as
n n n dq H,
*
…(3)
n
kl k l dq k n dq n n l dq
* *
n
H ml l n ml H l n
l
0 H ml l n 0ml H l n
l
aa - -
- - aa
a a -
a - a
- a a
a a -
a - a
- a a
Thus in B-E statistics, there is a greater tendency for the particles to bunch
together in the same states in comparison to M-B statistics.
On the other hand, in the F-D statistics, there is a greater relative tendency to
particles to remain apart in different states than there is in classical statistics.
5.6 BOSE-EINSTEIN STATISTICS:
Conditions:
𝑔𝑖 𝑛𝑖 + 𝑔𝑖 − 1 !
𝑛𝑖 !
Similarly the permutations of cells among themselves will not give rise to
different arrangements. Hence equation (1) must also be divided by 𝑔𝑖 !
𝑔𝑖 𝑛𝑖 + 𝑔𝑖 − 1 ! 𝑛𝑖 + 𝑔𝑖 − 1 !
=
𝑛𝑖 ! 𝑔𝑖 ! 𝑛𝑖 ! 𝑔𝑖 − 1 !
𝐺
𝑛𝑖 + 𝑔𝑖 − 1 !
= . … (2)
𝑛𝑖 ! 𝑔𝑖 − 1 !
𝑛𝑖 + 𝑔𝑖 !
𝐺=
𝑛𝑖 ! 𝑔𝑖 !
The probability Ω of the system is proportional to the total number of eigen states.
𝛺 = 𝐺 × 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡
𝑛 𝑖 +𝑔 𝑖 !
𝛺= × 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡 . … (3)
𝑛 𝑖 !𝑔 𝑖 !
+ 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡
1 1
𝛿 log 𝛺 = 𝑛𝑖 + 𝑔𝑖 𝛿𝑛𝑖 + log 𝑛𝑖 + 𝑔𝑖 𝛿𝑛𝑖 − 𝑛𝑖 𝛿𝑛𝑖 −
𝑛 𝑖 +𝑔 𝑖 𝑛𝑖
log𝑛𝑖 𝛿𝑛𝑖+0
=
𝑛𝑖
− 𝑙𝑜𝑔 𝛿𝑛𝑖 . … (5)
𝑛 𝑖 +𝑔 𝑖
𝛿(log 𝛺) = 0.
𝑛𝑖
− 𝑙𝑜𝑔 𝛿𝑛𝑖 = 0
𝑛 𝑖 +𝑔 𝑖
𝑛𝑖
𝑙𝑜𝑔 𝛿𝑛𝑖 = 0 … . (6)
𝑛 𝑖 +𝑔 𝑖
𝑛𝑖 = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡
𝛿𝑛𝑖 = 0 ….(7)
𝜀𝑖 𝑛𝑖 = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡
𝜀𝑖 𝛿𝑛𝑖 = 0 . … (8)
Multiply the equation (7) by 𝛼 and (8) by β and then adding to equation (6),we
get,
𝑛𝑖
log + 𝛼 + 𝛽𝜀𝑖 𝛿𝑛𝑖 = 0 ….(9)
𝑛 𝑖 +𝑔 𝑖
𝛿𝑛𝑖 ≠ 0
𝑛𝑖
∴ log + 𝛼 + 𝛽𝜀𝑖 = 0
𝑛 𝑖 +𝑔 𝑖
𝑛𝑖
log + 𝛼 + 𝛽𝜀𝑖 = 0
𝑛 𝑖 +𝑔 𝑖
𝑛𝑖
log = −(𝛼 + 𝛽𝜀𝑖 )
𝑛 𝑖 +𝑔 𝑖
𝑛𝑖
= 𝑒 − 𝛼+𝛽 𝜀 𝑖
𝑛 𝑖 +𝑔 𝑖
𝑛 𝑖 +𝑔 𝑖
= 𝑒 𝛼+𝛽 𝜀 𝑖
𝑛𝑖
𝑔𝑖
1+ = 𝑒 𝛼+𝛽 𝜀 𝑖
𝑛𝑖
𝑔𝑖
𝑛𝑖 = 𝛼 +𝛽 𝜀 𝑖 −1 ….(10)
𝑒
This equation represents the most probable distribution for a system obeying Bose
- Einstein statistics and known as Bose - Einstein distribution law.
Conditions:
𝑔𝑖 !
𝑛𝑖 ! 𝑔𝑖 − 𝑛𝑖 !
Considering all the available groups such that 𝑛1 particles with energy 𝜀1 , 𝑛2
particles with energy 𝜀2 and so on. Then the total number of possible
arrangements is given by
𝑔𝑖 !
𝐺= . … (2)
𝑛 𝑖 ! 𝑔 𝑖 −𝑛 𝑖 !
𝛺 = 𝐺 × 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡
𝑔𝑖 !
𝛺= × 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡 . … (3)
𝑛 𝑖 ! 𝑔 𝑖 −𝑛 𝑖 !
+ 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡
1 𝑔𝑖 − 𝑛𝑖
𝛿(log 𝛺) = −𝑛𝑖 𝛿𝑛𝑖 − log 𝑛𝑖 𝛿𝑛𝑖 − −𝛿𝑛𝑖 − log 𝑔𝑖 − 𝑛𝑖 (−𝛿𝑛𝑖 )
𝑛𝑖 𝑔𝑖 − 𝑛𝑖
𝑛𝑖
=− log 𝛿𝑛𝑖 . … (5)
𝑔 𝑖 −𝑛 𝑖
𝛿(log 𝛺) = 0
𝑛𝑖
− log 𝛿𝑛𝑖 = 0.
𝑔 𝑖 −𝑛 𝑖
𝑛𝑖
log 𝛿𝑛𝑖 = 0 . … (6)
𝑔 𝑖 −𝑛 𝑖
𝑛𝑖 = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡
𝛿𝑛𝑖 = 0 . … (7)
𝜀𝑛𝑖 = 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡.
𝜀𝛿𝑛𝑖 = 0 . … (8)
Multiplying equation (7) by 𝛼 and (8) by 𝛽 and then adding to equation (6), we
get
𝑛𝑖
log + 𝛼 + 𝛽𝜀𝑖 𝛿𝑛𝑖 = 0
𝑔 𝑖 −𝑛 𝑖
𝑛𝑖
log + 𝛼 + 𝛽𝜀𝑖 = 0
𝑔 𝑖 −𝑛 𝑖
𝑛𝑖
log + 𝛼 + 𝛽𝜀𝑖 = 0
𝑔 𝑖 −𝑛 𝑖
𝑛𝑖
log = − 𝛼 + 𝛽𝜀𝑖
𝑔 𝑖 −𝑛 𝑖
𝑛𝑖
= 𝑒 − 𝛼+𝛽 𝜀 𝑖
𝑔 𝑖 −𝑛 𝑖
𝑔 𝑖 −𝑛 𝑖
= 𝑒 𝛼+𝛽 𝜀 𝑖
𝑛𝑖
𝑔𝑖
− 1 = 𝑒 𝛼+𝛽 𝜀 𝑖
𝑛𝑖
Manonmaniam Sundaranar University, Directorate of Distance & Continuing Education, Tirunelveli.
173
𝑔𝑖
= 𝑒 𝛼+𝛽𝜀 𝑖 + 1
𝑛𝑖
𝑔𝑖
𝑛𝑖 =
𝑒 𝛼+𝛽𝜀 𝑖 + 1
This equation represents the most probable distribution for a system obeying
Fermi-Dirac statistics and known as Fermi-Dirac distribution law.
5.8 AN IDEAL GAS IN A QUANTUM MECHANICAL MICRO CANONICAL
ENSEMBLE:
We consider a gaseous system of N non-interacting, indistinguishable particles
confined to a space of volume V and energy E.
Let 𝛺(𝑁, 𝑉, 𝐸) be the number of distinct microstates accessible to the system
under the macro state 𝛺(𝑁, 𝑉, 𝐸).
Let 𝜀𝑖 denote the average energy of a level and 𝑔𝑖 be the number of levels in the
𝑖 𝑡 cell.𝑔𝑖 ≫ 1.
We have n1 particles in the first cell, n 2 particles in the second cell and so on. The
distribution set 𝑛𝑖 must confirm to the conditions,
𝑖 𝑛𝑖 = 𝑁 ….(1)
𝑖 𝑛𝑖 𝜀𝑖 = 𝐸 ….(2)
Then 𝛺 𝑁, 𝑉, 𝐸 = 𝛴 ′ 𝑊{𝑛𝑖 } ….(3)
𝑊{𝑛𝑖 } is the number of distinct microstates associated with the distribution
set{𝑛𝑖 }. The primed summation goes over all distribution sets that confirm to
conditions (1) and (2),
𝑊{𝑛𝑖 } = 𝑖 𝑤(𝑖) ….(4)
𝑤(𝑖)is the number of distinct ways in which the ni identical and indistinguishable
In the F-D case, no single level can accommodate more than one particle.
𝑔𝑖 !
𝑤𝐹𝐷 𝑖 = 𝑛
𝑖! 𝑔 𝑖 −1 !
𝑔𝑖 !
and hence 𝑊𝐹𝐷 𝑛𝑖 = 𝑛 ....(6)
𝑖 ! 𝑔 𝑖 −1 !
Manonmaniam Sundaranar University, Directorate of Distance & Continuing Education, Tirunelveli.
174
In M-B case the particle are distinguishable.
𝑔𝑖 𝑛 𝑖
𝑊𝑀𝐵 𝑛𝑖 = ….(7)
𝑛𝑖!
Our condition for determining the most probable distribution set 𝑛𝑖∗ now turns
out to be,
𝛿𝑙𝑛𝑊 𝑛𝑖 − 𝛼 𝑖 𝛿𝑛𝑖 + 𝛽 𝑖 𝜀𝑖 𝛿𝑛𝑖 =0 ….(9)
𝑙𝑛𝑊 𝑛𝑖 = 𝑖 ln 𝑤(𝑖)
In B-E case,
(𝑛 𝑖 +𝑔 𝑖 −1)!
𝑊𝐵𝐸 𝑛𝑖 =
𝑛 𝑖 ! 𝑔 𝑖 −1 !
𝑛 𝑖 +𝑔 𝑖 !
≈ 𝑛 𝑖 !𝑔 𝑖 !
ln 𝑊𝐵𝐸 𝑛𝑖 = ln 𝑛𝑖 + 𝑔𝑖 ! − ln 𝑛𝑖 ! − ln 𝑔𝑖 !
= 𝑛𝑖 + 𝑔𝑖 𝑙𝑛 𝑛𝑖 + 𝑔𝑖 − 𝑛𝑖 + 𝑔𝑖 − 𝑛𝑖 ln 𝑛𝑖 + 𝑛𝑖 −
𝑔𝑖ln𝑔𝑖+𝑔𝑖
= 𝑛𝑖 + 𝑔𝑖 𝑙𝑛 𝑛𝑖 + 𝑔𝑖 − 𝑛𝑖 𝑙𝑛 𝑛𝑖 −
𝑔𝑖ln𝑔𝑖
= 𝑛𝑖 ln 𝑛𝑖 + 𝑔𝑖 − 𝑛𝑖 ln 𝑛𝑖 + 𝑔𝑖 ln 𝑛𝑖 + 𝑔𝑖 − 𝑔𝑖 ln 𝑔𝑖
𝑛 𝑖 +𝑔 𝑖 𝑛 𝑖 +𝑔 𝑖
= 𝑛𝑖 ln + 𝑔𝑖 ln
𝑛𝑖 𝑔𝑖
𝑔𝑖 𝑛𝑖
= 𝑛𝑖 ln + 1 + 𝑔𝑖 ln 1 − 𝑎
𝑛𝑖 𝑔𝑖
In general,
𝑔𝑖 𝑔𝑖 𝑛𝑖
ln 𝑊 𝑛𝑖 = 𝑛𝑖 ln −𝑎 − ln 1 − 𝑎 ….(10)
𝑛𝑖 𝑎 𝑔𝑖
a= -1 BE case
a=+1 FD case
a=0 MB case
ln𝑛𝑖𝛿𝑛𝑖
= ln 𝑛𝑖 + 𝑔𝑖 − ln 𝑛𝑖 𝛿𝑛𝑖
𝑛 𝑖 +𝑔 𝑖
= ln 𝛿𝑛𝑖
𝑛𝑖
𝑔𝑖
= ln + 1 𝛿𝑛𝑖
𝑛𝑖
In F-D case,
𝑔𝑖 !
𝑊𝐹𝐷 = 𝑛 𝑖 ! 𝑔 𝑖 −𝑛 𝑖 !
ln 𝑊𝐹𝐷 = ln 𝑔𝑖 ! − ln 𝑛𝑖 ! − ln 𝑔𝑖 − 𝑛𝑖 !
= 𝑔𝑖 ln 𝑔𝑖 − 𝑔𝑖 − 𝑛𝑖 𝑙𝑛 𝑛𝑖 + 𝑛𝑖 − 𝑔𝑖 − 𝑛𝑖 ln 𝑔𝑖 − 𝑛𝑖 +
𝑔𝑖−𝑛𝑖
= 𝑔𝑖 ln 𝑛𝑖 − 𝑛𝑖 ln 𝑛𝑖 − 𝑔𝑖 − 𝑛𝑖 ln(𝑔𝑖 − 𝑛𝑖 )
𝛿𝑛 𝑖 −𝛿𝑛 𝑖
𝛿 ln 𝑊𝐹𝐷 = −𝑛𝑖 − ln 𝑛𝑖 𝛿𝑛𝑖 − 𝑔𝑖 − 𝑛𝑖 − ln 𝑔𝑖 −
𝑛𝑖 𝑔 𝑖 −𝑛 𝑖
𝑛𝑖−𝛿𝑛𝑖
In M-B case
𝑔𝑖 𝑛 𝑖
𝑊𝑀𝐵 𝑛𝑖 = 𝑛𝑖!
ln 𝑊𝑀𝐵 = 𝑛𝑖 ln 𝑔𝑖 − ln 𝑛𝑖 !
= 𝑛𝑖 ln 𝑔𝑖 − 𝑛𝑖 𝑙𝑛 𝑛𝑖 + 𝑛𝑖
𝛿𝑛 𝑖
𝛿 ln 𝑊𝑀𝐵 = 𝑛𝑖 . 0 +ln 𝑔𝑖 𝛿𝑛𝑖 − 𝑛𝑖 − ln 𝑛𝑖 𝛿𝑛𝑖 + 𝛿𝑛𝑖
𝑛𝑖
𝑔𝑖 𝑒 𝛼 +𝛽 𝜀 𝑖
= 𝑖 𝑛𝑖 ∗ 𝛼 + 𝛽𝜀𝑖 − ln
𝑎 𝑒 𝛼 +𝛽 𝜀 𝑖 +𝑎
𝑔𝑖 1
= 𝑖 𝑛𝑖 ∗ 𝛼 + 𝛽𝜀𝑖 − ln
𝑎 1+𝑎𝑒 − 𝛼 +𝛽 𝜀 𝑖
𝑔𝑖
= 𝑛𝑖 ∗ 𝛼 + 𝛽𝜀𝑖 + ln 1 + 𝑎𝑒 −𝛼−𝛽𝜀 𝑖 ….(13)
𝑎
The first sum on RHS of (13) is αN while the second sum is βE. For the third
sum, we have
1 𝑆
𝑔𝑖 ln 1 + 𝑎𝑒 𝛼+𝛽 𝜀 𝑖 = 𝑘 − 𝛼𝑁 − 𝛽𝐸 ….(14)
𝑎
𝜇 1
𝛼 = − 𝑘𝑇 and 𝛽 = 𝑘𝑇
1 ′ 𝑁! 𝑛𝜀
= 𝑁! {𝑛 𝜀 } 𝜀 𝑒 −𝛽𝜀 ….(8)
𝜀 𝑛𝜀 !
The number of single particle states with energies lying between ε and ε+dε is,
𝑄1 𝑉, 𝑇 = 𝜀 𝑒 −𝛽𝜀
2𝜋𝑉 3/2 ∞ −𝛽𝜀 1/2
≈ 2𝑚 𝑒 𝜀 𝑑𝜀 ….(10)
3 0
2𝜋𝑉 3/2 1 𝜋
= 2𝑚
3 2 𝛽3
∞ ′ 𝑛𝜀
= 𝑁=0 {𝑛 𝜀 } 𝜀 𝑧𝑒 −𝛽𝜀 ….(14)
The double summation in (14) over the number 𝑛𝜀 constrained by a fixed value of
the total number N and then over all possible values of N, which is equivalent to a
summation over all possible values of 𝑛𝜀 independently of one another.
Hence,
𝑛0 𝑛1
𝐐 𝑧, 𝑉, 𝑇 = 𝑛 0 ,𝑛 1 ,…… 𝑧𝑒 −𝛽 𝜀 0 𝑧𝑒 −𝛽𝜀 1 …
𝑛0 𝑛1
= 𝑛0 𝑧𝑒 −𝛽𝜀 0 𝑛1 𝑧𝑒 −𝛽𝜀 1 ….(15)
In B-E case 𝑛𝜀 can be either 0 or 1 or 2 or ……
In F-D case 𝑛𝜀 can be only 0 or 1.
1
𝐐 𝑧, 𝑉, 𝑇 = 𝜀 (1−𝑧𝑒 −𝛽𝜀 ) in B-E case with 𝑧𝑒 −𝛽𝜀 <1.
=∓ 𝜀 (1 ∓ 𝑧𝑒 −𝛽𝜀 )
ln ….(17)
with 𝑔𝑖 = 1
𝑃𝑉 1
𝑞 𝑧, 𝑉, 𝑇 = 𝑘𝑇 = 𝑎 𝜀 (1 + 𝑎𝑧𝑒 −𝛽𝜀 )
ln ….(18)
𝜕𝑞
𝑁≡𝑧 𝜕𝑧 𝑉,𝑇
1 𝑎𝑒 −𝛽𝜀
= 𝑧𝑎 1+𝑎𝑧 𝑒 −𝛽𝜀
𝑧
= 𝑒 𝛽𝜀 +𝑎𝑧
1
= ….(20)
𝑧 −1 𝑒 𝛽𝜀 +𝑎
𝜕𝑞
𝐸≡− 𝜕𝛽 𝑧,𝑉
1 𝑎𝑧 𝑒 −𝛽𝜀 (−𝜀)
= −𝑎 1+𝑎𝑧 𝑒 −𝛽𝜀
𝜀
= ….(21)
𝑧 −1 𝑒 𝛽𝜀 +𝑎
At the same time, the mean occupation number 𝑛𝜀 of level ε turns out to be,
1 1 𝜕𝐐
𝑛𝜀 = 𝐐 − 𝛽 𝜕𝜀 𝑧,𝑇,𝑎𝑙𝑙 𝑜𝑡𝑒𝑟 𝜀
1 𝜕𝑞
= −𝛽 𝜕𝜀 𝑧,𝑇,𝑎𝑙𝑙 𝑜𝑡𝑒𝑟 𝜀
1 1 𝑎𝑧 𝑒 −𝛽𝜀 (−𝛽)
= −𝛽 𝑎 1+𝑎𝑧 𝑒 −𝛽𝜀
𝑧𝑒 −𝛽𝜀
= 1+𝑎𝑧 𝑒 −𝛽𝜀
1
= 𝑧 −1 𝑒 𝛽𝜀 +𝑎 ….(22)
Thus the mean value <n> and the most probable value n* of the occupation number n
of a single particle state are identical.
Finally in Maxwell-Boltzmann case (a= 0), the mean occupation number takes the
form
n MB
exp ( ) kT exp( kT ) ….(2)
We note here that the distinction between the quantum statistics and the classical
statistics becomes imperceptible when for all values of ε ,
exp ( ) kT>>1. ….(3)
Condition (4) implies that the probability of any of the nε being greater than unity
is quite negligible.
The distinction between the classical treatment and the quantum-mechanical
treatment then becomes rather insignificant.
Correspondingly for large values of ( ) kT the quantum curves 1 and 2
essentially merge into the classical curve 3.
Manonmaniam Sundaranar University, Directorate of Distance & Continuing Education, Tirunelveli.
181
Condition (3) also implies that μ, the chemical potential of the system must be
negative and large in magnitude. This means that the fugacity z exp( kT ) of
1 1
2
n 2
Q ….(6)
Q
z ,T ,all other
It follows that
1 2
n n
2
2
ln Q
z ,T ,all other
1
n ….(7)
z ,T
For the relative mean-square fluctuation we obtain
n2 n
2
1 1 1
z e ….(8)
n
2
n
The actual value of this quantity will depend on the statistics of the particles
because for a given particle density (N/V) and a given temperature T, the value of
z will be different for different statistics.
Equation (8) can be written in the form
n2 n
2
1
2
a ….(9)
n n
to vanish as n 1.
p (n) F .D ( ze ) n 1 ze 1
= 1 n for n0
n for n 1 ….(11)
On normalization we get
( ze ) n / n! ( n ) n
n
p ( n) M . B e ….(12)
exp( ze ) n!
Equation (12) is clearly a Poisson distribution for which the mean square
deviation of the variable is equal to the mean value itself.
It also resembles the distribution of the total particle number N in a grand
canonical ensemble consisting of ideal, classical systems.
Here we seen that the ratio p (n) p (n 1) varies inversely with n, which is a
“normal” statistical behaviour of uncorrelated events.
The distribution in the Bose-Einstein case is geometric with a common ratio
n n 1 .
This means that the probability of a state ε acquiring one more particle for itself is
independent of the number of particles already occupying the state.
In comparison with “normal” statistical behaviour, bosons exhibit a special
tendency of “bunching” together. That is a positive statistical correlation among
them. Fermions exhibit negative statistical correlation.
5.11 IDEAL BOSE-EINSTEIN GAS:
𝑔𝑖 −𝜇 1
𝑛𝑖 = here 𝛼 = and 𝛽 = 𝑘𝑇
𝑒 𝛼 +𝛽 𝜀 𝑖 −1 𝑘𝑇
𝑔𝑖
= 1 𝛽𝜀 where, 𝐴 = 𝑒 −𝛼
𝑒 𝑖 −1
𝐴
Since the particles in a box are normal size and the translational levels are closely
spaced and hence now the summation is replaced by integration.
The number of particles states 𝑔 𝑝 𝑑𝑝 between momentum p and p+dp is given
by
4𝜋𝑝 2 𝑑𝑝 4𝜋𝑉 𝑝 2 𝑑𝑝
= 𝑔𝑠 = 𝑔𝑠
3 /𝑉 3
𝑔𝑠 = degenaracy factor = 1
4𝜋𝑉 𝑝 2 𝑑𝑝
Now 𝑔 𝑝 𝑑𝑝 = ….(2)
3
4𝜋𝑉 𝑝 2 𝑑𝑝 1
= . … (3)
3 𝑒 𝛼 +𝛽 𝜀 𝑖 −1
𝑝2
𝜀 = 𝜀𝑖 = 2𝑚
𝑝2 = 2𝑚𝜀
2𝑝𝑑𝑝 = 2𝑚𝑑𝜀
𝑚
= 𝑑𝜀
2𝑚𝜀
𝑚 1/2
= 𝑑𝜀
2𝜀
4𝜋𝑉 𝑚 1/2 1
𝑛 𝜀 𝑑𝜀 = (2𝑚𝜀) 𝑑𝜀 𝑒 𝛼 +𝜀/𝑘𝑇 −1
3 2𝜀
4𝜋𝑚𝑉 1/2 𝑑𝜀
= 2𝑚𝜀 . … (4)
3 𝑒 𝛼 +𝜀/𝑘𝑇 −1
𝜀
=𝑥
𝑘𝑇
𝜀 = 𝑘𝑇𝑥
𝑑𝜀 = 𝑘𝑇𝑑𝑥
𝑉 1/2 𝑑𝑥
3/2 2 𝑥
= 3 2𝜋𝑚𝑘𝑇 . … (5)
𝜋 𝑒 +𝑥 −1
𝛼
∞ 𝑉 1/2 𝑑𝑥
3/2 2 𝑥
𝑛= 0 3
2𝜋𝑚𝑘𝑇 𝜋 𝑒 +𝑥 −1
𝛼
1/2 𝑑𝑥
𝑉 3/2 2 ∞ 𝑥
= 3 2𝜋𝑚𝑘𝑇 𝜋 0 𝑒 +𝑥 −1
𝛼
𝑉 3/2
𝑛 = 3 2𝜋𝑚𝑘𝑇 𝑓1 𝛼 . … (6)
2 ∞ 𝑥 1/2 𝑑𝑥
where 𝑓1 𝛼 = . … (7)
𝜋 0 𝑒 𝛼 +𝑥 −1
∞ 𝑉 1/2 𝑑𝑥
3/2 2 𝑥
= 𝑘𝑇𝑥 3 2𝜋𝑚𝑘𝑇
0 𝜋 𝑒 +𝑥 −1
𝛼
3 𝑉 3/2 4 ∞ 𝑥 3/2 𝑑𝑥
= 2 3 𝑘𝑇 2𝜋𝑚𝑘𝑇 3 𝜋 0 𝑒 𝛼 +𝑥 −1
3 𝑉 3/2
= 2 3 𝑘𝑇 2𝜋𝑚𝑘𝑇 𝑓2 𝛼 . … (8)
4 ∞ 𝑥 3/2 𝑑𝑥
𝑓2 𝛼 = 3 0 𝑒 𝛼 +𝑥 −1
. … (9)
𝜋
2 ∞ 𝑥 1/2 𝑑𝑥
𝑓1 𝛼 = 𝜋 0 𝑒 𝛼 +𝑥 −1
2 ∞ 𝑥 1/2 𝑑𝑥
= 𝑥
𝜋 0 𝑒 −1
𝐴
2 ∞ 1/2 𝑒𝑥 −1
= 𝑥 −1 𝑑𝑥
𝜋 0 𝐴
2 ∞ 1/2
= 𝑥 𝐴𝑒 −𝑥 1 − 𝐴𝑒 −𝑥 −1
𝑑𝑥
𝜋 0
2 ∞ 1/2
= 𝑥 𝐴𝑒 −𝑥 1 + 𝐴𝑒 −𝑥 + 𝐴2 𝑒 −2𝑥 + ⋯ 𝑑𝑥
𝜋 0
2 𝜋 𝐴2 𝐴3
= 𝐴+ + + ⋯..
𝜋 2 23/2 33/2
𝐴2 𝐴3
= 𝐴 + 23/2 + 33/2 + ⋯ … … …
∞ 𝐴𝑟
𝑓1 𝛼 = 𝑟=1 𝑟 3/2 . … (10)
4 ∞ 𝑥 3/2 𝑑𝑥
Similarly 𝑓2 𝛼 = 3 −(𝛼 +𝑥 ) −1
𝜋 0 𝑒
𝐴2 𝐴3
= 𝐴 + 25/2 + 35/2 + ⋯
∞ 𝐴𝑟
𝑓2 𝛼 = 𝑟=1 𝑟 5/2 . … (11)
𝑉 3/2 𝐴2 𝐴3
= 3 2𝜋𝑚𝑘𝑇 𝐴 + 23/2 + 33/2 + ⋯ . … (12)
3 𝑉 3/2
𝐸 = 2 𝑘𝑇 2𝜋𝑚𝑘𝑇 𝑓2 (𝛼).
3
3 𝑉 3/2 𝐴2 𝐴3
𝐸 = 2 𝑘𝑇 2𝜋𝑚𝑘𝑇 𝐴 + 25/2 + 35/2 + ⋯ . … 13
3
−1
𝐸 3 𝐴2 𝐴3 𝐴2 𝐴3
= 2 𝑘𝑇 𝐴 + 25/2 + 35/2 + ⋯ 𝐴 + 23/2 + 33/2 + ⋯
𝑛
−1
3 𝐴 𝐴2 𝐴 𝐴2
= 2 𝑘𝑇 1 + 25/2 + 35/2 + ⋯ 1 + 23/2 + 33/2 + ⋯
3 𝐴 𝐴2
= 2 𝑘𝑇 1 − 25/2 + 35/2 + ⋯
3 𝐴 𝐴2
𝐸 = 2 𝑛𝑘𝑇 1 − 25/2 + 35/2 + ⋯ . … (14)
𝑉 3/2
𝑛 = 3 2𝜋𝑚𝑘𝑇 𝑓1 (𝛼)
𝑛 3
𝑓1 𝛼 = 𝑉 2𝜋𝑚𝑘𝑇 3/2
𝑛 3
𝐴=𝑉 . … (15)
2𝜋𝑚𝑘𝑇 3/2
B distribution.
For 𝐴 << 1, 𝑓1 𝛼 = 𝑓2 𝛼 = 𝐴
Manonmaniam Sundaranar University, Directorate of Distance & Continuing Education, Tirunelveli.
187
𝑉 3/2
Equation (12) becomes 𝑛 = 3 2𝜋𝑚𝑘𝑇 ×𝐴
3 𝑉 3/2
Equation (13) becomes 𝐸 = 2 𝑘𝑇 2𝜋𝑚𝑘𝑇 ×𝐴
3
𝐸 3 3
= 2 𝑘𝑇 ⇒ 𝐸 = 2 𝑛𝑘𝑇
𝑛
Thus the degree of degeneracy will be large when T is low, n/V is large and m is
small.
For low energy values A=1 and α=0.
1 1
Then 𝑓1 (𝛼) 𝑚𝑎𝑥 = 𝑓1 0 = 1 + 23/2 + 33/2 + ⋯
= 2.612
𝑛 2𝜋𝜋𝑚𝑘𝑇 3/2
Now = 2.612 . … (4)
𝑉 𝑚𝑎𝑥 3
𝑛 2𝜋𝑚𝑘𝑇 3/2
> 2.612 . … (5)
𝑉 3
Manonmaniam Sundaranar University, Directorate of Distance & Continuing Education, Tirunelveli.
188
n/v can be alternatively expressed in terms of critical temperature 𝑇𝑜 defined as
𝑛 2𝜋𝑚𝑘 𝑇0 3/2
= 2.612 . … (6)
𝑉 3
3/2 𝑛 1
2𝜋𝑚𝑘𝑇0 = 3 𝑉 2.612
𝑛 1 2/3
2𝜋𝑚𝑘𝑇0 = 2 𝑉 2.612
2 𝑛 1 2/3
𝑇0 = 2𝜋𝑚𝑘 𝑉 2.612
𝑔 𝜀 𝑑𝜀
𝑛 𝜀 𝑑𝜀 = 𝑒 𝛼 +𝜀/𝑘𝑇 −1 . … (8)
4𝜋𝑚𝑣 1/2
And 𝑔 𝜀 𝑑𝜀 = 2𝑚𝜀 𝑑𝜀 . … (9)
3
𝜀𝑖 = 𝜀0 = 0
𝑔𝑖 = 𝑔 𝜀𝑖 = 1
Manonmaniam Sundaranar University, Directorate of Distance & Continuing Education, Tirunelveli.
189
1
𝑛 = 𝑛0 = 𝑒 𝛼 −1 . … (10)
= 𝑛0 + 𝑛′ . … (11)
4𝜋𝑚𝑉 ∞ 𝜀 1/2 𝑑𝜀
𝑛′ = 2𝑚 1/2
….(12)
3 0 𝑒 𝛼 +𝜀/𝑘𝑇 −1
𝑉
𝑛′ = 3 2𝜋𝑚𝑘𝑇 3/2
𝑓1 𝛼 . . . . (13)
𝑛 2𝜋𝑚𝑘 𝑇0 3/2
From equation (6) = 2.612
𝑉 3
𝑉 3/2 1
2𝜋𝑚𝑘𝑇0 =𝑛
3 2.612
𝑇 3/2 𝑓1 (𝛼)
Now 𝑛′ = 𝑛 .…(14)
𝑇0 2.612
𝑛 = 𝑛0 + 𝑛′
𝑇 3/2
𝑛 = 𝑛0 + 𝑛 𝑇0
𝑇 3/2
𝑛𝑜 = 𝑛 − 𝑛 𝑇0
𝑇 3/2
𝑛0 = 𝑛 1 − 𝑓𝑜𝑟 𝑇 < 𝑇0 . … (15)
𝑇0
The frequencies of the trap are controlled by the applied magnetic field. One
can then lower the trap barrier using a resonant transition to remove the highest
energy atoms in the trap. If the atoms in the vapor are sufficiently coupled to one
other, then the remaining atoms in the trap are cooled by evaporation.
If the interactions between the atoms in the gas can be neglected, the energy of
each atom in the harmonic oscillator potential is
where 0 12 3 1/ 3 ; this assumes a single spin state per atom. The
, T
kT
4
x ln 1 e e dx
kT
4
g z
20 0
2 x
… (4)
3
0 3 4
where z = exp () is the fugacity. Volume is not a parameter in the
thermodynamic potential since the atoms are confined by the harmonic trap. The
average number of atoms in the excited states in the trap is
kT
N T g 3 z .… (5)
T 0
where (3) = g3(1) 1.202. While the spacing of the energy levels is of order 0 ,
the critical temperature for condensation is much larger than the energy spacing of
the lowest levels for N>> 1.
For T <Tc the number of atoms in the excited states is
3 3
N excited (3) kT T
.… (7)
N N 0 Tc
so the fraction of atoms that condense into the ground state of the harmonic
oscillator is
3
N0 T
1 ; .… (8)
N Tc
(a) Detection of the Bose-Einstein condensate:
The linear size of the ground state wave function in Cartesian direction is
a , .…(9)
m
while the linear size of the thermal distribution of the noncondensed atoms in that
direction is
kT kT
athermal a .…(10)
m2
At trap frequency f= 100 Hz and temperature T =100 nK, these sizes are about 1
m and 5m, respectively.
Instead of measuring the atoms directly in the trapping potential, experimenters
usually measure the momentum distribution of the ultracold gas by a time-of-
flight experiment.
At time t = 0, the magnetic field is turned off suddenly, eliminating the trapping
potential.
N0 3 r2
n0 r , t N o 0 r , t
1
.…(11)
2
3 / 2 1 a 1 2 t 2
exp
a 2 1 2 t 2
The atoms that are not condensed into the ground state can be treated semi
classically, that is, the position-momentum distribution function is treated
classically while the density follows the Bose-Einstein distribution function:
1
f (r , p,0) ….(12)
p m 2 2
2
exp 1 x 22 y 2 32 z 2 u 1
2m 2
After the potential is turned off at t = 0, the distribution evolves ballistically:
pt
f (r , p, 0) f r , p, 0 ….(13)
m
The spatial number density of atoms in the excited states is
1
nexcited r , t
pt
3
f r , p, t dp ….(14)
h m
which can be integrated to give
where = h / 2mkT is the thermal deBroglie wavelength. The integrals over the
condensed state and the excited states correctly count all the atoms:
N 0 no r , t dr , ….(16a)
Note that at early times t 1 both the condensed and the excited distributions
are anisotropic due to the anisotropic trapping potential.
However, at late times t 1 , the atoms from the excited states form a
spherically symmetric cloud because of the isotropic momentum dependence of
the t = 0 distribution function.
By contrast, the atoms that were condensed into the ground state expand
anisotropically due to the different spatial extents of the ground state
wavefunction at t = 0.
The direction that has the largest 0 is quantum mechanically squeezed the most at
t = 0; so, according to the uncertainty principle, it expands the fastest.
This is an important feature of the experimental data that confirms the onset of
Bose-Einstein condensation.
(b) Thermodynamic properties of the Bose-Einstein condensate :
The temperature, condensate fraction, and internal energy can all be observed
using time- of-flight measurements. The internal energy can also be written in
terms of the function gv(z):
U , T
3 1
d 3
kT g z
4
….(17)
0 20 0 3 4
3
e 1
U U U
C N T
T N T T T N
U N
U T T
T N
T ….(18)
Manonmaniam Sundaranar University, Directorate of Distance & Continuing Education, Tirunelveli.
195
Equations (5) and (6) can be used to determine the fugacity z numerically, as
shown in Figure (a). The fugacity can then be used in equation (17) to obtain the
scaled internal energy. [Figure (b)]
T 4 4
3 for T Tc ,
T 3
3 4
U
NkTc T g z
3 T 3 for T Tc;
4
c ….(19)
The scaled specific heat is given by
124 T 3
for T Tc ,
C N 3 Tc
Nk 1 T 9 g 2 z
12 g 4 z 3 for T Tc ,
3 Tc g 2 z
.(20)
and shown in figure (c).
CN 3 10.805 as T Tc
Nk 124 9 3
4.228 as T Tc
𝑔𝑖 −𝜇 1
𝑛𝑖 = here 𝛼 = and 𝛽 = 𝑘𝑇
𝑒 𝛼 +𝛽 𝜀 𝑖 +1 𝑘𝑇
𝑔𝑖
= 1 𝛽𝜀 where, 𝐴 = 𝑒 −𝛼
𝑒 𝑖 +1
𝐴
Since the particles in a box are normal size and the translational levels are closely
spaced and hence now the summation is replaced by integration.
The number of particle states 𝑔 𝑝 𝑑𝑝 between momentum p and p+dp is given by
4𝜋𝑝 2 𝑑𝑝 4𝜋𝑉 𝑝 2 𝑑𝑝
= 𝑔𝑠 = 𝑔𝑠
3 /𝑉 3
4𝜋𝑉 𝑝 2 𝑑𝑝
Now 𝑔 𝑝 𝑑𝑝 = 𝑔𝑠 ….(2)
3
𝑝2
𝜀 = 𝜀𝑖 = 2𝑚
𝑝2 = 2𝑚𝜀
2𝑝𝑑𝑝 = 2𝑚𝑑𝜀
𝑚
𝑑𝑝 = 𝑑𝜀
𝑝
𝑚
= 𝑑𝜀
2𝑚𝜀
𝑚 1/2
= 𝑑𝜀
2𝜀
4𝜋𝑉 𝑚 1/2 1
𝑛 𝜀 𝑑𝜀 = 𝑔𝑠 (2𝑚𝜀) 𝑑𝜀
3 2𝜀 𝑒 𝛼 +𝜀/𝑘𝑇 +1
4𝜋𝑚𝑉 1/2 𝑑𝜀
= 𝑔𝑠 2𝑚𝜀 . … (4)
3 𝑒 𝛼 +𝜀/𝑘𝑇 +1
𝜀
=𝑥
𝑘𝑇
𝜀 = 𝑘𝑇𝑥
𝑑𝜀 = 𝑘𝑇𝑑𝑥
𝑉 1/2 𝑑𝑥
3/2 2 𝑥
= 𝑔𝑠 3 2𝜋𝑚𝑘𝑇 . … (5)
𝜋 𝑒 𝛼 +𝑥 +1
∞ 𝑉 1/2 𝑑𝑥
3/2 2 𝑥
𝑛= 𝑔
0 𝑠 3
2𝜋𝑚𝑘𝑇 𝜋 𝑒 +𝑥 +1
𝛼
1/2 𝑑𝑥
𝑉 3/2 2 ∞ 𝑥
= 𝑔𝑠 3 2𝜋𝑚𝑘𝑇 𝜋 0 𝑒 𝛼 +𝑥 +1
𝑉 3/2
𝑛 = 𝑔𝑠 3 2𝜋𝑚𝑘𝑇 𝑓1 𝛼 . … (6)
∞ 𝑉 1/2 𝑑𝑥
3/2 2 𝑥
= 𝑘𝑇𝑥 𝑔𝑠 3 2𝜋𝑚𝑘𝑇
0 𝜋 𝑒 𝛼 +𝑥 +1
𝑉 3/2 2 ∞ 𝑥 3/2 𝑑𝑥
= 𝑔𝑠 3
𝑘𝑇 2𝜋𝑚𝑘𝑇 0 𝑒 𝛼 +𝑥 +1
𝜋
3 𝑉 3/2 4 ∞ 𝑥 3/2 𝑑𝑥
= 2 𝑔𝑠 3 𝑘𝑇 2𝜋𝑚𝑘𝑇 3 𝜋 0 𝑒 𝛼 +𝑥 +1
3 𝑉 3/2
= 2 𝑔𝑠 3 𝑘𝑇 2𝜋𝑚𝑘𝑇 𝑓2 𝛼 . … (8)
4 ∞ 𝑥 3/2 𝑑𝑥
𝑓2 𝛼 = 3 0 𝑒 𝛼 +𝑥 +1
. … (9)
𝜋
𝑛 𝜀
𝑓 𝜀 =𝑔 𝜀
1
= 𝑒 𝛼 +𝜀/𝑘𝑇 +1
1
= 𝑒 (𝜀−𝜀 𝐹 )/𝑘𝑇 +1 . … (10)
𝜇 𝜀
Where 𝛼 = − 𝑘𝑇 = − 𝑘𝑇𝐹
𝑛 𝜀 𝑑𝜀 = 𝑓 𝜀 𝑔 𝜀 𝑑𝜀
Manonmaniam Sundaranar University, Directorate of Distance & Continuing Education, Tirunelveli.
199
4𝜋𝑚𝑉 1/2 𝜀 1/2 𝑑𝜀
= 𝑔𝑠 2𝑚 (𝜀−𝜀 𝐹)/𝑘𝑇 . . . . (11)
3 𝑒 +1
At absolute zero, all states with 0 < 𝜀 < 𝜀𝐹 (0) are completely filled and all states
with 𝜀 > 𝜀𝐹 (0) are empty.
𝜀𝐹 at T=0 K i.e. 𝜀𝐹 (0) is determined by,
𝜀 𝐹 (0)
𝑛= 0
𝑛 𝜀 𝑑𝜀
𝜀 𝐹 (0)
= 0
𝑓(𝜀)𝑔 𝜀 𝑑𝜀
𝜀 𝐹 (0)
= 0
𝑔 𝜀 𝑑𝜀
3/2 3𝑛 3
𝜀𝐹 0 = 4𝜋𝑉 𝑔
𝑠 2𝑚 3/2
2 3𝑛 2/3
𝜀𝐹 0 = 2𝑚 … . (12)
4𝜋𝑉 𝑔𝑠
𝜀 𝐹 (0)
𝑇𝐹 = 𝑘
2 3𝑛 2/3
= 2𝑚𝑘 4𝜋𝑉 𝑔𝑠
𝑚𝑛 𝑛 𝜌
𝜌= ⇒𝑉=𝑚
𝑉
2 3𝜌 2/3
𝑇𝐹 = 2𝑚𝑘 4𝜋𝑚 𝑔𝑠
2 3𝜌 2/3
𝑇𝐹 = 2𝑚 5/3 𝑘 … . (13)
4𝜋𝑔𝑠
Manonmaniam Sundaranar University, Directorate of Distance & Continuing Education, Tirunelveli.
200
𝑇 << 𝑇𝐹 strong degeneracy α is negative
𝑇 > 𝑇𝐹 weak degeneracy α is positive
𝑇 >> 𝑇𝐹 non degeneracy α is positive
(A) Degeneracy:
i. Weak degeneracy: (High temperature and low density)
𝑇 > 𝑇𝐹 ; 𝜀𝐹 is negative ; α is positive and hence 𝐴 < 1.
For 𝐴 < 1, we can write
1 1 −1
1 = 𝑒𝑥 + 1
𝑒 𝑥 +1 𝐴
𝐴
= 𝐴𝑒 −𝑥 1 + 𝐴𝑒 −𝑥 −1
= 𝐴𝑒 −𝑥 1 − 𝐴𝑒 −𝑥 + 𝐴2 𝑒 −2𝑥 + ⋯
2 ∞ 𝑥 1/2 𝑑𝑥
𝑓1 𝛼 = 𝜋 0 𝑒 𝛼 +𝑥 +1
2 ∞ 𝑥 1/2 𝑑𝑥
= 𝜋 0 1 𝑒 𝑥 +1
𝐴
2 ∞ 1/2
= 𝑥 𝑑𝑥 𝐴𝑒 −𝑥 1 − 𝐴𝑒 −𝑥 + 𝐴2 𝑒 −2𝑥 + ⋯
𝜋 0
𝐴2 𝐴3
= 𝐴 − 25/2 + 35/2 + ⋯ . … (2)
𝑉 3/2
Now 𝑛 = 𝑔𝑠 3 2𝜋𝑚𝑘𝑇 𝑓1 𝛼
𝑉 3/2 𝐴2 𝐴3
= 𝑔𝑠 3 2𝜋𝑚𝑘𝑇 𝐴 − 23/2 + 33/2 + ⋯ … . (3)
3 𝑣 3/2
𝐸 = 2 𝑔𝑠 3 2𝜋𝑚𝑘𝑇 𝑘𝑇 𝑓2 𝛼
3 𝑣 3/2 𝐴2 𝐴3
= 2 𝑔𝑠 3 2𝜋𝑚𝑘𝑇 𝐴 − 25/2 + 35/2 + ⋯ … . (4)
𝐸 3 𝐴2 𝐴3 𝐴2 𝐴3 −1
= 2 𝑘𝑇 𝐴 − 25/2 + 35/2 + ⋯ 𝐴 − 23/2 + 33/2 + ⋯
𝑛
3 𝐴 𝐴2
= 2 𝑘𝑇 1 + 25/2 − 35/2 + ⋯ … . (5)
1
= 1 𝜀−𝜀 /𝑘𝑇
𝑒 𝐹 +1
𝐴
=1 𝑓𝑜𝑟 0 ≤ 𝜀 ≤ 𝜀𝐹 (0)
=0 𝑓𝑜𝑟 𝜀 > 𝜀𝐹
The total internal energy of Fermi-Dirac gas at T=0 K gives zero point energy.
𝜀 𝐹 (0)
𝐸0 = 0
𝜀 𝑛 𝜀 𝑑𝜀
4𝜋𝑚𝑉 5/2
1/2 𝜀 𝐹 (0)
= 𝑔𝑠 2𝑚
3 5/2
3𝑛 2 3𝑛 2/3
= 5 2𝑚 4𝜋𝑉 𝑔𝑠
3𝑛
𝐸0 = 𝜀𝐹 0 … . (7)
5
2 𝐸0
Zero point pressure 𝑝0 = 3 𝑉
2 1 3𝑛
=3 𝜀𝐹 (0)
𝑉 5
2 𝑛
𝑝0 = 5 𝜀𝐹 0 … . (8)
𝑉
Thus a strongly degenerate Fermi-Dirac gas possesses energy and pressure even
at absolute zero.
(b) Above absolute zero: 𝑨 >> 1 𝑎𝑛𝑑 𝑇 << 𝑻𝑭
The Fermi gas is strongly degenerate and 𝜀 is still positive.
The number of particles in the energy range between 𝜀 and 𝜀 + 𝑑𝜀 is given by,
1
1
4𝜋𝑚𝑉 𝜀 2 𝑑𝜀
𝑛 𝜀 𝑑𝜀 = 𝑔𝑠 2𝑚 2 … . (9)
3 𝑒 𝜀−𝜀 𝐹 /𝑘𝑇 +1
1/2
4𝜋𝑉 𝑔𝑠 2𝑚 3 𝜀 1/2 𝑑𝜀
= 3𝑛
3𝑛 6 𝑒 𝜀−𝜀 𝐹 /𝑘𝑇 +1
2 3𝑛 2/3
𝜀𝐹 0 = 2𝑚 4𝜋𝑣 𝑔𝑠
3𝑛 2/3 2𝑚 𝜀 𝐹 (0)
=
4𝜋𝑣 𝑔𝑠 2
3𝑛 2𝑚 𝜀 𝐹 (0) 3/2
=
4𝜋𝑣 𝑔𝑠 2
3/2 1/2
2 2𝑚 3 𝜀 1/2 𝑑𝜀
𝑛 𝜀 𝑑𝜀 = 3𝑛 2𝑚 𝜀 𝐹 (0) 6 𝑒 (𝜀−𝜀 𝐹 )/𝑘𝑇 +1
3 𝑛 𝜀 1/2 𝑑𝜀
= … . (10)
2 𝜀 𝐹 (0) 3/2 𝑒 (𝜀−𝜀 𝐹 )/𝑘𝑇 +1
∞
𝑛= 0
𝑛 𝜀 𝑑𝜀
3 𝑛 ∞ 𝜀 1/2 𝑑𝜀
=2 … . (11)
𝜀 𝐹 (0) 3/2 0 𝑒 (𝜀−𝜀 𝐹 )/𝑘𝑇 +1
∞
𝐸= 0
𝜀 𝑛 𝜀 𝑑𝜀
3 𝑛 ∞ 𝜀 1/2 𝑑𝜀
= 3/2 0 𝑒 (𝜀−𝜀 𝐹 )/𝑘𝑇 +1
… . (12)
2 𝜀 𝐹 (0)
To solve the integrals in equation (11) and (12), we consider the general integral
∞ 𝜑 𝜀 𝑑𝜀
𝐼= 0 𝑒 (𝜀−𝜀 𝐹 )/𝑘𝑇 +1
… . (13)
∞ 𝜑 𝜀 𝑑𝜀 𝜀𝐹 𝜋𝑘𝑇 2 7
0 𝑒 (𝜀−𝜀 𝐹 )/𝑘𝑇 +1
= 0
𝜑 𝜀 𝑑𝜀 + 𝜑′ 𝜀=𝜀 𝐹 + 360 𝜋𝑘𝑇 4
𝜑′′′ 𝜀=𝜀 𝐹 +⋯
6
Take 𝜑 𝜀 = 𝜀 1/2
1
𝜑′ 𝜀=𝜀 𝐹
= 2 𝜀𝐹 −1/2
1
𝜑 ′′ 𝜀=𝜀 𝐹
= − 4 𝜀𝐹 −3/2
3
𝜑 ′′′ 𝜀=𝜀 𝐹
= 8 𝜀𝐹 −5/2
3 𝑛 ∞ 𝜀 1/2 𝑑𝜀
𝑛=2 3/2 0 𝑒 (𝜀−𝜀 𝐹 )/𝑘𝑇 +1
𝜀 𝐹 (0)
3 𝑛 2 𝜋𝑘𝑇 2 1 7 43
=2 3/2
𝜀𝐹 3/2 + 𝜀𝐹 −1/2 + 360 𝜋𝑘𝑇 𝜀𝐹 −5/2 + ⋯
𝜀 𝐹 (0) 3 6 2 8
−2/3
𝜀𝐹 1 𝜋𝑘𝑇 2
= 1+8
𝜀 𝐹 (0) 𝜀𝐹
1 𝜋𝑘𝑇 2
= 1 − 12 … . (14)
𝜀𝐹
1 𝜋𝑘𝑇 2
𝜀𝐹 = 𝜀𝐹 0 1 − 12 … . (15)
𝜀𝐹
1 𝜋𝑘𝑇 2
𝜀𝐹 = 𝜀𝐹 0 1 − 12 … . (16)
𝜀𝐹 0
Take 𝜑 𝜀 = 𝜀 3/2
𝜀𝐹 2
0
𝜀 3/2 𝑑𝜀 = 5 𝜀𝐹 5/2
3
𝜑 ′ (𝜀) 𝜀=𝜀 𝐹
= 2 𝜀𝐹 1/2
3
𝜑 ′′ (𝜀) 𝜀=𝜀 𝐹
= 4 𝜀𝐹 −1/2
3
𝜑 ′′′ (𝜀) 𝜀=𝜀 𝐹
= − 8 𝜀𝐹 −3/2
3 𝑛 ∞ 𝜀 3/2 𝑑𝜀
𝐸=2 3/2 0 𝑒 (𝜀−𝜀 𝐹 )/𝑘𝑇 +1
𝜀 𝐹 (0)
3 𝑛 2 𝜋𝑘𝑇 2 3 7 43
= 𝜀𝐹 5/2 + 𝜀𝐹 1/2 − 𝜋𝑘𝑇 𝜀𝐹 −3/2 + ⋯
2 𝜀 𝐹 (0) 3/2 5 6 2 360 8
3 𝜀𝐹 5/2 5 𝜋𝑘𝑇 2
𝐸 = 5 𝑛𝜀𝐹 0 1+8 … . (17)
𝜀 𝐹 (0) 𝜀𝐹
3 𝜀𝐹 5/2 5 𝜋𝑘𝑇 2
𝐸 = 5 𝑛𝜀𝐹 0 1+8 … . (18)
𝜀 𝐹 (0) 𝜀𝐹 0
3 5 𝜋𝑘𝑇 2 5 𝜋𝑘𝑇 2
= 5 𝑛𝜀𝐹 0 1 − 24 1+8
𝜀 𝐹 (0) 𝜀 𝐹 (0)
3 5 𝜋𝑘𝑇 2 5 𝜋𝑘𝑇 2
= 5 𝑛𝜀𝐹 0 1 − 24 +8
𝜀 𝐹 (0) 𝜀 𝐹 (0)
3 5 𝜋𝑘𝑇 2
𝐸 = 5 𝑛𝜀𝐹 0 1 + 12 … . 19
𝜀𝐹 0
2 𝐸
𝑝= 3 𝑉
2 𝑛𝜀 𝐹 0 5 𝜋𝑘𝑇 2
=5 1 + 12 (20)
𝑉 𝜀𝐹 0
Equations (19) and (20) give the approximate energy and pressure of a strongly
degenerate Fermi gas and also known as equation of state of an ideal Fermi gas.
Prepared by
Dr. L. RAVISANKAR
Assistant Professor of Physics
Sri KGS Arts College, Srivaikuntam – 628 619