Read the material form the textbook-energy distributions and partition function

B=0.037;ωe=212;

B=10. ;ωe=2000;

(1) Show Fig representig rotational and the other vibrational energy distribution:

(2)

RotEnergy[a_,β_]:= ((-(∑_(J=0)^(J=700))(2J+1)*(∂_β) Exp[-a*J*(J+1)*β])/((∑_(J=0)^(J=700))(2J+1)*Exp[-a*J*(J+1)*β]))

RotHeatCapacity[β_]=-β^2*(∂_β) RotEnergy[B,β];

VibEnergy[a_,β_]:= ((-(∑_(v=0)^(v=100))(∂_β) Exp[-a*v*β])/((∑_(v=0)^(v=100))Exp[-a*v*β]));

VibHeatCapacity[β_]=-β^2*(∂_β) VibEnergy[ωe,β];

(3)

Plot[RotHeatCapacity[1/T],{T,0.01,100},PlotRange->Full,PlotLabel->"Rotational Heat Capacity"]

5. 202/2024

5N[ωe/k]

Plot[VibHeatCapacity[1/T],{T,0.01,1200},PlotRange->Full,PlotLabel->"Vibrational Heat Capacity"]

Plot[VibHeatCapacity[1/T]+RotHeatCapacity[1/T],{T,0.01,1200},PlotRange->Full,PlotLabel->"Total Heat Capacity"]