Consultas extra

Hoy viernes, de 16 a 18, pueden pasar por la oficina de Florencia (2144). Si no, pueden recurrir a la ayuda de nuestros representantes virtuales en el Campus. Importantes novedades, aquí mismo, más tarde.

Clase práctica del lunes pasado [actualizado]

Finalmente, las notas de la clase práctica del lunes quedaron pasadas en limpio. Pueden bajarlas [aquí] (las actualizamos con un breve comentario que se publicó en el Campus). En el límite de infinitas iteraciones, la red de este ejemplo es un fractal de dimensión igual a 2. Miren la figura del post anterior. La red completa puede subdividirse en cuatro subredes autosemejantes a la red original. Si se duplica el tamaño de estas subredes, se obtiene la red original. Es lo mismo que pasa con un cuadrado bidimensional. Si se divide al cuadrado en cuatro cuadrados de lado mitad, el cuadrado original puede obtenerse duplicando el tamaño de cualquiera de los cuadrados pequeños. En general, si un conjunto puede dividirse en N partes iguales autosemejantes al conjunto original, y el factor de escala es S, la dimensión del conjunto es tal que

El triángulo de Sierpinski es un ejemplo de fractal de dimensión no entera. Apliquen la fórmula anterior y calculen su dimensión fractal.

Últimas clases

La clase práctica de ayer, resolvimos un problema del libro de Pathria y Beale. Es el último problema del capítulo 14, pero no por ser el último es el más difícil. En realidad, no tiene ninguna cuenta complicada.

No hay (todavía) un apunte de clase, pero pueden ver [aquí] la presentación en pdf que usé ayer. Con esto tenemos una nueva familia de problemas para el parcial; ya veníamos repitiendo mucho.

Hemos empezado a corregir sus TP de Ising. Desde el momento en el que reciban el trabajo observado, van a tener una semana para enviar las correcciones.

Mañana hay clase de consulta desde las 17.

Clases prácticas de la semana

Mientras seguimos a la espera de sus prácticas computacionales, les recuerdo que las dos clases prácticas que tuvimos esta semana fueron muy importantes. El lunes vimos el grupo de renormalización aplicado a la cadena lineal. [Aquí] pueden bajar las notas de clase del cuatrimestre pasado. Encuentran este problema en la mayoría de los libros; miren, por ejemplo, los libros de Pathria y Beale, de Huang, de Dalvit et al.

El miércoles vimos la aproximación de Migdal-Kadanoff para la transformación del grupo de renormalización de la red cuadrada. [Aquí] pueden bajar la presentación de Bernardo, preparada para ver en pantalla completa; [aquí] pueden bajar las notas del cuatrimestre pasado. Es fácil inventar problemas parecidos con otras redes.

El lunes hay clase de práctica normal. El miércoles habrá clase de consultas desde las 17. El parcial es el lunes 8, con el mismo protocolo que el primer parcial. Estén atentos al viernes que viene, porque vamos a dar algunas palabras clave sobre los problemas. Previsiblemente, habrá un problema de fermiones, otro de bosones y uno que combine renormalización con Ising. En 2019, tomamos dos problemas de esta clase (parcial, recuperatorio). Con eso agotamos nuestra imaginación. El cuatrimestre pasado hubo un problema con una cadena decorada. Si quieren un problema nuevo de este tipo, tomen el problema de la escalera de Ising e inserten espines en el medio de todos los enlaces. La suma sobre esos espines puede hacerse explícitamente; encuentren la constante de acoplamiento renormalizada.

La semilla del mal

La clase práctica del miércoles pasado, vimos la relación entre la energía libre de Landau y la energía libre de Gibbs en la aproximación de Bragg-Williams para el modelo de Ising. Es decir, vimos un ejemplo concreto en donde la energía libre de Landau aparece naturalmente, y no como un artefacto a medida para explicar de manera fenomenológica las transiciones de fase. De nuevo vimos la importancia de tomar el límite termodinámico, sin el cual no habría ninguna singularidad a temperatura finita. En relación con esto, la animación de abajo muestra el valor medio del espín en función del campo externo para una temperatura mayor a la crítica, pero para un número finito de espines. Sólo tiene sentido hablar de magnetización espontánea cuando N tiende a infinito. 

[Aquí] pueden bajar las notas de clase, que corresponden al estreno mundial del problema 1 de la Guía 8, que sienta un precedente de ítem r.  [Aquí] pueden bajar la Guía 8.

Supongo que están muy concentrados en el trabajo computacional, así que no los molesto más.

Para los que entendieron

Es un gran hallazgo de la mecánica estadística el descubrimiento de que la transformación de cirujano en hechicero supremo es una transición de fase crítica de segundo orden.

La semana próxima empezamos con grupo de renormalización. Veremos los temas que doy usualmente en la materia, y algunos temas nuevos. Con instrucciones para que alguno de ustedes pueda convertirse en hechicera o hechicero. Y como estamos llegando a hacer la última guía de trabajos prácticos con tiempo, tengan en cuenta que grupo de renormalización va a entrar en el parcial.

Abre tu ojo


Propuse varias soluciones; todas, insuficientes. Las discutimos; al fin, Stephen Albert me dijo:
- En una adivinanza cuyo tema es el ajedrez ¿cuál es la única palabra prohibida? Reflexioné un momento y repuse:
- La palabra ajedrez.
- Precisamente —dijo Albert—. El jardín de los senderos que se bifurcan es una enorme adivinanza, o parábola, cuyo tema es el tiempo; esa causa recóndita le prohíbe la mención de su nombre.

Jorge Luis Borges, El jardín de los senderos que se bifurcan (1941).

Si este post fuera una adivinanza, no podríamos mencionar la palabra “autosemejanza”. Vamos a hablar de fractales y de senderos que se bifurcan. Los fractales son objetos matemáticos que, por construcción, son invariantes de escala (es decir, se prescriben con un conjunto de reglas, usualmente recursivas, que generan una figura o un conjunto autosemejante). Es importante notar que el fenómeno de autosemejanza que se observa en los sistemas físicos cerca del punto crítico no se genera de esta forma, con pasos que se repiten infinitamente. Y en este sentido, los fractales no nos pueden brindar una explicación a la causa de la invariancia de escala. Sin embargo, como objetos matemáticos, pueden servir para estudiar propiedades generales de sistemas que son invariantes de escala, para generar datos sintéticos que tengan esta propiedad (como para generar terrenos o texturas que parezcan realistas en videojuegos), o para crear métodos para cuantificar la posible autosemejanza de un conjunto de datos.

Uno de los ejemplos más sencillos y conocidos está dado por el conjunto de Cantor. Se construye tomando el segmento [0,1], partiéndolo en 3, y removiendo el segmento del medio. Esto nos deja dos nuevos segmentos: [0,1/3] y [2/3,1]. La operación se repite en cada uno de los nuevos segmentos. La figura a continuación muestra el resultado de repetir este procedimiento diez veces (hagan click en la imagen para ver un zoom):

En términos coloquiales, un fractal es una figura construida con pequeñas partes que son similares a la figura completa, en cualquier escala en la que se observe. La construcción recursiva del fractal (que puede ser determinista, o tener componentes aleatorias) asegura que la figura resultante sea autosemejante. Y su “fractalidad” puede cuantificarse de diferentes formas; por ejemplo, calculando funciones de correlación y sus exponentes críticos. O calculando la dimensión fractal o la dimensión de Hausdorff, que están relacionadas con el exponente crítico de la función de correlación a dos puntos.

El término “fractal” fue introducido por Benoit Mandelbrot, que formalizó varias ideas previas de otros matemáticos (especialmente, durante el siglo XX, de Lewis Fry Richardson, que también hizo contribuciones importantes a la meteorología y a la turbulencia). Y fueron usados por Mandelbrot para, entre otras aplicaciones, calcular el perímetro de regiones costeras. Aunque la costa irregular de un país no es generada por una persona que repite reglas como en el conjunto de Cantor (pero en The hitchhiker’s guide to the galaxy pueden opinar distinto), calcular la dimensión fractal de la costa permite obtener buenas estimaciones de la longitud de curvas muy rugosas, y en ciento sentido, autosemejantes. Los que estén interesados en los detalles pueden leer el paper (lindo, clásico, y muy breve) de Mandelbrot sobre este tema:

De la misma forma que conocer la longitud de correlación en el modelo de Ising nos permite inferir propiedades del tamaño de los dominios magnéticos, estimar la dimensión fractal le permitió a Mandelbrot resolver una aparente paradoja al intentar calcular la longitud de curvas autosemejantes: al medir la longitud de una costa, cuanto más detalle se tiene sobre su forma, más aumenta su longitud.

Los fractales también pueden generar imágenes visualmente interesantes, como el famoso conjunto de Mandelbrot:

Los que estén interesados en generar fractales con Python pueden ver los siguientes links con instrucciones paso a paso (recomiendo fuertemente el primero), y muchos ejemplos de códigos que pueden cortar y pegar en sus computadoras o en un Google Colab:

Como mencioné más arriba, los fractales pueden tener componentes aleatorias. Y aunque los fractales no brindan una explicación a la causa de la autosemejanza en ciertos sistemas naturales, pueden ser usados para caracterizarla. Además, cumplen teoremas muy interesantes que nos permiten descubrir relaciones sorprendentes entre procesos autosemejantes. Por ejemplo, los ceros de un camino al azar unidimensional de tiempo contínuo (es decir, cada vez que el caminante al azar vuelve a pasar por su punto de origen) forman un conjunto fractal. Esto tiene que ver con otro teorema muy extraño que se aplica a un proceso llamado evolución de Schramm-Loewner: una curva al azar en dos dimensiones que sea invariante conforme (una forma más fuerte de la invariancia de escala, en la que la curva no es solo invariante frente a cambios de escala, sino también invariante frente a transformaciones que preserven los ángulos localmente) tiene una relación directa con un proceso de movimiento browniano en una dimensión. Este teorema puede usarse para calcular exponentes críticos en modelos de Ising y de percolación en dos dimensiones, a partir de propiedades del movimiento browniano unidimensional que vimos al principio del curso. ¡Todo se conecta con todo! De pronto, un tema de esta materia viajó al pasado y tuvo un hijo.

Grupos para la práctica computacional

Les acabo de enviar un mail a todos con un link a una planilla para que anoten los grupos de la práctica computacional. Traten de resolver eso cuanto antes, así empiezan a trabajar. Hay una columna para los que están sin grupo. Si se organizan entre ustedes, mejor. Si no, veo cómo distribuir a los que quedaron sin grupo.

¡Se viene el Criticalpalooza!


Llegó el festival anual psicodélico de la criticalidad en la física y en áreas afines, cada año actualizado con papers más y más recientes. Tenemos entradas para ver ejemplos de fenómenos críticos en física de materiales, sistemas biológicos, fluidos, cosmología, y hasta en neurociencias. Los links de bandas que tocaron en un festival de nombre similar llevan a videos en YouTube. Y ahora sí, criticalidad y autosemejanza:

Comencemos con el lineup. En el escenario principal, después de Guns N’ Roses, tenemos dos ejemplos de física de materiales y materia condensada. Ya conocemos las transiciones de fase asociadas a los cambios de estado de la materia (sólido, líquido y gaseoso). ¡Pero existen muchos más cambios de fase! Así que veamos el show de bandas menos conocidas (y con ejemplos recientes de publicaciones en física). Primero tenemos a los cristales líquidos. Los cristales líquidos están formados por moleculas anisótropas (generalmente alargadas, por ejemplo con forma de largos cilindros), por lo que se comportan con propiedades intermedias entre los líquidos y los sólidos de acuerdo a en que dirección del material se aplican los esfuerzos. Y tienen al menos dos fases: en la fase nemática (a mayor temperatura) las moléculas están más desordenadas, pero se alinean a lo largo de sus ejes principales. En cambio, en la fase esméctica (a menor temperatura) las moléculas se acomodan en capas más ordenadas, y dentro de cada capa las moléculas están inclinadas con el mismo ángulo. El cambio entre ambas fases es una transición de fase con propiedades críticas. Los interesados pueden mirar un paper reciente (Gim, Beller & Yoon, Nat. Commun. 8, 15453, 2017), donde encontrarán esta figura alucinante con cambios morfológicos del material durante la transición (para temperatura creciente, de izquierda a derecha):

También pueden mirar este paper (el preprint de acceso libre está disponible acá) en el que los autores demuestran la existencia de una fase nemática en un modelo de cristales líquidos formados por barras en una red bidimensional. Y los que quieran leer otro ejemplo de transiciones de fase en materiales, pueden mirar este artículo sobre transiciones sólido-sólido que ocurren por cambios en la forma de partículas coloidales.

Vayamos a otro escenario del Criticalpalooza, y mientras escuchamos The Flaming Lips, veamos ejemplos recientes de transiciones de fase y autosemejanza observados en sistemas biológicos. Ciertas células cambian sus patrones de movimiento de acuerdo a la densidad de células en su entorno. A baja densidad muestran un movimiento desordenado, mientras que a alta densidad muestran patrones de movimiento colectivo y ordenado. La transición entre ambos comportamientos ocurre como una transición de fase. Pueden ver un ejemplo en Szabó et al., Phys. Rev. E 74, 061908 (2006) (el preprint está diponible en este link). La siguiente figura, de ese paper, muestra los patrones de movilidad al aumental la densidad de las células (de izquierda a derecha). Noten el cambio en el orden del sistema:

Los que tengan interés por ver más aplicaciones en biofísica, pueden mirar también un paper sobre fenómenos críticos en membranas lípidas (¡donde estiman exponentes críticos!).

En el mismo escenario donde toca Metallica tenemos una variedad de fenómenos autosemejantes que se observan en fluidos. El más conocido es el fenómeno de la turbulencia, que es heavy metal. Los fluidos más viscosos (donde la viscosidad se mide con un número adimensional, el número de Reynolds) fluyen en forma ordenada y laminar. Pero al aumentar el número de Reynolds (y reducirse la importancia de la viscosidad), generan flujos muy desordenados que tienen propiedades de invariancia de escala. En la siguiente figura, noten como zooms sucesivos en el flujo parecen repetir los patrones (algo similar a lo que vimos en el modelo de Ising en 2D):

Esta es una transición extraña, porque los flujos turbulentos tienen propiedades de criticalidad y autosemejanza, pero son un sistema fuera del equilibrio. Pueden ver más imágenes de flujos turbulentos aquí. Y pueden ver un ejemplo de una transición de fase en turbulencia en un condensado de Bose-Einstein en este paper.

Los que tengan interés por cosmología pueden ir al escenario donde tocan Empire of the Sun y Lee Smolin, y leer este artículo donde se discuten diversos problemas (como la formación de galaxias espirales, la estructura de gran escala del universo, o el universo temprano) desde el punto de vista de fenómenos críticos. Bastante más difícil de leer, pero que sirve como ilustración, es este paper donde se estudia una transición de fase al compactificar una dimensión en teorías de gravedad (compactificar es “enrollar” una dimensión sobre si misma, y hacer tender ese “rollo” a cero). Y en cosmología, algunas teorías predicen que otra transición de fase podría haber ocurrido en el universo temprano, cuando se formaron los primeros átomos. Luego de ese momento el universo se volvió transparente, pero antes de ese instante el universo puede haber sido opalescente en forma crítica (como la mezcla de un líquido y un gas cuando llega a la temperatura crítica y dejan de existir las diferencias entre ambas fases, como se ve en la foto del centro en esta imágen de una mezcla de etano líquido y gaseoso tomada de Wikipedia):

Ya que estamos en el festival, no nos olvidemos de ir a ver las bandas clásicas. El mecanismo de Higgs por el cual los bosones de gauge adquieren masa, es también un mecanismo de ruptura espontánea de la simetría.

Para cerrar el festival, este Criticalpalooza no tendrá a The Strokes, pero tiene ejemplos de fenómenos críticos en el cerebro (¿lo entendieron?). Las funciones cognitivas involucran procesos que van desde las neuronas individuales hasta regiones grandes del cerebro, y en los últimos años se encontró evidencia creciente de que el cerebro funciona en el borde entre el orden y el desorden, con propiedades de fenómenos críticos. Los interesados en ver ejemplos de criticalidad y autosemejanza en el cerebro pueden mirar este paper o este paper. O, si quieren mirar dos papers más recientes, pueden ver este paper publicado en Science Advances o este paper escrito por varios investigadores del Departamento de Física. Ambos tienen gráficos que muestran cómo se calculan los exponentes críticos en este sistema. La siguiente imagen, tomada de Cochi et al., Progress in Neurobiology 158, 132 (2017) (el primero de los cuatro papers) también es bastante sugerente:

Práctica computacional [Hoy hay práctica en el aula]

Para la práctica computacional, tienen que completar el notebook de Python que pueden encontrar [aquí]. Para los que no tienen la menor idea de Python, es recomendable que lean los siguientes tutoriales preparados por la FIFA (Federación Interestudiantil de Físicos Argentinos):

Aprovechen las herramientas de AI para generar código.

Para poder usar el notebook sin necesidad de tener Python instalado en sus computadoras, deben tener una cuenta de Google y bajarse la aplicación de Google Colaboratory. Hagan una copia del notebook en su propio drive. Después deberán compartirnos el enlace con el notebook completo.

La entrega consiste en el notebook de Python con los algoritmos completos y los resultados. La fecha límite de entrega es el  26 de junio a las 23:59, según el reloj de mi videocasetera. Envíennos el link a los cinco docentes de la práctica.

Formen grupos de tres personas, eventualmente de 2 y de 4. Si alguien no consiguió asociarse, comuníquense por mail conmigo y los pongo en contacto.

Cuando tengan el notebook funcionando, divídanse el trabajo para aprovechar el tiempo de computadora. De manera realista, tienen que poder hacer todo en menos de una semana a un ritmo tranquilo. Hay casi dos semanas entre la fecha de entrega y el segundo parcial. Aprovechen esta semana para hacer la práctica computacional. Después de hoy, no vamos a tener prácticas hasta el 24