Un mol de sólido son N_A átomos = 3N_A modos normales (el censo del artículo 03). La física clásica reparte el calor democráticamente — la equipartición asigna k_BT por oscilador (½ cinética + ½ potencial) — y produce una predicción sin parámetros:
Es la ley de Dulong-Petit (1819, empírica un siglo antes de explicarse): plomo 26,5, cobre 24,5, oro 25,4 — elementos sin nada en común clavando el mismo número, porque el número no es química: es un censo. Y luego está el diamante: 6,1 a temperatura ambiente. Un 75 % de calor desaparecido — el escándalo que este artículo liquida.
El diamante y la solución de Einstein
Problema. ¿Por qué el diamante incumple Dulong-Petit a 300 K — y el plomo no?
Solución. III.3 lo dejó dicho: un oscilador cuántico solo acepta energía en escalones ħω. Si ħω ≫ k_BT, el escalón no se puede pagar y el modo queda congelado — censado pero sin cobrar. El diamante es el caso extremo: átomos ligeros y muelles ferocísimos (el enlace C–C más rígido de la química) suben sus frecuencias hasta ~40 THz — escalones de ħω ≈ k_B×(2000 K): a temperatura ambiente, la mayoría de sus modos están cerrados por caros. El plomo — pesado y blando, modos de ~1 THz — tiene toda la plantilla activa desde −200 °C.
Resultado. Einstein publicó esta explicación en 1907: la segunda aplicación de los cuantos de la historia (tras el fotón), y la primera a la materia. Su modelo (todos los modos a una frecuencia) clavó la caída cualitativa; Debye la afinó en 1912 usando la dispersión del artículo 03 (los modos acústicos graves siempre quedan activos: de ahí el famoso C ∝ T³ a baja temperatura). El calor específico de un cristal es el censo de la cadena pasado por la taquilla de ħω — y fue una de las pruebas decisivas de que los cuantos no eran una rareza de la luz. El Termario heredará la taquilla completa.
El módulo, cerrado — y la escalera entera
Mira la escalera que acabas de subir. Dos carritos con álgebra (01) → una molécula con espectrómetro (02) → una cadena con dispersión y fonones (03) → un cristal con termodinámica (04). Es la misma matriz en cuatro tamaños — y la ruta histórica real por la que la mecánica de Lagrange desembocó en la física del estado sólido. Los cabos que quedan sueltos tienen dueño: la estadística de los modos (Termario), sus cuantos en serio (Cuantario), y el cristal completo con electrones (Cristalario). Este módulo es la estación de la que salen esos tres trenes. En Mecanario quedan dos paradas: el caos cuantitativo — donde la linealización de este módulo revienta — y los fluidos.
Ejercicios
La regla del gremio: un modo está activo si T supera aproximadamente ħω/k_B (su «temperatura característica»). Estima esa temperatura para los modos del hierro (~6 THz máx.) y del diamante (~40 THz), y predice quién cumple Dulong-Petit a temperatura ambiente.
Solución
ħω/k_B ≈ (48 K por THz): hierro ~290 K — justo activo a ambiente (C(Fe) = 25,1: cumple por los pelos); diamante ~1900 K — congelado hasta el rojo vivo (alcanza 3R por encima de ~1500 K, y se midió: la curva completa del diamante fue de las primeras verificaciones de Einstein-Debye). La regla «48 K por THz» — o su gemela óptica «un modo se ve en IR donde se congela en calor» — cose espectroscopía y termodinámica con un solo número. Guárdala: el Termario la usará sin avisar.
La equipartición aplicada a los gases explica un viejo misterio escolar: el aire (diatómico) tiene C_v = 5R/2 y no 3R/2 (monoatómico) ni 7R/2 (con vibración). Haz el censo de la molécula de N₂ a temperatura ambiente.
Solución
Tres traslaciones (½kT cada una) + dos rotaciones (los dos ejes transversales; girar sobre el propio eje de la molécula no cuenta — su momento de inercia es ridículo y su escalón, inalcanzable) = 5 términos: 5R/2 ✓. La vibración del enlace N≡N (2360 cm⁻¹ ≈ 3400 K) está congelada: sus 2 términos (cinético + potencial) no cobran hasta miles de kelvin. El calor específico del aire que respiras lleva impresa la taquilla cuántica — y fue OTRO de los escándalos (Maxwell ya lo señaló en 1875: «aquí hay algo que no entendemos») que solo 1907 resolvió. Los grados de libertad no se cuentan: se cuentan y se tarifican.
¿Por qué el helio líquido, el mejor conductor de calor conocido a baja temperatura, y el diamante, el mejor a temperatura ambiente, deben ambos su récord a los fonones — por caminos opuestos?
Solución
El calor de un aislante viaja en fonones, y conduce mejor quien los transporta más lejos sin chocar. Diamante: fonones rapidísimos (v ~ 18 000 m/s: muelles brutales, átomos ligeros — la fórmula del artículo 03) y red tan perfecta y anarmónica tan débil que vuelan micras sin dispersarse: ~2000 W/m·K, cinco veces el cobre — por eso los disipadores de lujo y los yunques de células de diamante. El helio superfluido va más allá: sus excitaciones viajan sin fricción ninguna (el «segundo sonido» — el calor propagándose como onda, no difundiéndose). Dos récords, una moraleja de cierre: el transporte de calor es dinámica de modos normales — la matriz del artículo 01, en su oficio más cotidiano.
Todo lo que el módulo deja utilizable, con dónde se dedujo cada cosa — y, cuando toca, qué deuda de un módulo anterior salda cada fila. Pensado para leerse dentro de seis meses sin releer nada: si para resolver un problema típico del módulo hay que volver al texto a buscar una constante, esta tabla ha fallado.
| Qué | Fórmula o valor | Dónde |
|---|---|---|
| Cerca del equilibrio, todo es cuadrático | T = ½ẋᵀTẋ y V = ½xᵀVx, con ambas matrices simétricas; V es el hessiano del potencial | art. 01 · las curvaturas de II.2, ahora en tabla |
| Problema de autovalores | V·a = ω²·T·a. Cada autovalor es una frecuencia al cuadrado; cada autovector, la forma del modo: quién se mueve, cuánto y con qué signo | art. 01 |
| Diagonalizar es resolver | En las coordenadas de los autovectores el sistema se desacopla en osciladores independientes | art. 01 · mismo teorema espectral que el tensor de III.4 |
| Los dos carritos, con matrices | T = m·1 y V con (k+k′) en la diagonal y −k′ fuera ⟹ ω₁² = k/m con a₁ = (1,1) y ω₂² = (k+2k′)/m con a₂ = (1,−1). Con k′ = k: ω₂/ω₁ = √3 = 1,732 | art. 01, ej. res. 1 · industrializa el «a ojo» de II.3 |
| Dos regalos del teorema espectral | Autovalores reales y autovectores ortogonales — la descomposición en modos es única (el Fourier de II.7 era su corolario). Y si un autovalor sale negativo, ω es imaginaria: el equilibrio era silla, no valle | art. 01 |
| El mismo cálculo es test de estabilidad | Versión matricial del criterio U″ > 0 de II.2; los modos de pandeo de una estructura son autovalores que cruzan el cero | art. 01 |
| El algoritmo del gremio | (1) encontrar el equilibrio · (2) expandir a segundo orden: matrices T y V · (3) diagonalizar · (4) escribir el movimiento como suma de modos con amplitudes de las condiciones iniciales | art. 01, callout |
| Péndulo doble linealizado | ω± = √((2 ∓ √2)·g/L). Lento: ambos al unísono. Rápido: en contrafase, con casi todo el movimiento en el segundo péndulo. El caos de I.5 vive fuera de esta aproximación | art. 01, ej. 1 |
| Ortogonalidad, con su producto escalar | Con masas distintas los autovectores se tuercen y la ortogonalidad prometida es respecto de la métrica T (aᵀTb = 0), no la euclídea | art. 01, ej. 3 |
| Los modos de la Tierra | Tras el terremoto de 2004 el planeta sonó semanas: ₀S₂ («balón de rugby») 53,9 minutos y el respirador radial 20,5. Se miden con gravímetros superconductores y pesan el núcleo | art. 01, ej. 2 |
| Censo de modos de vibración | 3N − 6 en general, 3N − 5 si la molécula es lineal. CO₂: 3×3 − 5 = 4 modos — dos estiramientos y una flexión doblemente degenerada | art. 02 |
| CO₂ longitudinal, calculado | Traslación (ω = 0) · simétrico ω² = k/m_O (el carbono no se entera) · antisimétrico ω² = k(1/m_O + 2/m_C). Cociente predicho √(1 + 2m_O/m_C) = √(1+32/12) = 1,91 | art. 02, ej. res. 1 |
| El veredicto del espectrómetro | 2349 cm⁻¹ / 1388 cm⁻¹ = 1,69. El desajuste mide química: estirar un enlace C=O cambia la rigidez del otro, y hace falta una constante de acoplo (campo de fuerzas de valencia) | art. 02, ej. res. 1 |
| Huella dactilar infrarroja | O–H ~3600 cm⁻¹ · C–H ~2900 · C=O ~1700; la región 500–1500 cm⁻¹ es la «huella dactilar». Cada pico es un autovalor y su intensidad, cuánto cambia el reparto de cargas ese autovector | art. 02 |
| Cuentas de modos | H₂O 3 · NH₃ 6 · CH₄ 9 · benceno 30. La simetría agrupa modos en degenerados y silencia en infrarrojo los que no mueven cargas: cuanto más simétrica, más muda | art. 02, ej. 1 |
| Marcado isotópico | ω = √(k/μ) con μ dominada por el átomo ligero: el estiramiento O–H de 3600 cm⁻¹ baja a ~2620 en O–D (factor 0,73). El enlace no cambia; solo la masa | art. 02, ej. 2 |
| Por qué N₂ y O₂ no calientan | Diatómica homonuclear: un único modo que por simetría jamás altera el reparto de cargas — infrarrojo-mudo. El agua, doblada, tiene 3 modos y todos activos: es el invernadero dominante | art. 02, ej. 3 y callout |
| Relación de dispersión de la cadena | ω(k) = 2√(K/m)·|sen(ka/2)|, para masas m unidas por muelles K a distancia a | art. 03 |
| Límite de onda larga | ka ≪ 1 ⟹ ω ≈ a√(K/m)·k: lineal, todas las ondas viajan a v = a√(K/m). Esa v es el sonido del material, y el medio se comporta como la cuerda continua de II.7 | art. 03 |
| Borde de Brillouin | k = π/a: vecinos en contrafase, λ = 2a — la onda más corta que la red puede representar. La velocidad de grupo v_g = a√(K/m)·cos(ka/2) se anula ahí (condición de Bragg en 1D) | art. 03, ej. 1 |
| El sonido del hierro, desde los átomos | a = 2,5 Å, m = 56 u = 9,3×10⁻²⁶ kg, K ≈ 30 N/m ⟹ v = a√(K/m) ≈ 4500 m/s (real ~5100: 10 % de error) y frecuencia máxima de la red ~6 THz | art. 03, ej. res. 1 · cierra el raíl del western de I.4 |
| Fonones | Modo normal de la cadena + cuanto de acción: las vibraciones llegan en paquetes ħω(k). Transportan el calor de los aislantes, frenan electrones y emparejan pares en los superconductores | art. 03 · aplica los escalones de III.3 |
| Cadena diatómica | Dos ramas con un hueco entre ellas: acústica (las dos masas a la par, ω → 0) y óptica (m contra M dentro de la celda, alta incluso a k = 0). La sal absorbe a ~5 THz por su rama óptica | art. 03, ej. 2 |
| Censo de un sólido | N masas, exactamente N modos; en tres dimensiones, 3N para N átomos. Número finito y conocido de osciladores donde repartir el calor | art. 03, ej. 3 |
| Ley de Dulong-Petit (1819) | C = 3N_A·k_B = 3R ≈ 25 J/(mol·K) para todo sólido. Plomo 26,5 · cobre 24,5 · oro 25,4 · hierro 25,1 — el número no es química, es censo | art. 04 |
| El escándalo del diamante | 6,1 J/(mol·K) a 300 K: falta el 75 % del calor previsto | art. 04 |
| Modos congelados (Einstein, 1907) | Si ħω ≫ k_BT el escalón no se puede pagar y el modo queda censado pero sin cobrar. Debye (1912) lo afinó con la dispersión del artículo 03: de ahí el C ∝ T³ a baja temperatura | art. 04, ej. res. 1 |
| La regla de conversión | ħω/k_B ≈ 48 K por THz. Hierro (~6 THz) → ~290 K: activo a temperatura ambiente por los pelos. Diamante (~40 THz) → ~1900 K: congelado hasta el rojo vivo | art. 04, ej. 1 |
| Equipartición en gases | N₂ a temperatura ambiente: 3 traslaciones + 2 rotaciones ⟹ C_v = 5R/2. La vibración del enlace (2360 cm⁻¹ ≈ 3400 K) está congelada y no cobra sus dos términos | art. 04, ej. 2 |
| Conducción por fonones | Diamante: v ~18 000 m/s y red casi perfecta ⟹ ~2000 W/(m·K), cinco veces el cobre. Helio superfluido: excitaciones sin fricción ninguna — el «segundo sonido», calor que se propaga como onda | art. 04, ej. 3 |