Loading [MathJax]/jax/output/HTML-CSS/jax.js
Skip to main content
Library homepage
 

Text Color

Text Size

 

Margin Size

 

Font Type

Enable Dyslexic Font
Engineering LibreTexts

Debye Model For Specific Heat

( \newcommand{\kernel}{\mathrm{null}\,}\)

The Debye model is a method developed by Peter Debye in 1912[7] for estimating the phonon contribution to the specific heat (heat capacity) in a solid[1]. This model correctly explains the low temperature dependence of the heat capacity, which is proportional to T3 and also recovers the Dulong-Petit law at high temperatures. However, due to simplifying assumptions, its accuracy suffers at intermediate temperatures.

Introduction to Phonons

In 1912 Debye realized that, inconsistent with the Einstein model, low-energetic excitations of a solid material were not oscillations of a single atom[2], but collective modes propagating through the material. Such vibrations can be considered to be sound waves, and their propagation speed is the speed of sound in the material[3]. Moreover, these modes only accept energy in discreet amounts.

phonon.gif
Figure 1: Collective motions called phonons (Public Domain; Sean Kelley via NIST)

Quantum theory uses the concepts of phonons, which are “quasi-particles” with definite energies and directions of motion, to treat the vibrations. The concept of phonon is analogous with photons of the electromagnetic wave. The relations between the energy of a phonon ε, the angular frequency ω and the wave vector q are:

ε=ω

ω=υs|q|

where vs is the velocity of the sound wave.

Phonons obey Bose–Einstein statistics. The expectation number of bosons in a state with energy E is[6]:

nE=1eE/kBT1=1eω/kBT1

where kB= 1.380 6504(24)×10−23J/K is the Boltzmann constant.

Debye frequency and Debye Temperature

Unlike electromagnetic radiation in a box, a phonon cannot have infinite frequency. Its frequency is bound by the medium of its propagation — the atomic lattice of the solid. If there are N primitive cells in the specimen, the total number of phonon modes are N. A cut-off frequency ωD, known as Debye frequency, is determined by the following manner[4]:

In the 3 dimensional reciprocal space, the volume for each allowed wave vector q is:

(2πL)3=8π3V

Since there is a cut-off wave vector qD=ωD/υs, all the modes are confined within a sphere with radius qD. Thus number of modes (not number of phonons) should be (5)

N=(43πqD3)/(8π3V)

or

qD=(6π2NV)13

ωD=υs(6π2NV)13

Debye temperature TD is defined as

TD=ωDkB=υskB(6π2NV)13

The significance of this physical term will be discussed below. For elements in the same group, heavier atoms have lower Debye temperatures, simply because the velocity of sound decreases as the density increases.[4] The Debye temperatures of several substances are listed in Table 1.[1]

Table 1: Debye Temperatures of several substances.
Aluminum 428K Iron 470K Silicon 645K Tungsten 400K
Cadmium 209K Lead 105K Silver 225K Zinc 327K
Chromium 630K Manganese 410K Tantalum 240K Carbon 2230K
Copper 343.5K Nickel 450K Tin(white) 200K Ice 192K
Gold 165K Platinum 240K Titanium 420K
Note

Heavier atoms have lower Debye temperatures, because the velocity of sound decreases as the density increases.

Derivation for Specific Heat

In the Debye approximation, the velocity of sound υs is taken as constant for each polarization type, as it would be for a classical elastic continuum. According to equation (7), the density of states is:

D(ω)=dNdω=Vω22π2υ3s

Thus, thermal energy for each polarization type is given by[4]:

U=dωD(ω)n(ω)ω=ωD0dωVω22π2υ3sωeω/kBT1

There are 3 polarization types, 1 longitude and 2 transverse. For brevity, we assume that phonon velocity is independent of polarization. Thus we multiply a factor 3 and use equation (7) to obtain the total energy of phonon:

U=3V2π2υ3sωD0dωω3eω/kBT1=3Vk4BT42π2υ3s3xD0dxx3ex1=9NkBT(TTD)3xD0dxx3ex1

where x=ω/kBT and xDTD/T.

The heat capacity is:

CV=UT=9NkB(TTD)3xD0dxx4ex(ex1)2

At the left of Figure 1[3] below, the experimental results of specific heats of four substances are plotted as a function of temperature and they look very different. But if they are scaled to T/TD, they look very similar and are very close to the Debye theory.

Screen Shot 2019-02-03 at 12.52.22 AM.png

Figure 1: Specific Heats of Lead, Silver, Aluminum and Diamond

High and Low Temperature Limits

The integral in Equation (12) cannot be evaluated in closed form. But the high and low temperature limits can be assessed.

High Temperature Limit

For the high temperature case where TTD, the value of x is very small throughout the range of the integral. This justifies using the approximation to the exponential ex1+x and reduces equation (11) and (12) to

U=9NkBT(TTD)3xD0x2dx=3NkBT

CV=3NkB

which is the classical Dulong-Petit result.

When the temperature is above the Debye temperature, the heat capacity is very close to the classical value 3NkBT. For temperatures below the Debye temperature, quantum effects become important and Cv decreases to zero. Note that diamond, with a Debye temperature of 1860K, is a “quantum solid” at room temperature. [8]

Low Temperature Limit

At very low temperature where TTD, only long wavelength acoustic modes are thermally excited. These are just the modes that can be treated as elastic continuum with macroscopic elastic constants. The energy of those short wavelength modes are too high to be populated significantly at low temperatures. We may approximate xDTD/T to infinity and make use of the standard integral

0dxx3ex1=π415

to obtain

U=3π4NkBT45T3D

CV=12π4NkBT35T3D324NkBT3T3D

For actual crystals, the temperatures at which the T3 approximation holds are quite low, even may be below TD/50. It is easy to understand T/TD3 with a simple argument. Only the modes with ω<kBT will be excited to any appreciable extent at a low temperature T. Define ωTkT/ as the largest frequency excited at this temperature. In the reciprocal space, the fraction occupied by the excited states is (qT/qD)3 or (wT/wD)3=(T/TD)3.

Extension: Einstein-Debye Specific Heat

This T dependence of the specific heat at very low temperatures agrees with experiment for nonmetals. For metals the specific heat of highly mobile conduction electrons is approximated by Einstein Model, which is composed of single-frequency quantum harmonic oscillators. The temperature dependence of Einstein model is just T. It becomes significant at low temperatures and is combined with the above lattice specific heat in the Einstein-Debye specific heat[3].

Cmetal=Celectron+Cphonon=π2Nk22EfT+12π4NkB5T3DT3

Finally, experiments suggest that amorphous materials do not follow the Debye T3 law even at temperatures below 0.01TD[8]. There is more yet to be learned.

References

  1. http://en.Wikipedia.org/wiki/Debye_model
  2. http://www.uam.es/personal_pdi/cienc...ivos/debye.pdf
  3. http://hyperphysics.phy-astr.gsu.edu...ds/phonon.html
  4. C. Kittel, introduction to solid state physics, pp.117-140. John Wiley & Sons, Inc., 1996
  5. L. Landau and E. Lifshitz, Statistical Physics, Part, pp.191 -217. Elsevier, 1984
  6. http://en.Wikipedia.org/wiki/Bose%E2...ein_statistics
  7. Zur Theorie der spezifischen Waerm, Annalen der Physik (Leipzig) 39(4)
  8. www.phys.unsw.edu.au/~gary/SM3_6.pdf

Contributors and Attributions

  • ContribMSE5317


Debye Model For Specific Heat is shared under a CC BY-SA license and was authored, remixed, and/or curated by LibreTexts.

  • Was this article helpful?

Support Center

How can we help?