Showing posts with label Física Estadística. Show all posts
Showing posts with label Física Estadística. Show all posts

Bosonización: dualidad fermiones/bosones

En diciembre Joseph Polchinski publicó un artículo recopilatorio sobre dualidades llamado simplemente Dualities (al que luego agregó 'of Fields and Strings'):

Dualities of Fields and Strings
Joseph Polchinski

Duality, the equivalence between seemingly distinct quantum systems, is a curious property that has been known for at least three quarters of a century. In the past two decades it has played a central role in mapping out the structure of theoretical physics. I discuss the unexpected connections that have been revealed among quantum field theories and string theories.

El artículo está muy bien escrito, de manera sencilla y sucinta, y creo que es accesible, al menos a nivel conceptual, incluso a estudiantes de licenciatura. Luego de introducir el tema de dualidades, Polchinski bosqueja algunas dualidades QFT/QFT, entre las que se encuentra la de bosonización (antes, otra cuestión que llama la atención es la de la dualidad onda/partícula como dos límites clásicos de QFT), que es básicamente una dualidad entre una teoría de fermiones de Dirac y una teoría de bosones en 1+1 dimensiones. Ésta dualidad es sorprendente y me pareció sumamente interesante por lo sencilla y porque con mi poco conocimiento de QFT supuse que podría trabajar los detalles; además descubrí que es prácticamente tema estándar en el área de Materia Condensada y seguramente tiene bastantes avances o generalizaciones.

Polchinski se limita a presentar la dualidad: las relaciones entre las constantes de acoplamiento y la correspondencia entre los campos que se deben satisfacer para que las amplitudes sean idénticas. Para obtener el resultado, original de Sidney Coleman (1975), uno tiene varias opciones, como intentar trabajar como sugiere Polchinski: empezar con teorías sin masa y sin interacción, $\mc{L}_\text{F}=i\bar\psi\gamma^\mu\partial_\mu\psi$ y $\mc{L}_\text{B}=-\frac{1}{2}\partial_\mu\phi\partial^\mu\phi$, calcular las amplitudes de dispersión, obtener la correspondencia que las vuelve idénticas y luego deformar nuevamente las teorías a unas con interacción. Y esto es básicamente lo que hace Coleman en su artículo original, calculando funciones de Green en ambos casos:

Quantum sine-Gordon equation as the massive Thirring model
Sidney Coleman

The sine-Gordon equation is the theory of a massless scalar field in one space and one time dimension with interaction density proportional to $\cos\beta\phi$, where $\beta$ is a real parameter. I show that if $\beta^2$ exceeds $8\pi$, the energy density of the theory is unbounded below; if $\beta^2$ equals $4\pi$, the theory is equivalent to the zero-charge sector of the theory of a free massive Fermi field; for other values of $\beta$, the theory is equivalent to the zero-charge sector of the massive Thirring model. The sine-Gordon soliton is identified with the fundamental fermion of the Thirring model.

El artículo es relativamente extenso y confieso que no lo he seguido paso a paso, lo que creo que quizá convendría únicamente si uno quiere entrarle en serio al tema de bosonización.

De cualquier modo, buscando un poco más sobre bosonización, además de encontrar varios títulos de libros sobre el tema enfocados a Física de Partículas y a Materia Condensada, fortuitamente también encontré otro artículo sobre el tema:

Bosonization as Duality
C.P. Burgess, F. Quevedo

We show that bosonization in two dimensions can be derived as a special case of the duality transformations that have recently been used to good effect in string theory. This allows the construction of the bosonic counterpart of any fermionic theory simply by `following your nose' using the standard duality transformation rules. We work through the bosonization of the Dirac fermion, the massive and massless Thirring models, and a fermion on a cylindrical spacetime as illustrative examples.

El artículo comienza con el comentario

Precise rules for various versions of this equivalence have been known for some time now [1][2][3].The statement of these rules generally suffers from the drawback of not being constructive in character, being instead couched in terms of the equivalence of a given pair of bose and fermi theories¹. This is particularly true for nonabelian bosonization, where a simple, systematic, constructive and unified presentation of the bosonization rules is still missing [4]. Our purpose in this note is to provide such a derivation, for the abelian case.

donde abeliano se refiere a que el grupo de simetría de la teoría en cuestión es conmutativo; los mismos autores tratan en otro artículo el caso no-abeliano (también tratado años antes por Witten). Este artículo me pareció bastante mono (al menos de un vistazo) por la sencillez y generalidad. Los autores mencionan un procedimiento llamado dualización que para entonces ya había sido útil para relacionar 'vacíos' aparentemente no relacionados en teoría de cuerdas y lo describen así:

This duality is a trick for constructing the equivalent theory by embedding the model of interest within a larger gauge theory, which reduces to the original version once the gauge potential, $A_\mu$, is set to zero. This larger gauge theory is typically constructed simply by gauging a global symmetry of the original model [7] [8]. The original model is then written as the gauged theory subject to some constraint that removes the gauge potential, and this is usually done by introducing a dummy variable –a Lagrange multiplier field, $\Lambda$– whose functional integral requires the gauge field strength to vanish, $F_{\mu\nu}=0$. (...) If $\Lambda$, and then $A_\mu$, are integrated out –in this order– of the gauged version of the theory, then the original model is retrieved. The dual version of the theory, on the other hand, is obtained by performing the various functional integrals in a different order, in particular by integrating out the original field, as well as $A_\mu$. The Lagrange multiplier, $\Lambda$, then plays the role of the new, dual, field.

La idea en realidad es sencilla, y como mencionan, puede entenderse mejor sólo con el caso concreto del fermión no masivo (por alguna razón lo mencionan como libre) en la sección 2.1). El procedimiento es básicamente plantear la función de partición de los fermiones
\begin{equation}Z_\text{F}[J]=\int\mc{D}\psi\,\exp\left[i\int{d^2x}\left(\bar\psi\gamma^\mu\partial_\mu\psi+\sum_iJ_i\mc{O}_i\right)\right]\end{equation} con los $J_i(x)$ y $\mc{O}_i(\psi)$ campos y operadores, respectivamente (que en el caso más sencillo, como hacen en el artículo, pueden tomarse como constantes de acoplamiento y términos de interacción) y encajarla en una teoría de norma más amplia (básicamente 'normando' una simetría global de la teoría original) sujeta a ciertas constricciones
\begin{equation}Z_\text{gauge}[J]=\int\mc{D}\psi\,\mc{D}A_\mu\,\mc{D}\Lambda\,\exp\left[i\int{d^2x}\left(\mc{L}_\text{F}(\psi,J)+\mc{L}_\text{gauge}(A_\mu,\Lambda)\right)\right]\end{equation} que permita recuperar $Z_\text{F}[J]$ al integrar en $\Lambda$ y en $A_\mu$, i.e. tal que $Z_\text{F}[J]={Z}_\text{gauge}[J]\bigg|_{\Lambda,A_\mu}$. Finalmente esto permite obtener la teoría dual en $\Lambda$ (que como se muestra en el artículo, corresponde a una de bosones) si se integra $Z_\text{gauge}[J]$ en $\psi$ y en $A_\mu$,
\begin{equation}Z_\text{B}[J]\equiv{Z}_\text{gauge}[J]\bigg|_{\psi,A_\mu}=\int\mc{D}\Lambda\exp\left[i\int{d^2x}\,\mc{L}_\text{B}(\Lambda,J)\right]\end{equation} con lo que se obtiene también la correspondencia entre las teorías identificando los términos en $J_i$.

La idea en realidad es sencilla, la parte que resulta complicada para un neófito como yo es la de elegir la teoría de norma y la de hacer la integración, ya sea con fuerza bruta o devorándose los artículos que mencionan los autores. De cualquier modo la parte importante es si quiera saber en primer lugar que esta sorprendente dualidad existe, siendo que el carácter de ambas teorías es tan distinta, y luego entender cómo funciona y el alcance que tiene.

Finalmente, regresando al artículo de Polchinski, en realidad no era consciente de la idea general que plantea acerca de una dualidad como un par de límites de una teoría más fundamental, así que definitivamente me ha causado una impresión importante su artículo (o la parte que he leído seriamente hasta ahora) y en general lo recomiendo ampliamente.

If the doors of perception were cleansed every thing would appear to man as it is, Infinite.

― William Blake, The Marriage of Heaven and Hell


Agujeros Negros en AdS

He estado ausente algún tiempo intentando terminar la licenciatura, y como mostré en la entrada anterior, he estado trabajando en la geometría (clásica) del espacio de Anti-de Sitter. Ahora quiero compartir algunos comentarios sobre Agujeros Negros en AdS, basándome en este ensayo de Peng Zhao.

El caso más sencillo es (A)dS-Schwarzschild, y la idea es que partiendo de una métrica esféricamente simétrica en forma de Schwarzschild,
\begin{equation}ds^2=-f(r)\,dt^2+f^{-1}(r)\,dr^2+r^2d\Omega^2\end{equation} uno puede tomar $f$ tal que cerca del origen se tenga un espacio de Schwarzschild y asintóticamente se tenga (A)dS (para una obtención más metódica de esta métrica véase e.g. la p.823 de las notas de M. Blau). Las métricas de AdS y de Schwarzschild, respectivamente, pueden escribirse como
\begin{align}ds_\text{AdS}^2&=-\left(1+\frac{r^2}{\alpha^2}\right)\,dt^2+\left(1+\frac{r^2}{\alpha^2}\right)^{-1}dr^2+r^2d\Omega^2\\
ds_\text{Sch}^2&=-\left(1-\frac{r_s}{r}\right)\,dt^2+\left(1-\frac{r_s}{r}\right)^{-1}dr^2+r^2d\Omega^2\end{align} donde $\alpha$ es una constante positiva (el radio de AdS) y $r_s$ es el radio de Schwarzschild. De aquí, se define entonces la métrica de AdS-Schwarzschild como
\begin{equation}ds^2=-\left(1+\frac{r^2}{\alpha^2}-\frac{r_s}{r}\right)\,dt^2+\left(1+\frac{r^2}{\alpha^2}-\frac{r_s}{r}\right)^{-1}dr^2+r^2d\Omega^2\end{equation} que en efecto se comporta como un agujero negro de Schwarzschild para $r$ pequeña y como AdS para $r$ grande. Lo propio puede hacerse empleando de Sitter.

***(Actualización)
La primera cuestión es a AdS se le puede asociar una energía finita y por tanto el análogo a una pared de potencial infinita en su infinito asintótico. Por teorema de Noether, en AdS la energía es una carga conservada, $E=-P_\mu{V}^\mu$ (el signo negativo es simplemente por consistencia, de cualquier modo siempre se le puede cargar a $E$ que es constante); en este caso nos interesa la energía $E_\text{loc}$ medida por un observador local (es decir, que localmente está en un sistema plano) estático, en donde tendríamos simplemente que $P=(E,\vec{0})$ y $V=(1,\vec{0})$. Lo que nos interesa es encontrar esta relación para este observador en el sistema global de coordenadas de AdS (esféricas). Considérense los vectores transformados $p$ y $v$. Primero empleemos la condición de normalización $v^\mu{v}_\mu=-1$ y el hecho de que el observador es estático también en el sistema global de AdS, i.e. $v=(v^t,\vec{0})$, esto implica que $v^t=1/\sqrt{-g_{tt}}$, entonces en este sistema coordenado,
\begin{equation}E_\text{loc}=-p_\mu{v}^\mu=-g_{\mu\nu}p^\mu{v}^\nu=\sqrt{-g_{tt}}p^t\end{equation} Ahora bien, en el sistema del observador local estático simplemente $p^t=m\dot{t}$ con $m$ la masa de la partícula (el observador), que podemos llevar a la forma en coordenadas de AdS sabiendo que $P^t=g_{tt}\dot{t}=-E$ en AdS, de modo que $p^t=-mE/g_{tt}\equiv-{E}_0/g_{tt}$ y así,
\begin{equation}E_\text{loc}(r)=\frac{E_0}{\sqrt{-g_{tt}(r)}}\label{1}\end{equation} de donde se lee que la energía medida por un observador local en algún punto $r$ tiene un corrimiento al rojo (debido al campo gravitacional) respecto a la energía $E_0$. En el ensayo de Zhao se escribe $E_0$ como $E_\infty$, sin embargo la notación se antoja un poco desafortunada, ya que ésta es la energía total medida por un observador local estático en $r=0$; la notación probablemente proviene de la interpretación de ser la energía medida por un observador alejado infinitamente del campo gravitacional del espacio, lo que resulta útil e.g. en Schwarzschild, donde $E$ en efecto es la energía (por unidad de masa) medida por un observador local estático en $r\to\infty$, aquí de hecho $E_\text{loc}\to0$ cuando $r\to\infty$, contrario a Schwarzschild, en donde $E_\text{loc}\to{mE}$ cuando $r\to\infty$ y $E_\text{loc}\to\infty$ cuando $r\to{r}_s$, esto hace la gran diferencia para ver a AdS como una caja. Nótese también que los vectores de momento y de velocidad (i.e. el vector de Killing transformado) se recuperan para el caso local en AdS con $r=0$ y en Schwarzschild con $r\to\infty$. Como se ve adelante, pueden construirse estados térmicos en AdS identificando periódicamente una coordenada de tiempo imaginaria con periodo $\beta_0=T^{-1}$. Esto significa entonces que (\ref{1}) puede escribirse también como $T_\text{loc}=\frac{T}{\sqrt{-g_{tt}}}$ y que la energía total de la radiación térmica en AdS es finita sin necesidad de encerrarla en una caja.
***

Luego entonces, tomando $r=r_+$ el horizonte de eventos de AdS-Schwarzschild, i.e. la raíz más grande de
\begin{equation}V(r)\equiv1-\frac{r_s}{r}+\frac{r^2}{\alpha^2}\end{equation} se fija el valor de $r_s$ a
\begin{equation}V(r_+)=0\,\Longrightarrow\,r_s=\frac{r_+(\alpha^2+r_+^2)}{\alpha^2}\end{equation} y entonces en efecto al ver la métrica cerca del horizonte, i.e. con $r=r_++\epsilon^2$ con $|\epsilon|\ll1$, para el orden dominante en $\epsilon$,
\begin{equation}V(r_++\epsilon^2)\approx\frac{\alpha^2+3r_+^2}{\alpha^2r_+}\epsilon^2\end{equation} y haciendo una rotación de Wick, $\tau=it$, se sigue para la métrica AdS-Schwarzschild,
\begin{equation}ds^2=\frac{4\alpha^2r_+}{\alpha^2+3r_+^2}\left[d\epsilon^2+\left(\frac{\alpha^2+3r_+^2}{2\alpha^2r_+}\right)^2\epsilon^2d\tau^2\right]+(r_++\epsilon^2)^2d\Omega^2\end{equation} Aquí puede resultar poco claro cómo ver que necesariamente $\tau$ es periódica. El punto es entender lo que significa que se tenga una singularidad cónica. La idea es básicamente que si $\tau$ tiene un cierto periodo $\beta_0$, el espacio se cubrirá completamente, lo que no sucederá necesariamente con cualquier otro valor; la analogía que utiliza Zhao es la de las coordenadas polares $ds^2=dr^2+r^2d\theta^2$, donde $\theta$ debe tener un periodo de $2\pi$ para corresponder a la métrica del plano, de otro modo se tiene un hueco del tipo de una rebanada de pastel, lo que lleva a algo así como un pico de un cono, lo que es la singularidad cónica.
En este caso es particularmente fácil verificar cuál es el periodo $\beta_0$ de $\tau$, dado que el término entre corchetes de la métrica AdS-Schwarzschild cerca del horizonte es análoga a la del plano en coordenadas polares. Se debe satisfacer entonces que
\begin{equation}\frac{\alpha^2+3r_+^2}{2\alpha^2r_+}\int_0^{\beta_0}d\tau\int_0^{r_+}\epsilon\,{d}\epsilon=\pi{r}_+^2\,\Longrightarrow\,\beta_0=\frac{4\pi\alpha^2r_+}{\alpha^2+3r_+^2}\end{equation}
En adelante la cosa es bastante clara, al menos antes de llegar a la parte de AdS/CFT, y el resultado $T=\beta^{-1}$ creo que es bastante majo, pues entonces la temperatura del agujero negro puede obtenerse únicamente con argumentos geométricos, sin invocar directamente la gravedad de superficie. Sólo hay detalles cuando se menciona la distribución canónica, aunque al final igual lo que se calcula es $S=\beta\langle{E}\rangle+\ln{Z}$. La intención es llegar a la dualidad fluido/gravedad, pero de eso tal vez escriba después  ;-)

Utilizando mecánica clásica en átomos: partículas idénticas y caos cuántico

Específicamente quiero responder por qué es cierto que cuando la longitud de onda de deBroglie en un gas a una temperatura dada es mucho más pequeña que la distancia entre partículas, el gas se comporta clásicamente (cada partícula o átomo puede modelarse con la mecánica de Newton o bien, el gas sigue una distribución de Maxwell-Boltzmann). En general el tema del límite clásico de la mecánica cuántica es uno mucho más amplio, pues ambas teorías son radicalmente distintas.

En este caso, fijémonos en una partícula; por el teorema de Ehrenfest, se tiene que
$$\frac{d\langle{p}\rangle}{dt}=-\left\langle{\frac{dV(x)}{dx}}\right\rangle$$ desarrollando en serie de potencias el lado derecho de la ecuación alrededor de ${\langle{x}\rangle}$,
$$\frac{dV(x)}{dx}=\frac{dV(\langle{x}\rangle)}{d\langle{x}\rangle}+\frac{dV(\langle{x}\rangle)^2}{d\langle{x}\rangle^2}(x-\langle{x}\rangle)+\frac{1}{2}\frac{dV(\langle{x}\rangle)^3}{d\langle{x}\rangle^3}(x-\langle{x}\rangle)^2+\mathcal{O}\left(\langle{x}\rangle^4\right)$$ ahora bien, ${\left\langle{x-\langle{x}\rangle}\right\rangle=0}$, y ${\left\langle(x-\langle{x}\rangle)^2\right\rangle=\sigma_x^2}$, de modo que si V varía lentamente en x, pueden considerarse únicamente los primeros términos de la expansión,
$$\frac{d\langle{p}\rangle}{dt}=-\frac{dV(\langle{x}\rangle)}{d\langle{x}\rangle}-\frac{1}{2}\sigma_x^2\frac{dV(\langle{x}\rangle)^3}{d\langle{x}\rangle^3}$$ a continuación al despreciar la varianza ${\sigma_x}$, lo que se está haciendo es precisamente decir que el tamaño de la función de onda es muy pequeña comparada con la variación del potencial V, pues la varianza se interpreta como qué tan extendida o dispersa es la distribución, en este caso espacial, de la partícula. Esto entonces da como resultado
$$\frac{d\langle{p}\rangle}{dt}=-\frac{dV(\langle{x}\rangle)}{d\langle{x}\rangle}$$ precisamente la ecuación para la segunda ley de Newton, lo que es de esperarse, pues significa que la probabilidad de hallar a la partícula en x estará abismalmente concentrada en su valor medio.

Esto entonces se puede interpretar como que en una descripción clásica de cada partícula (átomos si se quiere), la extensión de la función de onda (prácticamente la longitud de onda de deBroglie), es mucho más pequeña que la distancia entre partículas, que son las fuentes de potencial V. Esto es lo que suele decirse en textos básicos de física estadística o de dinámica molecular.

Imagen (si ves este texto, recarga la página)Hablando en general, en nuestro universo todas las partículas conocidas entran únicamente en dos tipos: bosones y femiones, que se conocen como partículas idénticas. Las partículas idénticas en el mundo cuántico son además indistinguibles. Es interesante entonces que en este límite, si se quiere, en el que partículas idénticas se comportan como partículas clásicas, éstas dejarán de ser indistinguibles. Lo dicho antes entonces puede traducirse en qué tanto interfieren las funciones de onda de las partículas, como en el caso del gas ideal en condiciones estándar, en que las moléculas están suficientemente separadas de modo que pueden considerarse distinguibles, en cambio dos electrones en el mismo nivel energético de un átomo siempre serán indistinguibles. La imagen mostrada a la izquierda pertenece a skeptics play, en donde puede leerse más a nivel divulgación sobre partículas idénticas.

Finalmente, una parte relacionada con este tema y por la que he estado muy interesado últimamente (por supuesto a nivel fanático estudiante de licenciatura en física) es la del caos cuántico, que básicamente busca la relación del caos (que únicamente se concibe como algo clásico) con la mecánica cuántica, es decir, cuando se observa caos clásicamente, ¿qué ocurre en la descripción cuántica?. Lo que he leído hasta ahora ha sido sobre todo acerca de la llamada dinámica de billar, que como el nombre sugiere, básicamente describe el movimiento de una partícula encerrada por una frontera rígida, en la cual se refleja.

Imagen (si ves este texto, recarga la página)
Estadio de Bunimovich: descripción clásica
Imagen (si ves este texto, recarga la página)
Estadio de Bunimovich: descripción cuántica

Imagen (si ves este texto, recarga la página)Las imágenes muestran el llamado estadio de Bunimovich, un estadio caótico, del cual puedes leer más en esta entrada de Terence Tao. Acá también pongo una animación del estadio de Bunimovich clásico. A la derecha se muestra un llamado corral cuántico. Esta última imagen la puedes encontrar en la página del American Institute of Physics, así como leer más detalle en este artículo. En general creo que esta es un área de investigación muy activa recientemente y hay resultados interesantísimos por parte de muchos investigadores. Si ya eres estudiante de alguna carrera relacionada a la física, o tienes inclinación hacia ella, aquí hay una gran motivación para que te decidas de lleno. Por parte de México también me he encontrado con trabajos interesantísimos, como éste, en que se generalizan los estadios dos dimensionales clásicos. Sea del lado clásico o del cuántico, el estudio de sistemas dinámicos y por supuesto el caos, me parece sumamente atractivo.

Fluctuaciones de energía en el ensamble canónico

La consistencia entre los formalismos microcanónico y canónico en el límite termodinámico, suele ser demostrada luego de estudiar la dispersión de los valores de energía de un sistema en el ensamble canónico (el ensamble microcanónico considera sistemas aislados, en los que la energía está restringida a un intervalo muy pequeño). Esto es, luego de analizar las fluctuaciones de energía en el ensamble canónico.

Consideremos que el hamiltoniano es simplemente la energía, ${\mathcal{H}=E}$. La forma común de obtener una expresión sencilla para las fluctuaciones de energía, es notar que éstas están dadas por
$$\left\langle\left(\Delta{E}\right)^2\right\rangle\equiv\left\langle{E^2}\right\rangle-\left\langle{E}\right\rangle^2=-\frac{\partial^2\ln{z}}{\partial\beta^2}=-\frac{\partial{U}}{\partial\beta}$$ donde z es la función de partición y ${\beta\equiv(\mathrm{k}T)^{-1}}$ en este caso, además se ha empleado la definición
$$U\equiv\langle{E}\rangle=-\frac{\partial\ln{z}}{\partial\beta}$$ (demuestra la última igualdad). La razón de que esto sea cierto es sencilla, ya que con la distribución canónica en el espacio fase $\Gamma$,
$$\langle{\mathcal{H}}\rangle=\frac{\displaystyle{\int{\mathcal{H}\,\mathrm{e}^{-\beta{\mathcal{H}}}}\,\mathrm{d}\Gamma}}{\displaystyle{\int{\mathrm{e}^{-\beta{\mathcal{H}}}}\,\mathrm{d}\Gamma}}$$ (calcula la derivada parcial y compruébalo). Hace poco, un problema de examen se me ha complicado bastante por no notar lo que he dicho antes en tan pocas líneas. Estás ya cansado de todo el periodo curricular; aún te falta esta prueba y te piden demostrar que
$$\left\langle{(\Delta{E})^3}\right\rangle=\mathrm{k}^2\left\{T^4\left(\frac{\partial{C_v}}{\partial{T}}\right)_V+2T^3C_v\right\}$$ Bueno pues lo que hay que hacer es prácticamente lo anterior, primero identifíquese
\begin{align*}\left\langle{E^2}\right\rangle-U^2&=\left\langle{E^2}\right\rangle-2U^2+U^2\\[0.1in]&=\left\langle{E^2-2EU+U^2}\right\rangle\\[0.1in]&=\left\langle{(E-U)^2}\right\rangle\\[0.1in]&\equiv\left\langle\left(\Delta{E}\right)^2\right\rangle\end{align*} lo que significa entonces que
$$\left\langle{(\Delta{E})^3}\right\rangle=\left\langle{E^3}\right\rangle-3\left\langle{E^2}\right\rangle{U}+2U^3$$ y curiosamente te tropiezas con que
\begin{align*}\frac{\partial^2U}{\partial\beta^2}&=\frac{\int{E}^3\mathrm{e}^{-\beta{E}}\,\mathrm{d}\Gamma}{\int{\mathrm{e}}^{-\beta{E}}\,\mathrm{d}\Gamma}-3\frac{\int{E}^2\mathrm{e}^{-\beta{E}}\,\mathrm{d}\Gamma}{\int\mathrm{e}^{-\beta{E}}\,\mathrm{d}\Gamma}\cdot\frac{\int{E}\mathrm{e}^{-\beta{E}}\,\mathrm{d}\Gamma}{\int\mathrm{e}^{-\beta{E}}\,\mathrm{d}\Gamma}+2\left[\frac{\int{E}\mathrm{e}^{-\beta{E}}\,\mathrm{d}\Gamma}{\int\mathrm{e}^{-\beta{E}}\,\mathrm{d}\Gamma}\right]^3\\[0.1in]&=\left\langle{E^3}\right\rangle-3\left\langle{E^2}\right\rangle\,U+2U^3\end{align*} (compruébalo) lo que lleva a completar la demostración haciendo la termodinámica correspondiente.

Es interesante esta relación, que probablemente pueda generalizarse más, aunque lo desconozco. Nota también la elección del signo para la definición de ${\Delta{E}}$; si has seguido hasta acá la entrada, da un argumento que justifique tal elección.

Energía sobre la pared

kinetic theory of gasesUn resultado impresionante, y con el que seguido se introduce en general a la física estadística, es la conexión entre la dinámica de un sistema microscópico con la termostática del sistema macroscópico vía la teoría cinética de Maxwell. Además de hacer las hipótesis de que las moléculas que conforman un gas son esferas duras, que colisionan elásticamente, que las condiciones iniciales del sistema (posición y velocidad) están distribuidas al azar (ésto es precisamente lo que vuelve probabilística la naturaleza del problema), que el gas está en equilibrio termodinámico y que las velocidades de las partículas son estadísticamente independientes de la misma distribución de probabilidades, Maxwell considera la energía para hacer la conexión con la termodinámica. Una forma de hacer la conexión es directamente encontrar la distribución que describen las partículas y de ahí determinar el valor esperado de la energía $\displaystyle{\frac{1}{2}m\left\langle{v^2}\right\rangle}$ de cada partícula, para finalmente hacer la conexión con la energía interna de un gas ideal. Sin embargo otra forma más interesante y que reafirma que la descripción de Maxwell va por buen rumbo es hacer la conexión por medio de la presión.

Supongamos que hicimos la tarea y encontramos que el número de moléculas que golpean la unidad de área de una superficie plana por segundo, con velocidades entre v y ${v+dv}$ y ángulos entre $\theta$ y ${\theta+d\theta}$ es
$$\frac{1}{2}vnf(v)\,dv\sin\theta\cos\theta\;d\theta$$ donde
$$f(v)\propto{v^2e^{-\frac{m}{2kT}v^2}}$$ es la distribución de Maxwell. Cada molécula que golpea la pared tiene un cambio en el momento de ${2mv\cos\theta}$, entonces para la presión p,
\begin{align*}dp&=\left(\frac{1}{2}vnf(v)\,dv\sin\theta\cos\theta{d\theta}\right)(2mv\cos\theta)\\[0.1in]&=mnv^2f(v)\,dv\cos^2\theta\sin\theta{d\theta}\end{align*} es decir
$$p=mn\int_0^\infty{v^2f(v)\,dv}\int_0^{\pi/2}\cos^2\theta\sin\theta\,d\theta=\frac{1}{3}mn\langle{v^2}\rangle$$ que, comparando con la ley de gas ideal, ${p=nkT}$ se obtiene
$$\frac{1}{2}m\langle{v^2}\rangle=\frac{3}{2}kT$$ que es un resultado bastante impresionante, pues aunque desde niños se nos dice que la temperatura es en realidad una manifestación de qué tan rápido se mueven dadas partículas, aquí lo estamos probando matemáticamente.

Una forma más directa, pero más tortuosa de obtener el mismo resultado es simplemente calcular el cociente,
$$\frac{1}{2}m\langle{v^2}\rangle=\frac{1}{2}m\frac{\displaystyle{\int_0^\infty{v^4e^{-\frac{m}{2kT}v^2}}dv}}{\displaystyle{\int_0^\infty{v^2e^{-\frac{m}{2kT}v^2}}dv}}$$ y el resultado es exactamente el mismo.

Ahora bien, al preguntarnos qué pasa únicamente con aquellas 'esferitas' que golpean la superficie, el problema se vuelve exactamente el mismo que el de la llamada efusión. La efusión "selecciona" moléculas más rápidas (por ley de Graham, más ligeras), pues éstas tienen una mayor probabilidad de alcanzar el orificio. La misma situación aplica si queremos determinar la energía promedio de las moléculas que de hecho golpean la superficie. Esto significa que la distribución de las moléculas que golpean la pared, ${\tilde{f}(v)}$, justo como con la efusión, tomarán un factor v extra, esto es.
$$\tilde{f}(v)\propto{v^3}e^{-\frac{m}{2kT}v^2}$$ Ahora bien, vemos que podemos escribir la siguiente integral usando la función Gamma,
$$\int_0^\infty{v^n}e^{-\alpha{v^2}}dv=\frac{1}{2}\alpha^{-\frac{n+1}{2}}\left(\frac{n-1}{2}\right)!$$ de este modo siempre podemos calcular fácilmente esta integral, sobre todo para n impares, y así, calculando directamente la energía cinética para esta distribución,
\begin{align*}\frac{1}{2}m\langle{v^2}\rangle&=\frac{1}{2}m\frac{\displaystyle{\int_0^\infty{v^5e^{-\frac{m}{2kT}v^2}}dv}}{\displaystyle{\int_0^\infty{v^3e^{-\frac{m}{2kT}v^2}}dv}}\\[0.1in]&=\frac{1}{2}m\left(\frac{\displaystyle{\frac{8k^3T^3}{m^3}}}{\displaystyle{\frac{2k^2T^2}{m^2}}}\right)\\[0.1in]&=2kT\end{align*} que es un resultado desconcertante sin el argumento de la efusión, sin embargo se vuelve coherente entendiendo lo anterior.

El origen de la vida, complejidad en el universo y la segunda ley de la termodinámica

Me he matriculado en un curso básico de física biológica y ha resultado una experiencia sumamente enriquecedora; incluso reconociéndome aficionado al uso de herramientas matemáticas hasta en la sopa y hasta cierto grado de formalismo con éstas, debo reconocer que la biología me resulta tan excitante como la física y las matemáticas, y considero su estudio tan indispensable para cualquier humano como el estudio de las últimas. Una discusión en clase se ha tornado al origen de la vida teniendo en consideración principalmente las leyes de la física. Seguro alguna vez, así te dediques o no a la ciencia -simplemente por ser humano-, te has preguntado cómo es que la maquinaria de lo que llamamos vivo es tan perfecta, y espero te lo sigas preguntando (aún más, que te lo sigas preguntando después de leer esto).

En termodinámica clásica generalmente se consideran tres leyes (o se puede decir que cuatro: las leyes cero, uno, dos y suele llamarse tercera ley al postulado de Nernst) empíricas. La segunda ley en particular es clave en la discusión del origen de la vida. Imagina que estás en una solitaria playa y decides hacer una escultura de arena, ¿qué pasará con ésta una vez que la dejes ahí sola? Bueno, todos sabemos por experiencia que los granitos de arena poco a poco se separarán hasta que quede sólo un bulto feo de arena de lo que alguna vez fue tu escultura. ¿Por qué sucede esto?, ¿por qué no pueden los granos de arena preservar o formar esculturas por sí solos? Pues estas preguntas no son triviales y hay muchos expertos en física trabajando en el área llamada Cosmología para resolver esta cuestión. Una descripción que, digamos funge como respuesta satisfactoria, se llama flecha del tiempo, y está relacionada íntimamente con el concepto de entropía.

La entropía suele describirse popularmente como desorden, sin embargo a los físicos les gusta describirla como configuraciones posibles de un sistema. Regresando al ejemplo de la escultura arenosa: tu escultura tiene muy pocas configuraciones posibles, i.e. baja entropía, ya que, como seguramente eres un artista muy exigente, hay quizás una o dos o tres formas precisas de acomodar cada granito de arena de modo que obtengas la escultura que deseas. Por otro lado, hay muchas configuraciones posibles para un bulto feo de arena, es decir, el bulto posee una entropía mucho mayor que la de tu escultura. La segunda ley de la termodinámica describe el comportamiento de la entropía (de un sistema cerrado), diciendo que ésta jamás disminuirá espontáneamente respecto a cualquier otra variable termodinámica.

Imagen (si ves este texto, recarga la página)Ludwig Boltzmann hizo la formulación de la segunda ley para el enfoque estadístico de la termodinámica y encontró que de hecho es posible que la entropía disminuya, simplemente es muy poco probable. Bueno, pues aquí empieza lo divertido de esta publicación. Nota primero que para la flecha del tiempo, lo que llamamos antes y después corresponde precisamente a baja y alta entropía respectivamente, ¿por qué entonces, en primer lugar, el universo comenzaría en un estado de baja entropía? ¡Además, pareciera que los seres vivos contradecimos la segunda ley de la termodinámica clásica y somos esa minúscula probabilidad de la termodinámica estadística! Intenta llenar una cubeta con carbono, calcio, hierro y demás elementos que nos forman... y jamás saldrá un humano de la cubeta. Se podría decir incluso que somos agentes anti-entropía, desde las proteínas y el DNA, hasta todo lo que concebimos como nosotros. Es simplemente impresionante la perfección de las máquinas vivas y deben haber condiciones muy específicas que las generen, no sólo elementos químicos. En este momento de la discusión suelen atribuirse inmediatamente poderes sobrenaturales, dioses, etc...

Imagen (si ves este texto, recarga la página)A estas cuestiones se debe decir primero que los seres vivos y nuestro entorno constituimos lo que se llama un sistema termodinámico abierto, es decir, podemos intercambiar masa y energía con los alrededores. Imagina un recipiente cerrado totalmente lleno de agua, para el cual pones por abajo una estufa y por arriba un bloque de hielo. Lo que sucederá es que (ley cero) el agua buscará el equilibrio, haciendo circular el calor recibido a través del recipiente (transferencia de calor por convección). Básicamente ésto es lo que ocurre cuando la Tierra recibe energía de baja entropía (i.e. energía a alta temperatura) desde el Sol.

Lo siguiente es un extracto del libro Theoretical physics and philosophical problems de Ludwig Boltzmann del capítulo Second Law of Thermodynamics escrito (o pronunciado) en 1886 (desconozco si existe alguna edición en español).

Imagen (si ves este texto, recarga la página)Between the earth and sun there is a colossal temperature difference; between these two bodies energy is thus not at all distributed according to the laws of probability. The equalization of temperature, based on the tendency towards greater probability, takes millions of years, because the bodies are so large and are so far apart. The intermediate forms assumed by solar energy, until it falls to terrestrial temperatures, can be fairly improbable, so that we can easily use the transition of heat from sun to earth for the performance of work, like the transition of water from the boiler to the cooling instillation. The general struggle for existence of animate beings is not a struggle for raw materials – these, for organisms, are air, water and soil, all abundantly available – nor for energy, which exists in plenty in any body in the form of heat (albeit unfortunately not transformable), but of a struggle for entropy, which becomes available through the transition of energy from the hot Sun to the cold Earth. In order to exploit this transition as much as possible, plants spread their immense surface of leaves and force the Sun's energy, before it falls to the earth's temperature, to perform in ways yet unexplored certain chemical syntheses of which no one in our laboratories has so far the least idea. The products of this chemical kitchen constitute the object of struggle of the animal world.

De cualquier modo esto da pie a algunos problemas Cosmológicos como la formación de estructura que dio pie a la formación del Sol en primer lugar y el hecho de que, respetando la 2a ley, el universo en general haya partido de una baja entropía en primer lugar. La proposición también atribuida principalmente a Boltzmann (relacionada con el llamado cerebro antrópico de Boltzmann) a partir de su enfoque estadístico, se resume a que en un universo que parte del equilibrio, la entropía puede fluctuar, disminuyendo de vez en cuando y aumentando de nuevo, es decir, quizá nuestro universo está en medio de una fluctuación de entropía, lo que a su vez permite en un sentido más profundo el desarrollo y evolución de la vida. Esto, como se ha dicho, es probable, aunque lo sea poco.

Ahora bien, aunque ésta es una idea genial y consistente que se deshace de dioses, poderes sobrenaturales o universos finamente ajustados, sigue siendo una propuesta incompleta; es casi como decir que las cosas son como son porque sí. Lo más importante es que la propuesta se verifique y pueda hacer predicciones sobre el comportamiento y naturaleza del universo. Más detalles sobre esto pueden leerse en esta entrada del blog de Sean Carroll sobre el cerebro antrópico de Boltzmann.