Physics · Libro 5 · Bachelor Year 3

Física universitaria — tercer año

Física universitaria — tercer año · Bachelor Year 3

18El colectivo macrocanónico

Una gota de agua en aire húmedo ni crece ni se encoge; el oxígeno se une a la hemoglobina en los pulmones y se suelta en los músculos; una pila empuja electrones a través de un teléfono porque son “más caros” en un electrodo que en el otro. En cada caso, la moneda que se equilibra no es la temperatura, sino el potencial químico μ\mu: el coste en energía libre de una partícula más. Este capítulo abre la segunda válvula del reservorio — partículas además de energía — y construye el colectivo macrocanónico, cuya recompensa es una eficacia asombrosa: la ocupación de un único nivel cuántico, sea fermiónico o bosónico, se sigue en dos líneas, y con ella las fórmulas maestras del capítulo siguiente. Sus aplicaciones van de las isotermas de adsorción de las máscaras antigás y los catalizadores a la ecuación de Saha, con la que una estudiante de doctorado de Harvard leyó en 1925 las temperaturas de las estrellas y descubrió, contra toda expectativa, que el universo está hecho de hidrógeno.

18.1 Las partículas entran en la negociación

Teorema 18.1 (Distribución macrocanónica)

Un sistema que intercambia energía y partículas con un reservorio a temperatura TT y potencial químico μ\mu ocupa su estado ss (de energía EsE_s y número de partículas NsN_s) con probabilidad

Ps=eβ(EsμNs)Ξ,Ξ=seβ(EsμNs),\mathcal P_s = \frac{\eu^{-\beta(E_s - \mu N_s)}}{\Xi} , \qquad \Xi = \sum_s\eu^{-\beta(E_s - \mu N_s)} ,

la función de partición macrocanónica. De Ξ\Xi: N=kBTlnΞ/μ\langle N\rangle = k_{\text{B}}T\,\partial\ln\Xi/\partial\mu, y el potencial macrocanónico Φ=kBTlnΞ\Phi = -k_{\text{B}}T\ln\Xi cumple Φ=PV\Phi = -PV.

Demostración parcial. Como en el Teorema 17.1, desarróllese la entropía del reservorio en las dos retiradas: Sres(EEs,NNs)SresEs/T+μNs/TS_{\text{res}}(E - E_s, N - N_s) \approx S_{\text{res}} - E_s/T + \mu N_s/T, usando S/N=μ/T\partial S/\partial N = -\mu/T (Teorema 16.6). Los promedios se siguen por derivación; Φ=PV\Phi = -PV se admite (se sigue de la extensividad).

Definición 18.2 (Qué significa μ\mu)

μ\mu es el coste, en energía libre, de añadir una partícula a TT y VV fijos. Las partículas fluyen espontáneamente de μ\mu alto a μ\mu bajo — la difusión, la ósmosis, la evaporación, las reacciones químicas y las corrientes eléctricas son todas gradientes de μ\mu que se relajan. En el equilibrio entre dos fases o dos lugares cualesquiera, μ\mu es igual a ambos lados — la tercera igualdad (junto con TT y PP) que gobierna la coexistencia. Para un gas ideal clásico,

μ=kBTln ⁣(nλT3)<0(gas diluido):\mu = k_{\text{B}}T\,\ln\!\big(n\lambda_T^3\big) < 0 \quad\text{(gas diluido)} :

añadir una partícula a un gas diluido baja la energía libre, porque la ganancia de entropía supera cualquier coste energético — el signo que impulsa la evaporación hacia el aire seco. Y una pila es una máquina de μ\mu: su tensión es la caída de potencial químico por electrón entre sus electrodos, U=Δμ/eU = \Delta\mu/e — los voltios son el tipo de cambio de la química del electrón, y por eso las pilas entregan voltajes de una sola cifra.

18.2 Un nivel cada vez: el truco maestro

Teorema 18.3 (Ocupación de un solo nivel cuántico)

Tómese como “el sistema” un único nivel de una partícula, de energía ϵ\epsilon, en contacto macrocanónico con el gas de todos los demás. Fermiones (ocupación 00 o 11):

Ξ=1+eβ(ϵμ)    nˉFD=1eβ(ϵμ)+1.\Xi = 1 + \eu^{-\beta(\epsilon - \mu)} \;\Longrightarrow\; \bar n_{\text{FD}} = \frac{1}{\eu^{\beta(\epsilon - \mu)} + 1} .

Bosones (ocupación 0,1,2,0, 1, 2, \dots):

Ξ=11eβ(ϵμ)    nˉBE=1eβ(ϵμ)1,\Xi = \frac{1}{1 - \eu^{-\beta(\epsilon - \mu)}} \;\Longrightarrow\; \bar n_{\text{BE}} = \frac{1}{\eu^{\beta(\epsilon - \mu)} - 1} ,

con la exigencia μ<ϵmin\mu < \epsilon_{\min}. Ambas se reducen al factor de Boltzmann eβ(ϵμ)\eu^{-\beta(\epsilon - \mu)} cuando las ocupaciones son escasas. Estas dos fórmulas, obtenidas en dos líneas cada una, son todo el material de partida de la estadística cuántica: el capítulo siguiente es su cosecha.

Demostración. Fermiones: dos términos. Bosones: una serie geométrica, convergente solo para μ<ϵ\mu < \epsilon. En ambos casos, nˉ=kBTlnΞ/μ\bar n = k_{\text{B}}T\, \partial\ln\Xi/\partial\mu da las formas enunciadas.

Las tres leyes de ocupación frente a ( - )/k_ BT: los fermiones saturan en uno por estado (un escalón suavizado en =), los bosones se amontonan sin límite cuando , y ambas se unen a la cola diluida de Boltzmann por la derecha — el régimen clásico es aquel en el que un estado apenas llega a visitarse.
Las tres leyes de ocupación frente a (ϵμ)/kBT(\epsilon - \mu)/k_{\text{B}}T: los fermiones saturan en uno por estado (un escalón suavizado en ϵ=μ\epsilon = \mu), los bosones se amontonan sin límite cuando ϵμ\epsilon \to \mu, y ambas se unen a la cola diluida de Boltzmann por la derecha — el régimen clásico es aquel en el que un estado apenas llega a visitarse.

18.3 Adsorción: sitios en alquiler

Proposición 18.4 (La isoterma de Langmuir)

Una superficie ofrece sitios independientes, cada uno vacío o con una molécula de gas de energía de enlace ϵ0-\epsilon_0; el gas de encima, a presión PP, es el reservorio. Cada sitio es un sistema macrocanónico de dos estados, y su ocupación es

θ=1eβ(ϵ0μ)+1=PP+P0(T),P0(T)=kBTλT3eϵ0/kBT,\theta = \frac{1}{\eu^{\beta(-\epsilon_0 - \mu)} + 1} = \frac{P}{P + P_0(T)} , \qquad P_0(T) = \frac{k_{\text{B}}T}{\lambda_T^3}\, \eu^{-\epsilon_0/k_{\text{B}}T} ,

usando el μ\mu del gas ideal. La cobertura crece linealmente a baja presión y satura en una monocapa — la isoterma de Langmuir, la curva de trabajo de los catalizadores de automóvil, el carbón activo, las máscaras antigás y la fijación del oxígeno por la sangre (Ejercicio 18.6). La presión de cruce P0P_0 cae exponencialmente con la energía de enlace y sube de forma abrupta con la temperatura: calentar una superficie sacude a sus huéspedes hasta soltarlos — por eso hornear el material de vidrio al vacío es la única manera de limpiarlo de verdad.

Demostración. La ocupación del sitio es la del Teorema 18.3 (en forma fermiónica: un sitio alberga como mucho uno) con ϵ=ϵ0\epsilon = -\epsilon_0; sustitúyase eβμ=nλT3=PλT3/kBT\eu^{\beta\mu} = n\lambda_T^3 = P\lambda_T^3/k_{\text{B}}T.

Isotermas de Langmuir: la cobertura frente a la presión, saturando en una monocapa. Caliéntese la superficie y la misma cobertura exigirá mucha más presión — la adsorción es una negociación de  que la superficie caliente no deja de perder.
Isotermas de Langmuir: la cobertura frente a la presión, saturando en una monocapa. Caliéntese la superficie y la misma cobertura exigirá mucha más presión — la adsorción es una negociación de μ\mu que la superficie caliente no deja de perder.

18.4 Equilibrios de ionización: la ecuación de Saha

Proposición 18.5 (La ecuación de Saha)

En un gas caliente, la reacción Hp+e\text{H} \rightleftharpoons \text{p} + \text{e} se equilibra cuando μH=μp+μe\mu_{\text{H}} = \mu_{\text{p}} + \mu_{\text{e}}. Escribir cada μ\mu con la forma del gas ideal (con el átomo acreditado con su energía de enlace EI-E_{\text{I}}) da

nenpnH=1λT,e3eEI/kBT(λT,e=h/2πmekBT)\frac{n_{\text{e}}\,n_{\text{p}}}{n_{\text{H}}} = \frac{1}{\lambda_{T,\text{e}}^3}\, \eu^{-E_{\text{I}}/k_{\text{B}}T} \qquad\Big(\lambda_{T,\text{e}} = h/\sqrt{2\pi m_{\text{e}}k_{\text{B}}T}\Big)

(salvo factores de degeneración de espín del orden de la unidad). De nuevo dos monedas que compiten: el factor de Boltzmann grava la ionización con EI=13.6eVE_{\text{I}} = 13.6\,\mathrm{eV}, mientras que el prefactor — la entropía del electrón liberado, con sus λT3\lambda_T^{-3} estados accesibles — la subvenciona. La subvención es tan grande a las densidades estelares que el hidrógeno se ioniza apreciablemente ya cerca de los 10000K10\,000\,\mathrm{K}, donde kBTk_{\text{B}}T vale solo 0.9eV0.9\,\mathrm{eV}: los equilibrios se disputan en energía libre, no en energía — y los espectros de las estrellas son el libro de cuentas (Problema 18.1).

Demostración parcial. Hágase μH=μp+μe\mu_{\text{H}} = \mu_{\text{p}} + \mu_{\text{e}} con μi=kBTln(niλT,i3)+εi\mu_i = k_{\text{B}}T\ln(n_i\lambda_{T,i}^3) + \varepsilon_i (εH=EI\varepsilon_{\text{H}} = -E_{\text{I}}, y los demás nulos); exponénciese. Las longitudes de onda térmicas del protón y del átomo casi se cancelan (masas iguales); los factores de degeneración se admiten.

La ionización de Saha del hidrógeno a densidad fotosférica (esquemática): una transición brusca cerca de los 10\,000\, K, muy por debajo de los ingenuos E_ I/k_ B 160\,000\, K — la entropía del electrón liberado paga la mayor parte de la factura.
La ionización de Saha del hidrógeno a densidad fotosférica (esquemática): una transición brusca cerca de los 10000K10\,000\,\mathrm{K}, muy por debajo de los ingenuos EI/kB160000KE_{\text{I}}/k_{\text{B}} \approx 160\,000\,\mathrm{K} — la entropía del electrón liberado paga la mayor parte de la factura.

Método 18.6 (El oficio macrocanónico)

(1) Todo lo que intercambie partículas con un socio grande: dese al socio un μ\mu y pésense los estados con eβ(EμN)\eu^{-\beta(E - \mu N)}. (2) Los sitios o niveles independientes se factorizan: resuélvase uno y multiplíquese. (3) Para un reservorio de gas clásico, eβμ=nλT3\eu^{\beta\mu} = n\lambda_T^3 convierte μ\mu en presión medible. (4) Reacciones y cambios de fase: equilíbrense los μ\mu, uno por especie, con las energías de enlace incluidas — espérese que los prefactores entrópicos desplacen los equilibrios muy lejos de las conjeturas ingenuas de Boltzmann. (5) Electrones en sólidos: su μ\mu es el nivel de Fermi, y las diferencias de μ\mu son tensiones — el puente hacia el Capítulo 24.

18.5 Ejercicios

Ejercicio 18.1

(a) Se conectan dos cajas del mismo gas a la misma TT pero con densidades distintas: ¿en qué sentido fluyen las partículas, en lenguaje de μ\mu? (b) Muéstrese, a partir del μ\mu del gas ideal, que el flujo cesa con densidades iguales. (c) Hay agua bajo aire húmedo: enúnciese la condición de equilibrio en términos de μ\mu. (d) ¿Por qué se seca la ropa incluso por debajo de 100100\,^\circC cuando el aire está seco (compárense los μ\mu, no las temperaturas)?

Solución

Solución de Ejercicio 18.1.

(a) De la caja más densa a la menos densa: de μ=kBTln(nλT3)\mu = k_{\text{B}}T\ln(n\lambda_T^3) alto a bajo. (b) Con la misma TT, λT\lambda_T es común: igual μ\mu es igual nn. (c) μlıˊquido=μvapor\mu_{\text{líquido}} = \mu_{\text{vapor}} — la definición de humedad saturante. (d) En aire seco, el μ\mu del vapor (con nn bajo) queda por debajo del del líquido: el agua se evapora cuesta abajo en μ\mu por más que nunca hierva.

Ejercicio 18.2

El potencial químico del aire. (a) Calcúlese λT\lambda_T para el N2_2 a 300K300\,\mathrm{K}. (b) Calcúlese nλT3n\lambda_T^3 a 1bar1\,\mathrm{bar} y el μ\mu resultante en eV. (c) ¿Por qué no es paradójico μ<0\mu < 0 (qué término de F=ETSF = E - TS domina al añadir una molécula)? (d) ¿A qué densidad alcanzaría μ\mu el valor cero, y qué anuncia nλT31n\lambda_T^3 \sim 1 (Ejercicio 18.11)?

Solución

Solución de Ejercicio 18.2.

(a) λT=1.9×1011m\lambda_T = 1.9 \times 10^{-11}\,\mathrm{m}. (b) nλT3=1.7×107n\lambda_T^3 = 1.7 \times 10^{-7}: μ=kBTln(1.7×107)=0.40eV\mu = k_{\text{B}}T\ln(1.7 \times 10^{-7}) = -0.40\,\mathrm{eV}. (c) El término entrópico: una molécula nueva abre 1/nλT3\sim1/n\lambda_T^3 estados, y TS-TS gana a cualquier coste cinético. (d) n=λT31.5×1032m3n = \lambda_T^{-3} \approx 1.5 \times 10^{32}\,\mathrm{m}^{-3} — miles de veces la densidad del líquido para el aire; nλT31n\lambda_T^3 \sim 1 es la frontera de la degeneración cuántica (Ejercicio 18.11).

Ejercicio 18.3

Dedúzcase la ocupación de Fermi–Dirac a partir de su Ξ\Xi de dos términos, paso a paso; compruébense los límites ϵμ\epsilon \ll \mu y ϵμ\epsilon \gg \mu, y el valor en ϵ=μ\epsilon = \mu. ¿Dónde ha aparecido ya esta función en este libro (Ejercicio 16.8)?

Solución

Solución de Ejercicio 18.3.

Ξ=1+eβ(ϵμ)\Xi = 1 + \eu^{-\beta(\epsilon-\mu)}; nˉ=kBTμlnΞ=eβ(ϵμ)/(1+eβ(ϵμ))\bar n = k_{\text{B}}T\, \partial_\mu\ln\Xi = \eu^{-\beta(\epsilon-\mu)}/(1 + \eu^{-\beta(\epsilon-\mu)}), la forma enunciada. Límites: 11 muy por debajo de μ\mu, cola de Boltzmann muy por encima, y exactamente 12\tfrac12 en ϵ=μ\epsilon = \mu. Es la ocupación térmica de dos niveles del sólido de Schottky — la misma curva, con otra etiqueta.

Ejercicio 18.4

Los voltios son química. (a) Una pila de litio da 3.0V3.0\,\mathrm{V}: exprésese la diferencia de potencial químico por electrón en eV y por mol en kJ. (b) ¿Por qué tienen todas las pilas cotidianas voltajes de una sola cifra (qué fija la escala de las diferencias químicas de μ\mu)? (c) Una batería de plomo y ácido de coche apila seis celdas: ¿por qué apilar en vez de buscar una reacción de 12eV12\,\mathrm{eV}? (d) En lenguaje de μ\mu, ¿qué hace una pila cuando está “gastada”?

Solución

Solución de Ejercicio 18.4.

(a) 3.0eV3.0\,\mathrm{eV} por electrón; ×NAe\times N_{\text{A}}e: 290kJ/mol\approx290\,\mathrm{kJ}/\mathrm{mol}. (b) Los potenciales químicos difieren en energías de enlace — electronvoltios: las pilas heredan la escala de energías de la química. (c) Ningún par químico ofrece un escalón de 12eV12\,\mathrm{eV} (destrozaría el electrolito); apilar en serie suma caídas de μ\mu sin peligro. (d) Los potenciales químicos electrónicos de sus electrodos se han igualado: sin gradiente no hay empuje, por mucha materia que quede.

Ejercicio 18.5 ★★

Bosones en un nivel. (a) Súmese la serie geométrica de Ξ\Xi y dedúzcase nˉBE\bar n_{\text{BE}}. (b) ¿Por qué debe quedarse μ\mu por debajo del nivel más bajo — qué diverge si no, y qué anticipa esa divergencia (Capítulo 19)? (c) Calcúlese nˉ\bar n para (ϵμ)/kBT=0.1(\epsilon - \mu)/k_{\text{B}}T = 0.1, 11 y 33, y compárese con Boltzmann. (d) Los fotones tienen μ=0\mu = 0: dese la razón física (su número no se conserva — las paredes los absorben y los emiten sin trabas).

Solución

Solución de Ejercicio 18.5.

(a) Serie geométrica de razón eβ(ϵμ)<1\eu^{-\beta(\epsilon-\mu)} < 1; derívese. (b) Cuando μϵmin\mu \to \epsilon_{\min}, la serie diverge: el nivel puede absorber población sin límite — el amontonamiento macroscópico que se convierte en la condensación de Bose–Einstein. (c) nˉ=9.51\bar n = 9.51, 0.580.58 y 0.0520.052 frente a los 0.900.90, 0.370.37 y 0.0500.050 de Boltzmann: idénticos en la cola diluida y enormemente realzados cerca de μ\mu. (d) Como las paredes crean y destruyen fotones a su antojo, el equilibrio no puede depender de su número: el coste en energía libre de un fotón más es nulo, μ=0\mu = 0.

Ejercicio 18.6 ★★

Langmuir en acción. (a) A partir de la isoterma, ¿qué presión logra θ=0.5\theta = 0.5? ¿Y 0.90.9? (b) La mioglobina une O2_2 en un solo sitio: la saturación sigue a Langmuir en la presión parcial de oxígeno — ¿por qué se observa una curva hiperbólica (y no sigmoidal) para la mioglobina y una curva cooperativa en S para la hemoglobina, de cuatro sitios? (c) El monóxido de carbono se une a los mismos sitios con un ϵ0\epsilon_0 mayor en 0.13eV\approx0.13\,\mathrm{eV}: a presiones parciales iguales y 310K310\,\mathrm{K}, calcúlese la razón de ocupaciones CO:O2_2 y explíquese el envenenamiento por CO. (d) ¿Por qué funciona la oxigenoterapia hiperbárica, en lenguaje de isotermas?

Solución

Solución de Ejercicio 18.6.

(a) P=P0P = P_0; P=9P0P = 9P_0. (b) Un solo sitio independiente da la hipérbola P/(P+P0)P/(P + P_0) — la mioglobina; los cuatro sitios que interactúan de la hemoglobina se cargan de forma cooperativa y empinan la curva hasta la forma en S que se carga por completo en los pulmones y se descarga en los tejidos. (c) e0.13/0.0267130\eu^{0.13/0.0267} \approx 130: a presiones parciales iguales, el CO gana los sitios por un factor cien — una traza de CO desaloja al oxígeno de la sangre. (d) Subir la presión parcial de O2_2 sube su μ\mu y vuelve a inclinar la competencia: la isoterma se disputa en μ\mu, y el oxígeno hiperbárico simplemente puja más alto.

Ejercicio 18.7 ★★

Fluctuaciones del número. (a) Muéstrese que VarN=kBTN/μ\operatorname{Var}N = k_{\text{B}}T\,\partial\langle N\rangle/\partial\mu. (b) Para el gas ideal clásico, dedúzcase VarN=N\operatorname{Var}N = \langle N\rangle: Poisson. (c) Reconcíliese con las fluctuaciones de densidad del Ejercicio 16.12. (d) Para un solo nivel fermiónico, calcúlese Varn=nˉ(1nˉ)\operatorname{Var}n = \bar n(1 - \bar n): ¿dónde están más silenciosos los números de fermiones, y por qué?

Solución

Solución de Ejercicio 18.7.

(a) Derívese N=kBTμlnΞ\langle N\rangle = k_{\text{B}}T\,\partial_\mu \ln\Xi una vez más. (b) Neβμ\langle N\rangle \propto \eu^{\beta\mu}: μN=βN\partial_\mu\langle N\rangle = \beta\langle N\rangle, luego VarN=N\operatorname{Var}N = \langle N\rangle. (c) La gaussiana de varianza nˉ\bar n de aquel ejercicio es este enunciado de Poisson, visto en un subvolumen. (d) Varn=nˉ(1nˉ)\operatorname{Var}n = \bar n(1 - \bar n): nulo cuando el nivel está con certeza vacío o con certeza lleno — el bloqueo de Pauli acalla las fluctuaciones en el fondo de un mar de Fermi.

Ejercicio 18.8 ★★

La ósmosis a partir de μ\mu. Una membrana deja pasar el agua pero no el soluto (de concentración cc por volumen, diluido). (a) Igualando el μ\mu del agua a los dos lados de la membrana, muéstrese que el lado de agua pura debe estar sobrepresurizado en Π=ckBT\Pi = ck_{\text{B}}T (van ’t Hoff) — acéptese que el soluto disuelto baja el μ\mu del agua en kBTc/nw-k_{\text{B}}T\,c/n_{\text{w}} por molécula de agua. (b) Calcúlese Π\Pi para una solución salina fisiológica (c3×1026m3c \approx 3 \times 10^{26}\,\mathrm{m}^{-3} de partículas totales). (c) Glóbulos rojos en agua pura: ¿en qué sentido fluye el agua y con qué resultado? (d) ¿Por qué no necesitan los árboles ni bombas ni milagros para subir la savia decenas de metros (gradientes osmóticos y capilares de μ\mu)?

Solución

Solución de Ejercicio 18.8.

(a) El soluto diluye la entropía del agua y baja su μ\mu; la presión sube μ\mu: el balance fija Π=ckBT\Pi = ck_{\text{B}}T. (b) Π=3×1026×1.38×1023×3101.3MPa\Pi = 3 \times 10^{26} \times 1.38 \times 10^{-23} \times 310 \approx 1.3\,\mathrm{MPa} — trece atmósferas dentro de las venas, compensadas a diario. (c) El μ\mu del agua pura supera al del interior celular: el agua entra en tromba y las células revientan (hemólisis) — por eso las infusiones son salinas y nunca de agua pura. (d) Los gradientes de μ\mu de la raíz a la hoja (solutos, evaporación en la hoja) tiran del agua hacia arriba como diferencias de presión osmótica y capilar: decenas de metros sin una sola pieza móvil.

Ejercicio 18.9 ★★

Equilibrio en un campo. En un campo gravitatorio, el equilibrio de una columna de gas exige que el potencial químico total μ(h)=kBTln(nλT3)+mgh\mu(h) = k_{\text{B}}T\ln(n\lambda_T^3) + mgh sea uniforme. (a) Dedúzcase de aquí la ley barométrica en una línea. (b) Generalícese a una centrifugadora de velocidad angular ω\omega (con potencial 12mω2r2-\tfrac12 m\omega^2r^2) y obténgase el perfil radial. (c) Las centrifugadoras de enriquecimiento de uranio giran a ωr600m/s\omega r \sim 600\,\mathrm{m}/\mathrm{s}: calcúlese la razón de equilibrio del enriquecimiento de 238^{238}UF6_6 frente a 235^{235}UF6_6 por etapa, exp[Δmω2r2/2kBT]\exp[\Delta m\,\omega^2 r^2/2k_{\text{B}}T] con Δm=5×1027kg\Delta m = 5 \times 10^{-27}\,\mathrm{kg} y T=320KT = 320\,\mathrm{K}. (d) ¿Por qué se encadenan muchas etapas?

Solución

Solución de Ejercicio 18.9.

(a) μ\mu uniforme: kBTlnn(h)+mgh=k_{\text{B}}T\ln n(h) + mgh = const., es decir, la exponencial barométrica. (b) n(r)e+mω2r2/2kBTn(r) \propto \eu^{+m\omega^2r^2/2k_{\text{B}}T}: las especies pesadas se agolpan en el borde. (c) exp[Δmω2r2/2kBT]=exp[5×1027×3.6×105/(2×4.4×1021)]e0.21.2\exp[\Delta m\,\omega^2r^2/2k_{\text{B}}T] = \exp[5 \times 10^{-27} \times 3.6 \times 10^{5}/(2 \times 4.4 \times 10^{-21})] \approx \eu^{0.2} \approx 1.2 por etapa. (d) Pasar del 0.7%0.7\% natural al grado de reactor o al militar exige muchos factores 1.21.2 compuestos: de ahí las cascadas de miles de centrifugadoras — y por eso contar cascadas es diplomacia nuclear.

Ejercicio 18.10 ★★★

Saha, a mano. (a) Realícese la demostración con el balance de μ\mu de la Proposición 18.5. (b) Defínase la fracción ionizada xx a densidad total de hidrógeno nn y muéstrese que x2/(1x)=eEI/kBT/nλT,e3x^2/(1 - x) = \eu^{-E_{\text{I}}/k_{\text{B}}T}/n\lambda_{T,\text{e}}^3. (c) Fotosfera solar: T=5800KT = 5800\,\mathrm{K}, n1023m3n \approx 10^{23}\,\mathrm{m}^{-3}: muéstrese que la fracción ionizada de hidrógeno es solo de 104\sim10^{-4} — la superficie del Sol es neutra. (d) A T=10000KT = 10\,000\,\mathrm{K} (con la misma nn): recalcúlese y coméntese lo abrupto del cambio.

Solución

Solución de Ejercicio 18.10.

(a) Como en el texto: exponénciese el balance de μ\mu con potenciales químicos de gas ideal y el desplazamiento EI-E_{\text{I}} del átomo. (b) Con ne=np=xnn_{\text{e}} = n_{\text{p}} = xn y nH=(1x)nn_{\text{H}} = (1-x)n: el cociente enunciado. (c) A 5800K5800\,\mathrm{K}: λT,e1.0nm\lambda_{T,\text{e}} \approx 1.0\,\mathrm{nm}, e27.2=1.5×1012\eu^{-27.2} = 1.5 \times 10^{-12}, 1/nλ31.1×1041/n\lambda^3 \approx 1.1 \times 10^{4}: x2/(1x)1.6×108x^2/(1-x) \approx 1.6 \times 10^{-8} y x1.3×104x \approx 1.3 \times 10^{-4} — la fotosfera solar es hidrógeno neutro al 99.99%99.99\%. (d) A 10000K10\,000\,\mathrm{K} (con la misma nn): x6%x \approx 6\% — dos órdenes y medio en un factor 1.71.7 de temperatura; a las densidades electrónicas menores de la fotosfera la subida llega aún antes. Las transiciones de Saha son acantilados.

Ejercicio 18.11 ★★★

Cuando los gases se vuelven cuánticos. La degeneración se instala en el punto nλT31n\lambda_T^3 \sim 1. (a) Para el aire a 300K300\,\mathrm{K}: ¿cuánto por debajo del umbral (Ejercicio 18.2)? (b) Para los electrones de conducción del cobre (n=8.5×1028m3n = 8.5 \times 10^{28}\,\mathrm{m}^{-3}): hállese la temperatura por debajo de la cual nλT3>1n\lambda_T^3 > 1 — y conclúyase que los electrones de un metal están degenerados a cualquier temperatura de laboratorio. (c) Para átomos de rubidio a n=1020m3n = 10^{20}\,\mathrm{m}^{-3}: la temperatura umbral cuántica (la escala del condensado de Bose–Einstein del Capítulo 19). (d) Ordénense los tres y enúnciese la regla: las partículas ligeras y las densidades altas se vuelven cuánticas primero.

Solución

Solución de Ejercicio 18.11.

(a) 1.7×1071.7 \times 10^{-7}: siete órdenes por debajo del umbral — el aire es clásico sin discusión. (b) Hágase nλT3=1n\lambda_T^3 = 1: Tdeg=h2n2/3/2πmkB6×104KT_{\text{deg}} = h^2n^{2/3}/2\pi mk_{\text{B}} \approx 6 \times 10^{4}\,\mathrm{K} para los electrones del cobre — todos los metales que se han tocado alguna vez llevan un gas electrónico profundamente cuántico. (c) Para el rubidio a 1020m310^{20}\,\mathrm{m}^{-3}: Tdeg107KT_{\text{deg}} \sim 10^{-7}\,\mathrm{K} — la frontera de los cien nanokelvin del Capítulo 19. (d) Electrones en metales: siempre cuánticos; átomos de laboratorio: solo a nanokelvin; aire: nunca — la degeneración va como n2/3/mn^{2/3}/m.

Ejercicio 18.12 ★★★

El universo pierde sus positrones. En el universo primitivo caliente, e++e2γ\eu^+ + \eu^- \rightleftharpoons 2\gamma mantenía los pares en equilibrio con μe++μe=0\mu_{e^+} + \mu_{e^-} = 0 (los fotones: μ=0\mu = 0). (a) Para un plasma simétrico esto da μe±0\mu_{e^\pm} \approx 0: muéstrese que la densidad de pares es entonces n±λT3emec2/kBTn_\pm \approx \lambda_T^{-3}\eu^{-m_{\text{e}}c^2/k_{\text{B}}T} una vez que kBTmec2k_{\text{B}}T \ll m_{\text{e}}c^2. (b) Calcúlense la temperatura a la que mec2=kBTm_{\text{e}}c^2 = k_{\text{B}}T y la escala de tiempo cósmico que le corresponde (acéptese t1st \sim 1\,\mathrm{s} a 1010K10^{10}\,\mathrm{K}). (c) Muéstrese que un enfriamiento adicional modesto (un factor de unas pocas unidades) desploma la densidad de pares en muchos órdenes de magnitud: los positrones se aniquilan y dejan solo el diminuto exceso primordial de electrones. (d) Los fotones liberados recalientan el baño de radiación: ¿qué radiación reliquia lleva el recibo (Capítulo 27)?

Solución

Solución de Ejercicio 18.12.

(a) Con μ=0\mu = 0, la densidad de pares es la λT3eβmec2\lambda_T^{-3}\eu^{-\beta m_{\text{e}}c^2} suprimida por Boltzmann, una vez que los pares dejan de ser relativistas. (b) mec2/kB=5.9×109Km_{\text{e}}c^2/k_{\text{B}} = 5.9 \times 10^{9}\,\mathrm{K}, alcanzada alrededor de un segundo tras el comienzo. (c) Enfriar de 6×109K6 \times 10^{9}\,\mathrm{K} a 2×109K2 \times 10^{9}\,\mathrm{K} multiplica el exponente por tres: la densidad se desploma como e2\eu^{-2} \to e6\eu^{-6} solo en la exponencial, y en órdenes más con el prefactor — en cuestión de segundos los positrones desaparecen, aniquilados contra todos los electrones salvo una parte en mil millones. (d) Los fotones de aniquilación se sumaron al baño de radiación y lo recalentaron un poco; su resto estirado es el fondo cósmico de microondas — el recibo a 2.7K2.7\,\mathrm{K} del suceso.

18.6 Problema: leer las temperaturas de las estrellas

Problema 18.1

Problema de fin de semana — Saha, Payne y de qué está hecho el universo

Los espectros estelares se catalogaron a lo largo de la década de 1890 en las clases O, B, A, F, G, K y M según la intensidad de sus líneas de absorción — con las líneas de Balmer del hidrógeno máximas en la clase A y misteriosamente débiles tanto en las estrellas más calientes como en las más frías. En 1925, Cecilia Payne aplicó la flamante ecuación de Saha y demostró que todo el zoo era una misma sustancia a distintas temperaturas — y que las estrellas son abrumadoramente hidrógeno, conclusión tan herética que tuvo que suavizarla en su tesis. Este problema recorre su razonamiento. Datos: la absorción de Balmer parte de n=2n = 2, a E2E1=10.2eVE_2 - E_1 = 10.2\,\mathrm{eV} sobre el estado fundamental; EI=13.6eVE_{\text{I}} = 13.6\,\mathrm{eV}; densidades electrónicas fotosféricas ne1020m3n_{\text{e}} \sim 10^{20}\,\mathrm{m}^{-3}.

Parte I — Dos factores en competencia.

  1. La absorción de Balmer necesita hidrógeno neutro en el nivel n=2n = 2: escríbase el factor de Boltzmann para la población de n=2n = 2 (degeneraciones: g2/g1=4g_2/g_1 = 4).
  2. Evalúese a 50005000, 1000010000 y 20000K20\,000\,\mathrm{K}: ¿en qué sentido empuja este factor la intensidad de la línea al subir TT?
  3. Ahora el factor de Saha: al subir TT, ¿qué le ocurre a la fracción neutra y, por tanto, a las líneas de Balmer?
  4. Explíquese cualitativamente por qué el producto — excitado y todavía neutro — debe tener un máximo a alguna temperatura intermedia.
  5. Usando la ecuación de Saha con ne=1020m3n_{\text{e}} = 10^{20}\,\mathrm{m}^{-3}, estímese la temperatura de semiionización del hidrógeno.
  6. Júntelo todo en el máximo: ¿cerca de qué temperatura deberían ser más intensas las líneas de Balmer, y a qué clase espectral (las estrellas A, 10000K\sim10\,000\,\mathrm{K}) señala eso?

Parte II — La secuencia espectral, descifrada.

  1. Las estrellas M (3000K3000\,\mathrm{K}) muestran bandas moleculares intensas (TiO) e hidrógeno débil: explíquense ambas cosas con los dos factores.
  2. Las estrellas O (40000K40\,000\,\mathrm{K}) muestran líneas de helio y de nuevo hidrógeno débil: explíquese (qué le ha pasado al hidrógeno y por qué le llega ahora el turno al helio, con EI=24.6eVE_{\text{I}} = 24.6\,\mathrm{eV}).
  3. El Sol (G, 5800K5800\,\mathrm{K}) muestra líneas intensas de calcio ionizado (EI,Ca=6.1eVE_{\text{I,Ca}} = 6.1\,\mathrm{eV}) y, sin embargo, hidrógeno débil: muéstrese con la lógica de Saha por qué el calcio ya está ionizado mientras que el hidrógeno sigue neutro y sobre todo en n=1n = 1.
  4. Los astrónomos anteriores a Payne leían la intensidad de línea como abundancia: ¿qué les hicieron concluir sobre la composición del Sol sus poderosas líneas de calcio y sus enclenques líneas de Balmer?
  5. La intuición de Payne: corregir la intensidad de cada línea con su “factor de visibilidad” de Boltzmann y Saha. Cuando lo hizo, ¿qué jerarquía de abundancias verdaderas emergió?
  6. Su resultado — hidrógeno un millón de veces más abundante que el calcio, y estrellas hechas sobre todo de hidrógeno y helio — ¿por qué se le opuso resistencia (qué sugería la química centrada en la Tierra) y cuándo se la reivindicó?

Parte III — Números sobre el máximo.

  1. Calcúlese el factor de Boltzmann de n=2n = 2 exactamente a 9600K9600\,\mathrm{K}.
  2. Con la estimación de Saha del apartado 5, tómese la fracción neutra a 96009600, 1500015000 y 25000K25\,000\,\mathrm{K} como aproximadamente 0.50.5, 10310^{-3} y 10510^{-5}: fórmese el producto (fracción excitada ×\times fracción neutra) a esas tres temperaturas y localícese su máximo.
  3. ¿Por qué es abrupto el máximo de Balmer por el lado caliente (qué exponencial gana) y gradual por el lado frío?
  4. La fotosfera de una enana blanca tiene una nen_{\text{e}} mil veces mayor: ¿en qué sentido se desplaza la temperatura de semiionización, y por qué (Le Chatelier en lenguaje de μ\mu)?
  5. Dos estrellas de idéntica temperatura pero distinta presión muestran, así, intensidades de línea distintas: ¿qué propiedad estelar (gigante frente a enana — la gravedad superficial) permite eso leer a los espectroscopistas?
  6. Los cartografiados modernos fijan las temperaturas estelares con un ±1%\pm1\% a partir de razones de líneas: enúnciese la cadena postulado \to Boltzmann \to Saha \to termómetro.

Parte IV — Qué quedó zanjado.

  1. La secuencia OBAFGKM había sido un caos alfabético reordenado empíricamente: tras Saha, ¿qué única variable física la ordena?
  2. Dígase por qué el mismo elemento puede dominar un espectro a una temperatura y desvanecerse a otra sin ningún cambio de abundancia — la lección central que Payne enseñó a la astronomía.
  3. El hidrógeno, con tres cuartas partes de toda la masa bariónica: nómbrense otros dos pilares de este libro que descansan en esa abundancia (la fusión estelar, Capítulo 27; el cielo de 21 cm, Capítulo 12).
  4. Otto Struve llamó a la tesis de Payne “la tesis doctoral más brillante jamás escrita en astronomía”: a partir de este problema, justifíquese el elogio en una frase.
  5. La misma física de Saha, aplicada a la inversa a escala cósmica, predice la temperatura a la que se combinó por primera vez todo el hidrógeno del universo (3000K\sim3000\,\mathrm{K}, a la densidad cósmica mucho menor): nómbrense el suceso y su reliquia (Capítulo 27).
  6. Los equilibrios de ionización gobiernan también las lámparas fluorescentes y el procesado por plasma: dígase en una frase por qué puede existir un estado ionizado “de tipo 10000K10\,000\,\mathrm{K}” en un tubo cuyo vidrio permanece frío (¿qué dos subsistemas no llegan a equilibrarse?).
  7. Resúmase el resultado con nombre propio: una sola condición de equilibrio, μH=μp+μe\mu_{\text{H}} = \mu_{\text{p}} + \mu_{\text{e}}, convierte el alfabeto espectral en un termómetro que va de los 30003000 a los 40000K40\,000\,\mathrm{K}, sitúa el máximo de Balmer en las estrellas A y revela un universo hecho de hidrógeno.
Solución

Solución de Problema 18.1.

1. n2/n1=4e10.2eV/kBTn_2/n_1 = 4\,\eu^{-10.2\,\text{eV}/k_{\text{B}}T}. 2. 2×10102 \times 10^{-10}, 2.9×1052.9 \times 10^{-5} y 1.1×1021.1 \times 10^{-2}: suben abruptamente — cuanto más caliente, mejor para excitar n=2n = 2. 3. Saha ioniza el hidrógeno y lo hace desaparecer: por encima de 10000K\sim10\,000\,\mathrm{K} la fracción neutra se desploma y se lleva consigo las líneas de Balmer. 4. Un factor sube con TT y el otro baja: su producto tiene que alcanzar un máximo entre ambos regímenes. 5. Imponer λT,e3eEI/kBT=ne=1020m3\lambda_{T,\text{e}}^{-3}\eu^{-E_{\text{I}}/ k_{\text{B}}T} = n_{\text{e}} = 10^{20}\,\mathrm{m}^{-3} da T9600KT \approx 9600\,\mathrm{K}. 6. Cerca de 10000K10\,000\,\mathrm{K} — exactamente las estrellas A, cuyos espectros son el escaparate del hidrógeno. 7. A 3000K3000\,\mathrm{K} nada llega a n=2n = 2 (e39\eu^{-39}), mientras que moléculas frágiles como el TiO sobreviven y absorben en bandas: las estrellas frías son moleculares y oscuras en hidrógeno. 8. El hidrógeno está ionizado hasta la invisibilidad; el helio, con su escalera de 24.6eV24.6\,\mathrm{eV}, solo ahora dispone del kBTk_{\text{B}}T y de las condiciones de Saha para poblar y retener sus estados absorbentes. 9. La ionización del calcio, de 6.1eV6.1\,\mathrm{eV}, es mucho más barata que la del hidrógeno, de 13.6eV13.6\,\mathrm{eV}: a 5800K5800\,\mathrm{K}, Saha ioniza el calcio casi por completo (las líneas del Ca+^+ absorben entonces desde su estado fundamental — sin impuesto de Boltzmann), mientras que el hidrógeno neutro se queda en n=1n = 1, ciego a Balmer. 10. Que el Sol es sobre todo calcio y hierro, parecido a la corteza terrestre — la hipótesis reinante que Payne heredó. 11. Al dividir cada intensidad por su factor de visibilidad, las diferencias descomunales se vinieron abajo: casi todos los elementos salieron en proporciones parecidas — salvo el hidrógeno (y el helio), que resultó ser un millón de veces más abundante. 12. Se daba por hecho que la química cósmica reflejaba la roca terrestre; Russell la presionó para que calificara el resultado del hidrógeno de “casi con certeza irreal” — y publicó él mismo la misma conclusión, con reconocimiento, cuatro años después. Payne tenía razón. 13. 4e10.2/0.827=4e12.31.8×1054\eu^{-10.2/0.827} = 4\eu^{-12.3} \approx 1.8 \times 10^{-5}. 14. Productos 9×106\approx 9 \times 10^{-6}, 1.5×1061.5 \times 10^{-6} y 3.5×1073.5 \times 10^{-7}: el máximo cae en 9600K\approx9600\,\mathrm{K} — el pico de las estrellas A, reproducido con dos exponenciales. 15. Por el lado caliente, la exponencial de Saha aniquila de golpe la población neutra; por el lado frío, el factor de Boltzmann apenas mengua: acantilado arriba, ladera abajo. 16. Una nen_{\text{e}} mayor empuja la recombinación (hay más parejas que capturar): la semiionización se desplaza hacia lo caliente — se aprieta la reacción hacia el lado ligado, Le Chatelier en lenguaje de μ\mu. 17. La gravedad superficial: las enanas compactas (presión alta) y las gigantes hinchadas (baja) se separan a una misma temperatura — clases de luminosidad leídas de la física de las líneas. 18. Estados equiprobables a priori \to poblaciones de Boltzmann \to ionización de Saha \to razones de intensidad de línea invertidas para dar TT: un termómetro sin piezas móviles a cualquier distancia. 19. La temperatura — un solo eje, de los 3000K3000\,\mathrm{K} de las M a los 40000K40\,000\,\mathrm{K} de las O, sustituye a todo el alfabeto. 20. La visibilidad son poblaciones, y las poblaciones son exponenciales en TT: la intensidad de una línea espectral mide primero la termodinámica de la estrella y solo después, tras corregirla, su química. 21. La fusión del hidrógeno como fuente de energía estelar y la cartografía de 21 cm de las galaxias: ambas se apoyan en la abundancia de Payne. 22. Con dos fórmulas de colectivos convirtió cuarenta años de taxonomía espectral en termometría y volvió a pesar la composición del universo — en una tesis doctoral. 23. La recombinación cósmica: la semiionización del propio universo a 3000K\sim3000\,\mathrm{K}, cuya luz liberada, estirada mil veces, es el fondo de microondas a 2.7K2.7\,\mathrm{K}. 24. Los electrones (calientes, de escala 104K10^{4}\,\mathrm{K}) y el gas más las paredes (frías) intercambian energía demasiado despacio para equilibrarse: dos temperaturas coexisten en un mismo tubo — un estado estacionario fuera del equilibrio, y por eso la lámpara brilla sin fundirse. 25. Una sola condición, μH=μp+μe\mu_{\text{H}} = \mu_{\text{p}} + \mu_{\text{e}}, ordena OBAFGKM solo por temperatura, aparca el máximo de Balmer en las estrellas A cerca de 9600K9600\,\mathrm{K} y — una vez corregidas las intensidades por visibilidad — revela un universo de hidrógeno.

Términos definidos en este capítulo

Ver los 431 términos del glosario