Modelo de Debye

Modelo de Debye

En la termodinámica y física del estado sólido, el modelo de Debye es un método desarrollado por Peter Debye en 1912[1] para la estimación de la contribución de los fonones al calor específico en un sólido. Trata las vibraciones de la red atómica (calor) como fonones en una caja, en contraste con el modelo de Einstein, que modeliza los sólidos como muchos osciladores armónicos cuánticos no interactuantes entre sí. El modelo de Debye predice correctamente la dependencia a temperaturas bajas de la capacidad calorífica, que es proporcional a T3. Al igual que el modelo de Einstein, también recupera la ley de Dulong-Petit a altas temperaturas; sin embargo, debido a la simplicidad de los supuestos sobre los que se apoya, resulta deficiente a la hora de explicar los fenómenos observables a temperaturas intermedias.

Contenido

Desarrollo teórico

El de Debye es un modelo de la física del estado sólido equivalente a la ley de Planck de la radiación del cuerpo negro, que trata la radiación electromagnética como un gas de fotones en una caja. El modelo de Debye trata las vibraciones atómicas como fonones en una caja (la caja es el sólido). La mayor parte del desarrollo teórico es idéntica.

Consideremos un cubo de lado L. Del artículo partícula en una caja sabemos que la resonancia de los modos de las perturbaciones sonoras en el interior de la caja (considerando por ahora sólo los alineados con alguno de los ejes) tienen longitudes de onda dadas por

\lambda_n = {2L\over n}\,,


donde n es un entero. La energía de un fonón es

E_n\ =h\nu_n\,,

donde h es la constante de Planck y νn es la frecuencia del fonón. Hacemos la aproximación de que la frecuencia es inversamente proporcional a la longitud de onda, de lo que resulta:

E_n=h\nu_n={hc_s\over\lambda_n}={hc_sn\over 2L}\,,

donde cs es la velocidad del sonido en el sólido. En tres dimensiones emplearemos:

E_n^2=E_{nx}^2+E_{ny}^2+E_{nz}^2=\left({hc_s\over2L}\right)^2\left(n_x^2+n_y^2+n_z^2\right)\,.

La aproximación de la frecuencia como inversamente proporcional a la longitud de onda (para una velocidad del sonido fija) funciona para fonones de baja energía, pero no para fonones de alta energía (véase el artículo sobre fonones ). Esta es una de las limitaciones del modelo, y corresponde a un fallo de las predicciones para temperaturas intermedias, mientras que tanto para bajas como para altas temperaturas son exactas.

Vamos ahora a calcular la energía total en la caja,

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

donde \bar{N}(E_n) es el número de fonones en la caja con energía En. En otras palabras, el total de la energía es igual a la suma de la energía multiplicada por el número de fonones con esa energía (en una dimensión). En 3 dimensiones tenemos:

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

Aquí es donde el modelo de Debye y la ley de Planck de la radiación del cuerpo negro difieren. Al contrario de lo que pasa con la radiación electromagnética en una caja, existe un número finito de estados de energía posibles para los fonones, pues un fonón no puede tener frecuencia infinita. Su frecuencia está limitada por el medio por el que se propaga (la red atómica del sólido). Observemos la siguiente ilustración de un fonón propagándose transversalmente.

Debye limit.svg

Es razonable asumir que la longitud de onda mínima de un fonón sea del doble de la separación entre átomos, como se observa en la figura inferior. Hay N átomos en un sólido. Nuestro sólido tiene forma cúbica, lo que significa que existen \sqrt[3]{N} átomos por lado. La separación entre átomos viene dada entonces por L/\sqrt[3]{N}, y la longitud de onda mínima será

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

tomando el modo de mayor orden n (que sería infinito en el caso de los fotones)

n_{\rm max} = \sqrt[3]{N}\,.

Este es el límite superior para la triple sumatoria de la energía

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)\,.

Para funciones suaves y de variación lenta, la suma puede ser reemplazada por una integral (esto se conoce como aproximación de Thomas-Fermi)

U \approx\int_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\,.

Hasta el momento, no ha habido mención alguna a \bar{N}(E), el número de fonones con energía E\,.. Los fonones obedecen la estadística de Bose-Einstein. Su distribución está dada por la famosa fórmula de Bose-Einstein

\langle N\rangle_{BE} = {1\over e^{E/kT}-1}\,.

Dado que un fonón tiene tres estados posibles de polarización (uno longitudinal y dos transversales que prácticamente no afectan a su energía) hemos de multiplicar la fórmula anterior por 3,

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

En realidad se utiliza una velocidad sónica efectiva cs: = ceff, es decir, que la temperatura de Debye Td (ver más adelante) es proporcional a ceff. Más rigurosamente, T_D^{-3}\propto c_{{\rm eff}}^{-3}:=(1/3)c_{{\rm long}}^{-3}+(2/3)c_{{\rm trans}}^{-3}, donde podemos distinguir las contribuciones longitudinal y transversal a la velocidad del sonido (1/3 y 2/3 respectivamente). La temperatura de Debye o la velocidad efectiva del sonido es una medida de la dureza del cristal.

Sustituyendo esto en la integral de la energía llegamos a

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

La facilidad con que se evalúan estas integrales para fotones se debe al hecho de que la frecuencia de la luz, al menos semiclásicamente, es independiente. Como ilustra la figura anterior, esto no es cierto para fonones. Para aproximar esta integral triple, Debye utilizó coordenadas esféricas

\ (n_x,n_y,n_z)=(n\cos \theta \cos \phi,n\cos \theta \sin \phi,n\sin \theta )

suponiendo valientemente que era lícito aproximar el cubo como una octava parte de esfera

U \approx\int_0^{\pi/2}\int_0^{\pi/2}\int_0^R E(n)\,{3\over e^{E(n)/kT}-1}n^2 \sin\theta\, dn\, d\theta\, d\phi\,,

donde R es el radio de esta esfera, que se halla mantieniendo invariante el número de partículaes en el cubo y en el octavo de esfera. El volumen del cubo es igual a N volúmenes unitarios,

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

y así llegamos a:

R = \sqrt[3]{6N\over\pi}\,.

La sustitución del dominio original por el esférico en la integral supone otra de las fuentes de error del modelo.

La integral de la energía se convierte en

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

haciendo el cambio el cambio de variable x = {hc_sn\over 2LkT},

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

Para simplificar el aspecto de esta expresión, definamos la temperatura de Debye TD (un resumen de algunas de las constantes y variables dependientes del material):

T_D\ \stackrel{\mathrm{def}}{=}\  {hc_sR\over2Lk} = {hc_s\over2Lk}\sqrt[3]{6N\over\pi} = {hc_s\over2k}\sqrt[3]{{6\over\pi}{N\over V}}

Llegamos así a la energía interna específica:

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

donde D3(x) es la (tercera) función de Debye.

Derivando con respecto a T obtenemos la capacidad calorífica adimensional:

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

Estas fórmulas tratan el modelo de Debye para todo el rango de temperaturas. Las fórmulas más elementales que se muestran más adelante se corresponden con el comportamiento asintótico en los límites de bajas y altas temperaturas. Como ya hemos mencionado, este comportamiento es exacto, a diferencia del comportamiento intermedio. La razón esencial para la exactitud en los rangos de bajas y altas temperaturas, respectivamente, es que el modelo de Debye da (i) la relación de dispersión correcta E(ν) para frecuencias bajas, y (ii) corresponde exactamente a la regla de suma (\int g(\nu ) \, {\rm d\nu}\equiv 3N)\,, sobre el número de vibraciones por intervalo de frecuencia.

Desarrollo de Debye

A decir verdad, Debye llegó a estos resultados de una manera ligeramente diferente y más simple. Utilizando la Mecánica de sólidos deformables encontró que el número de estados vibracionales con una frecuencia inferior a un cierto valor tendía asintóticamente a

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

donde V es el volumen y F es un factor que se calcula a partir de los coeficientes de elasticidad y densidad. Combinando esto con la energía esperable de un oscilador armónico a temperatura T (ya empleado anteriormente por Einstein en su modelo del sólido) daría una energía de

U = \int_0^\infty \,{h\nu^3 V F\over e^{h\nu/kT}-1}\, d\nu\,,

si las frecuencias vibracionales continuasen hasta el infinito. Esta fórmula da el exponente 3 de la temperatura, que es el comportamiento correcto a bajas temperaturas. Debye se dio cuenta de que no podía haber más de 3N estados vibracionales para N átomos. Lanzó la hipótesis de que en un sólido atómico el espectro de frecuencias de los estados vibracionales seguiría la ley anterior hasta una frecuencia máxima νm tal que el el número total de los estados fuese 3N:

 3N = {1 \over 3} \nu_m^3 V F \,.

Debye sabía que este supuesto no era realmente correcto (las frecuencias más altas están más estrechamente espaciadas de lo que él asumió), pero esto garantizaba un buen comportamiento para altas temperaturas (ley de Dulong-Petit). La energía venía dada entonces por:

U = \int_0^{\nu_m} \,{h\nu^3 V F\over e^{h\nu/kT}-1}\, d\nu\,,
 = V F kT (kT/h)^3 \int_0^{T_D/T} \,{x^3 \over e^x-1}\, dx\,,
donde TD es hνm / k.
 = 9 N k T (T/T_D)^3 \int_0^{T_D/T} \,{x^3 \over e^x-1}\, dx\,,
 = 3 N k T D_3(T_D/T)\,,

donde D3 es la función que luego se llamó función de Debye de tercer orden.

Límite de bajas temperaturas

La temperatura de un sólido de Debye se dice que es baja si T \ll T_D, lo que lleva a

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

Esta integral puede ser evaluada exactamente:

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

En el límite de bajas temperaturas, las limitaciones del modelo de Debye mencionadas anteriormente no se observan, y proporciona una relación correcta entre la capacidad calorífica (fonónica), la temperatura, los coeficientes elásticos y el volumen por átomo (las cantidades que están contenidas en la temperatura de Debye).

Límite de altas temperaturas

La temperatura de un sólido de Debye se dice que es alta si T > > TD. e^x - 1\approx  x si | x | < < 1, nos lleva a

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

Esta es la ley de Dulong-Petit, y es bastante exacta a pesar de que no tiene en cuenta la anarmonicidad, que causa que la capacidad calorífica continúe aumentando. Para describir la capacidad calorífica del sólido, si nos referimos a un conductor o a un semiconductor, deberíamos tener en cuenta también la nada desdeñable contribución de los electrones.

Debye contra Einstein

Debye vs. Einstein. Capacidad calorífica predicha como función de la temperatura.

¿Cuánto se ajustan realmente los modelos de Debye y Einstein a los hechos experimentales? Lo cierto es que sorprendentemente mucho, pero Debye gana en la región de bajas temperaturas, pues Einstein no predice con exactitud el comportamiento de los sólidos en estas condiciones.

¿Cómo de diferentes son los modelos? Para contestar a esta pregunta uno tendería simplemente a graficar ambos sobre el mismo conjunto de ejes... para darse cuenta a continuación de que no puede. Tanto el modelo de Einstein como el de Debye proveen de una forma funcional para la capacidad calorífica. Son modelos, y ningún modelo es tal sin una escala. Una escala relaciona el modelo con su contraparte en el mundo real. Uno puede ver que la escala del modelo de Einstein, dada por

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

es \epsilon/k. Y la escala del modelo de Debye es TD, la temperatura de Debye. Ambos suelen determinarse ajustando los modelos a los datos experimentales (La temperatura de Debye puede ser calculada teóricamente a partir de la velocidad del sonido y de las dimensiones del cristal). Dado que ambos métodos se aproximan al problema desde distintas direcciones y diferentes geometrías, las escalas de Debye y Einstein no son la misma, es decir,

{\epsilon\over k} \ne T_D\,,

lo que significa que carece de sentido representarlas sobre el mismo conjunto de ejes coordenados. Ambos son modelos de lo mismo, pero sus escalas son distintas. Si definiésemos una temperatura de Einstein como

T_E \ \stackrel{\mathrm{def}}{=}\  {\epsilon\over k}\,,

entonces podríamos decir

T_E \ne T_D\,,

y, para relacionar ambas, deberíamos buscar el ratio

\frac{T_E}{ T_D} = ?

El sólido de Einstein está compuesto por osciladores armónicos cuánticos de frecuencia única \epsilon = \hbar\omega = h\nu. Esta frecuencia, si de hecho existiese, estaría relacionada con la velocidad del sonido en el sólido. Si uno imagina la propagación del sonido como una secuencia de átomos chocando los unos con los otros, se convierte en obvio que la frecuencia de oscilación debe de hecho corresponderse con la longitud de onda mínima soportable por la red, λmin.

\nu = {c_s\over\lambda} = {c_s\sqrt[3]{N}\over 2L} = {c_s\over 2}\sqrt[3]{N\over V}

lo que hace la temperatura de Einstein

T_E = {\epsilon\over k} = {h\nu\over k} = {h c_s\over 2k}\sqrt[3]{N\over V}\,,

y la relación buscada resulta

{T_E\over T_D} = \sqrt[3]{\pi\over6}\,.

Ahora ambos modelos pueden ser comparados sobre la misma gráfica. Nótese que este coeficiente es la raíz cúbica del cociente del volumen de un octante de una esfera tridimensional con el volumen del cubo que la contiene, que es justo el factor de correción empleado por Debye para aproximar la integral de la energía.

Tabla de temperaturas de Debye

Aún cuando el modelo de Debye no es completamente correcto, da una buena aproximación para la capacidad calorífica a bajas temperaturas de los sólidos cristalinos aislantes, donde otras contribuciones (como la debida a los electrones de la banda de conducción) son despreciables. Para metales, la contribución al calor específico de los electrones es proporcional a T, lo que hace que para bajas temperaturas esta contribución domine al T3 predicho por Debye como resultado de las oscilaciones de la red. En este caso, el modelo de Debye solo puede emplearse para aproximar la contribución de la red al calor específico. La siguiente tabla muestra las temperaturas de Debye para varias sustancias:[2]

Aluminio 428 K
Cadmio 209 K
Cromo 630 K
Cobre 343.5 K
Oro 165 K
Hierro 470 K
Plomo 105 K
Manganeso 410 K
Níquel 450 K
Platino 240 K
Silicio 645 K
Plata 225 K
Tántalo 240 K
Estaño (blanco) 200 K
Titanio 420 K
Wolframio 400 K
Zinc 327 K
Carbono 2230 K
Hielo 192 K


Véase también

  • Gas de Bose
  • gas en una caja

Referencias

  1. 'Zur Theorie der spezifischen Waerme', Annalen der Physik (Leipzig) 39(4), p. 789 (1912)
  2. Kittel, Charles, Introduction to Solid State Physics, 7th Ed., Wiley, (1996)
  • CRC Handbook of Chemistry and Physics, 56th Edition (1975-1976)
  • Schroeder, Daniel V. An Introduction to Thermal Physics. Addison-Wesley, San Francisco, Calif. (2000). Section 7.5.
  • Kittel, Charles, Introduction to Solid State Physics, 7th Ed., Wiley, (1996)

Enlaces externos


Wikimedia foundation. 2010.

Игры ⚽ Поможем сделать НИР

Mira otros diccionarios:

  • Modelo de iones móviles — El modelo de iones móviles es un modelo para la reacción de las moléculas en donde se tiene en cuenta el potencial químico y el campo eléctrico; se pueden representar los solventes iónicos de forma explicita o no. Representaciones más utilizadas… …   Wikipedia Español

  • Longitud de Debye — En Física de plasmas la longitud de Debye, también llamada radio de Debye es la escala a través de la cual portadores de carga móviles por ejemplo, electrones generan un apantallamiento de los campos eléctricos en los plasmas y otros conductores …   Wikipedia Español

  • Peter Debye — Saltar a navegación, búsqueda Peter Debye (1912) Petrus (Peter) Josephus Wilhelmus Debye o Petrus Josephus Wilhelmus Debije (n. Maastricht, 24 de marzo de 1884 † Ithaca, Nueva York, 2 de noviembre de 1966) fue …   Wikipedia Español

  • Ecuación de Debye-Huckel — Saltar a navegación, búsqueda Peter Debye y Erich Hückel plantearon un modelo que describe una mezcla de iones (electrolitos), inmersos en un medio dieléctrico continuo de temperatura T, presión P y concentración molar ci, , en donde los iones de …   Wikipedia Español

  • Ley de Dulong-Petit — Figura 1. Dependencia del calor específico atómico de elementos sólidos con la temperatura …   Wikipedia Español

  • Capacidad calorífica — La capacidad calorífica de un cuerpo es el cociente entre la cantidad de energía calorífica transferida a un cuerpo o sistema en un proceso cualquiera y el cambio de temperatura que experimenta. En una forma menos formal es la energía necesaria… …   Wikipedia Español

  • Par de Cooper — Saltar a navegación, búsqueda En superconductividad se conoce como par de Cooper a la pareja de electrones que se hallan ligados debido a que en el estado superconductor ambas partículas se comportan como si se atrayeran, pese a tener ambos carga …   Wikipedia Español

  • Física del estado sólido — La física del estado sólido, rama de la física de la materia condensada, trata sobre el estudio de la materia rígida, o sólidos. Estudia las propiedades físicas de los materiales sólidos utilizando disciplinas tales como la mecánica cuántica, la… …   Wikipedia Español

  • Universidad Técnica de Aquisgrán — RWTH Aachen University Rheinisch Westfälische Technische Hochschule Aachen Lema «Zukunft beginnt bei uns[1] » « El futuro empieza con nosotros » Tipo Pública …   Wikipedia Español

  • Fotón — Saltar a navegación, búsqueda Para otros usos de este término, véase Fotón (desambiguación). Fotón (γ) Fotones emitidos en un rayo coherente por un láser …   Wikipedia Español

Compartir el artículo y extractos

Link directo
Do a right-click on the link above
and select “Copy Link”