Problem 1. Ultra-relativistic bosons in 2D (Points: 10+10+10 = 30)
In this problem, we consider non-interacting ultra-relativistic bosons in two dimensions (2D) in a 2D “volume” V, in contact with an external particle reservoir, and in thermal equilibrium with the surroundings. The system has a grand canonical partition function Zg given by
Zg =Y
k
1−e−β(εk−µ)−1 =eβpV.
Here, β = 1/kBT, wherekB is Boltzmanns’s constant andT is temperature, while µis the chemical potential of the system. Furthermore, pis the pressure of the system. For ultra-relativistic particles, the energy of a particle in a single-particle state specified by specifying k= (kx, ky), is given by
εk= ¯hck.
where cis the speed of light, ¯h=h/2π, his Planck’s constant, and k ≡ |k|. The Hamiltonian of the system is given by H =Pkεknk, with nk = 0,1,2, ....
a. Show that the pressure of the system and the average number of particles in the system,hNi, are given by
βpV = −
Z ∞
−∞dε g(ε) ln1−e−β(ε−µ), hNi = N0+
Z ∞
−∞
dε g(ε) eβ(ε−µ)−1,
where g(ε) = Pkδ(ε−εk) and δ(x) is the δ-function. Here, N0 is the number of particles in the lowest-energy state. Show also that the internal energy U =hHi is given by
U =
Z ∞
−∞
dε g(e) ε eβ(ε−µ)−1.
b. For this particular system, we have
g(ε) =V 2π
(hc)2 ε Θ(ε), where Θ(x) = 0, x <0, Θ(x) = 1, x >0. Show that
U =K pV, and determine the purely numerical constant K.
c. Let µ → 0−, and compute the temperature Tλ below which we may have a macroscopic oc- cupation of the ground state (Bose-Einstein condensation).
d. In a 2D non-relativistic system of particles with mass m, with εk = ¯h2k2/2m, one cannot have Bose-Einstein condensation at any temperature Tλ > 0. In problem 1 c, you found a Tλ >0.
Explain on physical grounds why an ultra-relativistic non-interacting 2D boson system can feature Bose-Einstein condensation at finite temperature, while a corresponding non-relativistic system can- not.
Useful formula:
hNi= ∂lnZg
∂(βµ).
Problem 2. Non-interacting classical spins (Points: 10+10+10+10=40)
Consider N classical spins Si = (Six, Siy, Siz) in a constant external magnetic field h = (0,0, h).
Let |Si| = S be the length of the spins, which we take to be equal on all lattice sites, and ignore couplings between spins. Consider a particular case of the model, namely a system of three-state Ising spins, defined by Si = (Six, Siy, Siz) = S(0,0, σi), σi = 0,±1.
The Hamiltonian H and the canonical partition function Z of this spin-system are given by H = −h·
N
X
i=1
Si,
Z = X
{Si}
e−βH.
Here, β = 1/kBT, wherekB is Boltzmann’s constant and T is temperature.
a. Show that the partition functions for this spin-model is given on the form Z = (F(Sβh))N,
where the functional form of F(x) is given by F(x) = 1 + 2 cosh(x).
b. Calculate the enthalpy,He, of the systems, where He =−∂lnZ
∂β .
c. Show that in general, we have for the specific heat at constant magnetic field, Ch ≡(∂He/∂T)h Ch =kBβ2hhH2i − hHi2i,
where hOi ≡(1/Z)P{Si}Oe−βH.
d. Compute Ch in the limit Sβh1. Explain on physical grounds the result you find.
Problem 3. Classical particles in 2D (Points: 5+5+10+10=30)
Consider N non-relativistic classical particles with mass m moving in the two-dimensional (2D) (x, y)-plane with “volume” V = πR2, where R is the radius of the circle to which the N particles are confined. These particles are non-interacting, but are subject to an external potential. The Hamiltonian of the system is given by
H =
N
X
i=1
"
p2i
2m +V0 ri2
#
.
Here, pi = (pxi, pyi) is the momentum of particle i, while ri = qx2i +yi2 is the distance from the center of the volume V of particle i. The canonical partition function Z for this system is given by
Z = 1
N!h2N
Z
..
Z
dΓ e−βH; dΓ =
N
Y
i=1
d2pi d2ri.
Here, h is a constant with dimension Js, while β = 1/kBT, kB is Boltzmann’s constant, and T is temperature. For later use, we also define the length R0 ≡1/√
βV0.
a. Show that the internal energyU =hHiin general is given by, in the canonical ensemble U =−∂lnZ
∂β .
b. Show that in this case, the canonical partition function is given
Z = 1
N! 1 λ2N QN1 ,
λ ≡ h
√2πmkBT, Q1 ≡
Z
d2r e−βV0r2.
c. Calculate the pressure p that the particles exert on the “walls” of the system (the perimeter of the circle) in the limits R/R0 1 andR/R0 1.
d. In polar coordinates, the Hamiltonian may be written H =
N
X
i=1
"
p2ri
2m + p2θi
2mr2i +V0ri2
#
Here, pri is the radial component of the linear momentum of particlei, and pθiis its angular momen- tum. Compute
˜
u≡ h p2θi
2mri2 +V0ri2i; R/R0 1.
Useful formula:
hOi = 1 Z
1 N!h2N
Z
..
Z
dΓ O e−βH.
Useful formulae:
X
k
F(εk) =
Z ∞
−∞ de g(e) F(e) g(e) ≡ X
k
δ(e−εk)
Z
dνr F(|r|) = Ων
Z
dr rν−1 F(r); Ων = 2πν/2 Γ(ν/2) Γ(z) ≡
Z ∞ 0
dx xz−1 e−x Γ(z+ 1) = z Γ(z)
ζ(z) ≡
∞
X
l=1
1 lz
Z a 0
dx xν−1 e−xν = 1 ν
Z aν 0
du e−u
Z ∞ 0
dx xz
ex−1 = ζ(z+ 1) Γ(z+ 1)
Z ∞ 0
dx x e−x = 1 Ch = ∂He
∂T
!
h
=−kBβ2 ∂He
∂β
!
h
Generalized Equipartition Principle:
Let the Hamiltonian of a system be given by H =α|q|ν +H0. Here q is a generalized coordinate or momentum which does not appear in H0. Let the partition function be given by
Z =
Z
dq
Z
dΓ0e−βH, such that we have
hα|q|νi= 1 Z
Z
dq
Z
dΓ0α|q|ν e−βH. Then we have
hα|q|νi= kBT ν .
Three-dimensional volume element in spherical coordinates:
d3r=dΩ r2dr; dΩ = dθ sinθ dφ Here, θ is a polar angle and φ is an azimuthal angle.