FACTOID # 27: Want your kids to stay in school? Send them to Norway.
 
 Home   Encyclopedia   Statistics   Countries A-Z   Flags   Maps   Education   Forum   FAQ   About 
 
WHAT'S NEW
RECENT ARTICLES
More Recent Articles »
 

SEARCH ALL

FACTS & STATISTICS    Advanced view

Search encyclopedia, statistics and forums:

 

 

(* = Graphable)

 

 


Encyclopedia > Debye model

In thermodynamics and solid state physics, the Debye model is a method developed by Peter Debye in 1912[1] for estimating the phonon contribution to the specific heat (heat capacity) in a solid. The Debye model treats the vibrations of the atomic lattice (heat) as phonons in a box, in contrast to Einstein model, which treats the solid as many individual, non-interacting quantum harmonic oscillators. This model correctly predicts the low temperature dependence of the heat capacity, which is proportional to T3. Just like the Einstein model, it also recovers the Dulong-Petit law at high temperatures. But due to simplifying assumptions, its accuracy suffers at intermediate temperatures. Thermodynamics (from the Greek thermos meaning heat and dynamis meaning power) is a branch of physics that studies the effects of changes in temperature, pressure, and volume on physical systems at the macroscopic scale by analyzing the collective motion of their particles using statistics. ... Solid-state physics, the largest branch of condensed matter physics, is the study of rigid matter, or solids. ... Peter Joseph William Debye (March 24, 1884 - November 2, 1966) (born Petrus Josephus Wilhelmus Debije) was a Dutch physical chemist. ... In physics, a phonon is a quantized mode of vibration occurring in a rigid crystal lattice, such as the atomic lattice of a solid. ... The specific heat capacity (symbol c or s, also called specific heat) of a substance is defined as heat capacity per unit mass. ... In jewelry, a solid gold piece is the alternative to gold-filled or gold-plated jewelry. ... See Oscillator (disambiguation) for particular types of oscillation and oscillators. ... Rose des Sables (Sand Rose), formed of gypsum crystals In mineralogy and crystallography, a crystal structure is a unique arrangement of atoms in a crystal. ... In physics, a phonon is a quantized mode of vibration occurring in a rigid crystal lattice, such as the atomic lattice of a solid. ... This article needs to be cleaned up to conform to a higher standard of quality. ... The quantum harmonic oscillator is the quantum mechanical analogue of the classical harmonic oscillator. ... This article needs to be cleaned up to conform to a higher standard of quality. ... The Dulong-Petit law, found in 1819 by Pierre Louis Dulong and Alexis Thérèse Petit, states the classical expression for the specific heat capacity of a crystal due to its lattice vibrations. ...

Contents


Derivation

The Debye model is a solid-state equivalent of Planck's law of black body radiation, where one treats electromagnetic radiation as a gas of photons in a box. The Debye model treats atomic vibrations as phonons in a box (the box being the solid). Most of the calculation steps are identical. Black body spectrum as a function of wavelength In physics, the spectral intensity of electromagnetic radiation from a black body at temperature T is given by the Plancks law of black body radiation: where: I(ν) is the amount of energy per unit time per unit surface area per unit... Electromagnetic radiation can be conceptualized as a self propagating transverse oscillating wave of electric and magnetic fields. ... The results of the quantum particle in a box can be used to look at the equilibrium situation for a quantum ideal gas in a box which is a box containing a large number particles which do not interact with each other except for instantaneous thermalizing collisions. ... In physics, a phonon is a quantized mode of vibration occurring in a rigid crystal lattice, such as the atomic lattice of a solid. ...


Consider a cube of side L. From the particle in a box article, the resonating modes of the sonic disturbances inside the box (considering for now only those aligned with one axis) have wavelengths given by In physics the particle in a box refers to a simple mathematical exercise performed in quantum mechanics. ...

lambda_n = {2Lover n}

where n is an integer. The energy of a phonon is

E_n =hnu_n

where h is Planck's constant and νn is the frequency of the phonon. We make the approximation that the frequency is inversely proportional to the wavelength, giving: A commemoration plaque for Max Planck on his discovery of Plancks constant, in front of Humboldt University, Berlin. ...

E_n=hnu_n={hc_soverlambda_n}={hc_snover 2L}

in which cs is the speed of sound inside the solid. In three dimensions we will use:

E_n^2=E_{nx}^2+E_{ny}^2+E_{nz}^2=left({hc_sover2L}right)^2left(n_x^2+n_y^2+n_z^2right)

The approximation that the frequency is inversely proportional to the wavelength (giving a constant speed of sound) is good for low-energy phonons but not for high-energy phonons. (See the article on phonons.) This is one of the limitations of the Debye model. In physics, a phonon is a quantized mode of vibration occurring in a rigid crystal lattice, such as the atomic lattice of a solid. ...


Let's now compute the total energy in the box

U = sum_n E_n,bar{N}(E_n)

where bar{N}(E_n) is the number of phonons in the box with energy En. In other words, the total energy is equal to the sum of energy multiplied by the number of phonons with that energy (in one dimension). In 3 dimensions we have:

U = sum_{n_x}sum_{n_y}sum_{n_z}E_n,bar{N}(E_n)

Now, this is where Debye model and Planck's law of black body radiation differ. Unlike electromagnetic radiation in a box, there is a finite number of phonon energy states because a phonon cannot have infinite frequency. Its frequency is bound by the medium of its propagation -- the atomic lattice of the solid. Consider an illustration of a transverse phonon below. Black body spectrum as a function of wavelength In physics, the spectral intensity of electromagnetic radiation from a black body at temperature T is given by the Plancks law of black body radiation: where: I(ν) is the amount of energy per unit time per unit surface area per unit... In physics, a phonon is a quantized mode of vibration occurring in a rigid crystal lattice, such as the atomic lattice of a solid. ... In physics, a phonon is a quantized mode of vibration occurring in a rigid crystal lattice, such as the atomic lattice of a solid. ...

It is reasonable to assume that the minimum wavelength of a phonon is twice the atom separation, as shown in the lower figure. There are N atoms in a solid. Our solid is a cube, which means there are sqrt[3]{N} atoms per side. Atom separation is then given by L/sqrt[3]{N}, and the minimum wavelength is Download high resolution version (800x800, 60 KB)Lattice vibrations. ... In physics, a phonon is a quantized mode of vibration occurring in a rigid crystal lattice, such as the atomic lattice of a solid. ...

lambda_{rm min} = {2L over sqrt[3]{N}}

making the maximum mode number n (infinite for photons) In physics, the photon (from Greek φως, phōs, meaning light) is the quantum of the electromagnetic field; for instance, light. ...

n_{rm max} = sqrt[3]{N}

This is the upper limit of the triple energy sum

U = sum_{n_x}^{sqrt[3]{N}}sum_{n_y}^{sqrt[3]{N}}sum_{n_z}^{sqrt[3]{N}}E_n,bar{N}(E_n)

For slowly-varying, well-behaved functions, a sum can be replaced with an integral (a.k.a Thomas-Fermi approximation)

U approxint_0^{sqrt[3]{N}}int_0^{sqrt[3]{N}}int_0^{sqrt[3]{N}} E(n),bar{N}left(E(n)right),dn_x, dn_y, dn_z

So far, there has been no mention of bar{N}(E), the number of phonons with energy E. Phonons obey Bose-Einstein statistics. Their distribution is given by the famous Bose-Einstein formula In statistical mechanics, Bose-Einstein statistics determines the statistical distribution of identical indistinguishable bosons over the energy states in thermal equilibrium. ...

langle Nrangle_{BE} = {1over e^{E/kT}-1}

Because a phonon has three possible polarization states (one longitudinal and two transverse) which do not affect its energy, the formula above must be multiplied by 3 Longitudinal waves, also referred to as compressional waves or pressure waves, are waves that have vibrations along or parallel to their direction of travel. ... A light wave is an example of a transverse wave. ...

bar{N}(E) = {3over e^{E/kT}-1}

Substituting this into the energy integral yields

U = int_0^{sqrt[3]{N}}int_0^{sqrt[3]{N}}int_0^{sqrt[3]{N}} E(n),{3over e^{E(n)/kT}-1},dn_x, dn_y, dn_z

The ease with which these integrals are evaluated for photons is due to the fact that light's frequency, at least semi-classically, is unbound. As the figure above illustrates, this is not true for phonons. In order to approximate this triple integral, Debye used spherical coordinates In physics, the photon (from Greek φως, phōs, meaning light) is the quantum of the electromagnetic field; for instance, light. ... In physics, a phonon is a quantized mode of vibration occurring in a rigid crystal lattice, such as the atomic lattice of a solid. ... Peter Joseph William Debye (March 24, 1884 - November 2, 1966) (born Petrus Josephus Wilhelmus Debije) was a Dutch physical chemist. ...

 (n_x,n_y,n_z)=(ncos theta cos phi,ncos theta sin phi,nsin theta )

and boldly approximated the cube by an eighth of a sphere

U approxint_0^{pi/2}int_0^{pi/2}int_0^R E(n),{3over e^{E(n)/kT}-1}n^2 sintheta, dn, dtheta, dphi

where R is the radius of this sphere, which is found by conserving the number of particles in the cube and in the eighth of a sphere. The volume of the cube is N unit-cell volumes,

N = {1over8}{4over3}pi R^3

so we get:

R = sqrt[3]{6Noverpi}

The substitution of integration over a sphere for the correct integral introduces another source of inaccuracy into the model.


The energy integral becomes

U = {3piover2}int_0^R ,{hc_snover 2L}{n^2over e^{hc_sn/2LkT}-1} ,dn

Changing the integration variable to x = {hc_snover 2LkT},

U = {3piover2} kT left({2LkTover hc_s}right)^3int_0^{hc_sR/2LkT} {x^3over e^x-1}, dx

To simplify the look of this expression, define the Debye temperature TD -- a shorthand for some constants and material-dependent variables.

T_Dequiv {hc_sRover2Lk} = {hc_sover2Lk}sqrt[3]{6Noverpi} = {hc_sover2k}sqrt[3]{{6overpi}{Nover V}}

We then have the specific internal energy:

frac{U}{Nk} = 9T left({Tover T_D}right)^3int_0^{T_D/T} {x^3over e^x-1}, dx = 3T D_3 left({T_Dover T}right)

where D3(x) is the (third) Debye function. In mathematics, the family of Debye functions is defined by Categories: Math stubs ...


Differentiating with respect to T we get the dimensionless heat capacity:

frac{C_V}{Nk} = 9 left({Tover T_D}right)^3int_0^{T_D/T} {x^4 e^xoverleft(e^x-1right)^2}, dx

These formulae give the Debye model at all temperatures. The more elementary formulae given further down give the asymptotic behavior in the limit of low and high temperatures.


Debye's derivation

Actually, Debye derived his equation somewhat differently and more simply. Using the solid mechanics of a continuous medium, he found that the number of vibrational states with a frequency less than a particular value was asymptotic to Solid mechanics is the branch of physics and mathematics that concern the behavior of solid matter under external actions (e. ...

n sim {1 over 3} nu^3 V F

in which V is the volume and F is a factor which he calculated from elasticity coefficients and density. Combining this with the expected energy of a harmonic oscillator at temperature T (already used by Einstein in his model) would give an energy of Elasticity has meanings in two different fields: In physics and mechanical engineering, the theory of elasticity describes how a solid object moves and deforms in response to external stress. ... This article needs to be cleaned up to conform to a higher standard of quality. ...

U = int_0^infty ,{hnu^3 V Fover e^{hnu/kT}-1}, dnu

if the vibrational frequencies continued to infinity. This form gives the T4 behavior which is correct at low temperatures. But Debye realized that there could not be more than 3N vibrational states for N atoms. He made the assumption that in an atomic solid, the spectrum of frequencies of the vibrational states would continue to follow the above rule, up to a maximum frequency νm chosen so that the total number of states is 3N:

3N = {1 over 3} nu_m^3 V F

Debye knew that this assumption was not really correct (the higher frequencies are more closely spaced than assumed), but it guarantees the proper behavior at high temperature (the Dulong-Petit law). The energy is then given by: The Dulong-Petit law, found in 1819 by Pierre Louis Dulong and Alexis Thérèse Petit, states the classical expression for the specific heat capacity of a crystal due to its lattice vibrations. ...

U = int_0^{nu_m} ,{hnu^3 V Fover e^{hnu/kT}-1}, dnu
= V F kT (kT/h)^3 int_0^{T_D/T} ,{x^3 over e^x-1}, dx
where TD is hνm / k.
= 9 N k T (T/T_D)^3 int_0^{T_D/T} ,{x^3 over e^x-1}, dx
= 3NkTD3(TD / T)

where D3 is the function later given the name of third-order Debye function. In mathematics, the family of Debye functions is defined by Categories: Math stubs ...


Low temperature limit

The temperature of a Debye solid is said to be low if T ll T_D, leading to

frac{C_V}{Nk} sim 9 left({Tover T_D}right)^3int_0^{infty} {x^4 e^xover left(e^x-1right)^2}, dx

This definite integral can be evaluated exactly:

frac{C_V}{Nk} sim {12pi^4over5} left({Tover T_D}right)^3

In the low temperature limit, the limitations of the Debye model mentioned above do not apply, and it gives a correct relationship between (phononic) heat capacity, temperature, the elastic coefficients, and the volume per atom (the latter quantities being contained in the Debye temperature).


High temperature limit

The temperature of a Debye solid is said to be high if T > > TD. e^x - 1approx x if | x | < < 1, leads to

frac{C_V}{Nk} sim 9 left({Tover T_D}right)^3int_0^{T_D/T} {x^4 over x^2}, dx
frac{C_V}{Nk} sim 3

This is the Dulong-Petit law, and is fairly accurate although it does not take into account anharmonicity, which causes the heat capacity to rise further. The total heat capacity of the solid, if it is a conductor or semiconductor, may also contain a non-negligible contribution from the electrons. The Dulong-Petit law, found in 1819 by Pierre Louis Dulong and Alexis Thérèse Petit, states the classical expression for the specific heat capacity of a crystal due to its lattice vibrations. ... Electrical conduction is the current (movement of charged particles) through a material in response to an electric field. ... A semiconductor is a material with an electrical conductivity that is intermediate between that of an insulator and a conductor. ...


Debye versus Einstein

Debye vs. Einstein. Predicted heat capacity as a function of temperature.
Debye vs. Einstein. Predicted heat capacity as a function of temperature.

So how closely do the Debye and Einstein models correspond to experiment? -- Surprisingly close, but Debye is correct at low temperatures whereas Einstein is not. Debye vs. ... Debye vs. ...


How different are the models? To answer that question one would naturally plot the two on the same set of axes... except one can't. Both the Einstein model and the Debye model provide a functional form for the heat capacity. They are models, and no model is without a scale. A scale relates the model to its real-world counterpart. One can see that the scale of the Einstein model, which is given by

C_V = 3Nkleft({epsilonover k T}right)^2{e^{epsilon/kT}over left(e^{epsilon/kT}-1right)^2}

is ε / k. And the scale of the Debye model is TD, the Debye temperature. Both are usually found by fitting the models to the experimental data. (The Debye temperature can theoretically be calculated from the speed of sound and crystal dimensions.) Because the two methods approach the problem from different directions and different geometries, Einstein and Debye scales are not the same, that is to say

{epsilonover k} ne T_D

which means that plotting them on the same set of axes makes no sense. They are two models of the same thing, but of different scales. If one defines Einstein temperature as

T_E equiv {epsilonover k}

then one can say

T_E ne T_D

and, to relate the two, we must seek the ratio

{T_Eover T_D} = ?

The Einstein solid is composed of single-frequency quantum harmonic oscillators, epsilon = hbaromega = hnu. That frequency, if it indeed existed, would be related to the speed of sound in the solid... even though there is no sound in Einstein solid. If one imagines the propagation of sound as a sequence of atoms hitting one another, then it becomes obvious that the frequency of oscillation must correspond to the minimum wavelength sustainable by the atomic lattice, λmin. This article needs to be cleaned up to conform to a higher standard of quality. ... The quantum harmonic oscillator is the quantum mechanical analogue of the classical harmonic oscillator. ...

nu = {c_soverlambda} = {c_ssqrt[3]{N}over 2L} = {c_sover 2}sqrt[3]{Nover V}

which makes the Einstein temperature

T_E = {epsilonover k} = {hnuover k} = {h c_sover 2k}sqrt[3]{Nover V}

and the sought ratio is therefore

{T_Eover T_D} = sqrt[3]{piover6}

Now both models can be plotted on the same graph. Note that this ratio is the cube root of the ratio of the volume of one octant of a 3-dimensional sphere to the volume of the cube that contains it, which is just the correction factor used by Debye when approximating the energy integral above.


Debye temperature table

Even though the Debye model is not completely correct, it gives a good approximation for the low temperature heat capacity of insulating, crystalline solids where other contributions (such as highly mobile conduction electrons) are negligible. For metals, the electron contribution to the heat capacity goes as T, which at low temperatures, dominates the debye T^3 result for lattice vibrations. In this case, the debye model can only be said to approximate the lattice contribution to the specific heat. The following table lists Debye temperatures for several substances:

Aluminum 426K
Cadmium 186K
Chromium 610K
Copper 344.5K
Gold 165K
α-Iron 464K
Lead 96K
α-Manganese 476K
Nickel 440K
Platinum 240K
Silicon 640K
Silver 225K
Tin (white) 195K
Titanium 420K
Tungsten 405K
Zinc 300K
Diamond 2200K
Ice 192K

See also

An ideal Bose gas is a quantum-mechanical version of a classical ideal gas. ... The results of the quantum particle in a box can be used to look at the equilibrium situation for a quantum ideal gas in a box which is a box containing a large number particles which do not interact with each other except for instantaneous thermalizing collisions. ...

References

  •  'Zur Theorie der spezifischen Warmen', Annalen der Physik 39(4), p. 789 (1912)
  • CRC Handbook of Chemistry and Physics, 56th Edition (1975-1976)

External links

  • Experimental determination of specific heat, thermal and heat conductivity of quartz using a cryostat.

  Results from FactBites:
 
Debye Theory of Specific Heat (428 words)
The density of states for these modes, which are called "phonons", is of the same form as the photon density of states in a cavity.
The full expression for the Debye specific heat must be evaluated by numerical procedures.
Explaining the drastic departure from the Law of Dulong and Petit was a major contribution of the Einstein and Debye models.
Debye model - Wikipedia, the free encyclopedia (1413 words)
The Debye model is a solid-state equivalent of Planck's law of fl body radiation, where one treats electromagnetic radiation as a gas of photons in a box.
Debye knew that this assumption was not really correct (the higher frequencies are more closely spaced than assumed), but it guarantees the proper behavior at high temperature (the Dulong-Petit law).
Even though the Debye model is not completely correct, it gives a good approximation for the low temperature heat capacity of insulating, crystalline solids where other contributions (such as highly mobile conduction electrons) are negligible.
  More results at FactBites »


 

COMMENTARY     


Share your thoughts, questions and commentary here
Your name
Your comments
Please enter the 5-letter protection code

Want to know more?
Search encyclopedia, statistics and forums:

 


Lesson Plans | Student Area | Student FAQ | Reviews | Press Releases |  Feeds | Contact
The Wikipedia article included on this page is licensed under the GFDL.
Images may be subject to relevant owners' copyright.
All other elements are (c) copyright NationMaster.com 2003-5. All Rights Reserved.
Usage implies agreement with terms.