1. Not finding help here? Sign up for a free 30min tutor trial with Chegg Tutors
    Dismiss Notice
Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Heat capacity vs temperature

  1. Dec 8, 2015 #1

    MMS

    User Avatar

    Hello,
    I'm reposting this 'cause I think it's more relevant here than where it was posted before. Not sure how to delete the other one...

    I'm trying to plot the heat capacity as a function of the temperature from some small temperature to Debye's using numerical summation over the energy (shown in the photo below).
    However, I'm struggling to determine the boundary of k (actually, of all three since there are 3 sums).


    M7dVt2S.png

    I'm using Maple for those wondering.

    I'd be more than happy if someone could help me out with that.

    Thank you in advance.
     
  2. jcsd
  3. Dec 9, 2015 #2
    In a periodic crystal the largest k of lattice vibrations is pi/a. (I am not sure I understand the question.)
     
  4. Dec 9, 2015 #3

    MMS

    User Avatar

    Yes, so I believe. I summed over the FBZ (-pi/a to pi/a) for each of k's products.

    I'm yet to be able to work this out on Maple though.

    Do you have any experience with it? Or if anyone reading this knows some Maple, how will I be able to insert the expression for C(T) above (I can upload what I've tried if needed) and actually have it work (I'm not getting anything)? :P
     
  5. Dec 10, 2015 #4
    It is sufficient to sum over one irreducible wedge of the Brillouin zone, that saves some computing time.
     
  6. Dec 10, 2015 #5

    MMS

    User Avatar

    Hi,

    I tried plotting it first in 1D (hence sin(k*a/2) where k has only one product) and i receive the expected graph of C(T) vs T.

    However, when I plot it in 2D and 3D it becomes this deadly noise-like signal (example shown below for what I got in 3D).
    hoHEj05.png
    Any idea where I could be mistaken? I literally just plugged the expression in Maple..
     
  7. Dec 10, 2015 #6
    So that is basically a constant large value. With numerical noise that increases linearly with temperature.
    And I assume that you know that it should start out as T^3, approaching the Dulong-Petit value asymptotically for hight T.

    Sorry, I do not know anything about Maple.
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook

Have something to add?
Draft saved Draft deleted