A model to describe the vibrations of atoms in a solid is to assume that the atoms are isotropic harmonic oscillators and that the vibrations are independent of the vibrations of the other atoms. We use this model to describe the entropy and heat capacity of Bohrium (B). The entropy is given as 9.797 J/(mol K) at T = 400K
What is the heat capacity of B at T = 600K?
Energy level of a harmonic oscillator: En = hw(1/2 + n) (1 dimension)
Z = the partition function = [tex]\sum[/tex]e^(-En/k*T) = 1/(2*sinh(hw/k*T))
F = -k*T*ln(Z)
The partition function above is valid for a single atom which is only moving in one dimention but to get the correct parition function for motion in 3D one simply does Zcorrect = Z³ which gives the free energy as
Fcorrect = -3*k*T*ln(Z)
S = -([tex]\partial[/tex]F/[tex]\partial[/tex]T)
Cv = T([tex]\partial[/tex]S/[tex]\partial[/tex]T)
The Attempt at a Solution
From the partition function i can derive the function for S in a straight forward way. I am given the entropy at T = 400 and what i need to know to derive the heat capacity Cv from S is the value for hw. But i cant figure out how to get hw from S given the entropy. I have tried solving the function by the Newton Rhapson method using MatLab but it doesnt want to work. So im stuck and im wondering if im even on the right track. Is there a simple way to find the energy levels of a vibration as in this case when given the entropy?