Homework Set 2

From phys813
Jump to navigationJump to search

Problem 1: Canonical partition function for a single non-interacting spin

The Hamiltonian for an electron spin degree of freedom in the external magnetic field 𝐁 is given by:

H^=gμBσ^𝐁

where μB is the Bohr magneton, g is the gyromagnetic ratio, and σ^=(σ^x,σ^y,σ^z) is the vector of the Pauli matrices:

σ^x=(0110),

σ^y=(0ii0),

σ^z=(1001).

(a) In the quantum canonical ensemble, evaluate the density matrix if 𝐁 is along the z axis.

(b) Repeat the calculation from (a) assuming that 𝐁 points along the x axis.

(c) Calculate the average energy in each of the above cases.

Problem 2: Canonical partition function for two-interacting spins

In some antiferromagnetic materials, such as TiCuCl3, spins S=1/2 are arranged in pairs. To first approximation, such dimers can be considered independently of each other. The Hamiltonian of a single dimer in the external magnetic field 𝐁=(0,0,B) is given by:

H^=J𝐒^1𝐒^2+2μBB(S^1z+S^2z)

where J>0 is the exchange coupling constant and μB is the Bohr magneton.

(a) How many energy eigenlevels does this Hamiltonian have? List all eigenergies explicitly.

(b) Using your result in (a), compute the canonical partition function and free energy of N dimers, as well as find their entropy.

NOTE: A more pedagogical expression for 𝐒^1𝐒^2 is S^1xS^2x+S^1yS^2y+S^1zS^2z.

Problem 3: Density matrix and canonical partition function for one-dimensional harmonic oscillator

Consider a single one-dimensional quantum harmonic oscillator described by the Hamiltonian:

H^=p^22m+mω2q22,

where p^=iddq.

(a) Find the partition function Z in the quantum canonical ensemble at temperature T.

(b) Using the result from (a), calculate the averge energy E=H^.

(c) Write down the formal expression for the canonical density operator ρ^ in terms of the eigenstates |n of the Hamiltonian and the corresponding energy levels εn=ω(n+1/2).

(d) Using the result in (c), write down the density matrix in the coordinate representation q|ρ^|q.

(e) In the coordinate representation, calculate explicitly q|ρ^|q in the high temperature limit T.

HINT: One approach is to utilize the following result

eβA^eβB^=eβ(A^+B^)+β2[A^,B^]/2+O(β3)

which you can apply to the Boltzmann operator:

eβH^=eβp^22mβmω2q22

while neglecting terms of order β2 and higher since β is very small in the high temperature limit.

(f) At low temperatures, ρ^ is dominated by low-energy states. Use the ground state wave function q|0 only, evaluate the limiting behavior of q|ρ^|q as T0.