Eigenstate Thermalization Hypothesis (ETH) for off-diagonal matrix elements in integrable spin chains

Este estudio investiga los elementos de matriz fuera de la diagonal en cadenas de espín integrables, demostrando mediante el Ansatz de Bethe algebraico que dichos elementos decaen exponencialmente con el tamaño del sistema tanto para estados dentro del mismo macroestado térmico como para estados en macroestados diferentes, siguiendo distribuciones de Gumbel pero con tasas de decaimiento distintas.

Lo esencial

¡Claro que sí! Imagina que este artículo es como un informe de detectives cuánticos que intentan entender cómo funciona el "caos" y el "orden" en el mundo de los átomos.

Aquí tienes la explicación de la investigación de Federico Rottoli y Vincenzo Alba, contada como una historia:

🎭 El Gran Misterio: ¿Cómo se "calienta" un sistema cuántico?

Imagina que tienes una habitación llena de bailarines (átomos) que nunca se tocan entre sí, pero que bailan siguiendo una música perfecta (la física cuántica). Si los dejas bailar solos durante mucho tiempo, ¿llegarán a un estado de "calma" donde todos se mueven de forma predecible, como si estuvieran en una fiesta desordenada?

En la física moderna, hay una regla famosa llamada Hipótesis de Termalización de Eigenestados (ETH). Esta regla dice que, en sistemas "caóticos" (como una fiesta real donde todos chocan y bailan desordenadamente), los bailarines individuales terminan comportándose como si estuvieran en una temperatura promedio. Es como si cada bailarín llevara un termómetro interno que siempre marca la misma temperatura, sin importar de quién sea.

🧩 El Problema: Los Sistemas "Integrables" (Los Bailarines Perfectos)

El problema es que existen sistemas especiales, llamados integrables (como la cadena de espines Heisenberg que estudian estos autores). Imagina que estos bailarines tienen una coreografía perfecta y matemática: nunca chocan de forma caótica, sino que se mueven en patrones predecibles.

En estos sistemas "perfectos", la regla ETH debería fallar. Se pensaba que los bailarines nunca se calentarían de la manera normal porque tienen demasiadas "reglas de baile" (cargas conservadas) que los mantienen en orden.

🔍 La Misión: Mirar los "Ojos" de los Bailarines

Los autores no miraron el baile completo (que es imposible de calcular para miles de bailarines). En su vez, decidieron mirar algo muy específico: los "ojos" que un bailarín tiene puestos en otro.

En física cuántica, esto se llama elementos de matriz. Imagina que cada bailarín tiene una tarjeta con un número escrito.

• Si miras al mismo bailarín (elemento diagonal), el número te dice su energía.
• Si miras de un bailarín a otro (elemento no diagonal o off-diagonal), el número te dice qué tan probable es que el bailarín A "salte" al estado del bailarín B.

La pregunta del artículo es: ¿Qué pasa con estos números cuando los sistemas son grandes?

📉 Los Descubrimientos: Dos Tipos de "Distancias"

Los autores usaron una herramienta matemática muy potente (el Ansatz de Bethe Algebraico) para calcular estos números en cadenas de hasta 240 átomos (algo que antes era imposible de calcular con ordenadores normales).

Descubrieron dos cosas fascinantes:

1. Cuando los bailarines son "vecinos" (Mismo Macroestado)

Si tomas dos bailarines que están en el mismo "estado térmico" (digamos, ambos están en una fiesta a 20°C), sus números de conexión (los elementos de matriz) desaparecen muy rápido a medida que la cadena crece.

• La analogía: Es como si intentaras susurrar un secreto a alguien en una multitud gigante. A medida que la multitud crece, tu voz se vuelve inaudible.
• El hallazgo: La probabilidad de que se conecten cae exponencialmente (como $e^{-L}$). Pero lo más curioso es la forma en que cae. No es una curva suave y aburrida (Gaussiana), sino que sigue una distribución llamada Gumbel.
  • ¿Qué es Gumbel? Imagina que lanzas 1000 monedas. La distribución normal te dice cuántas caras obtienes en promedio. La distribución Gumbel te dice: "¿Cuál fue el récord de caras seguidas que lograste?". Es una distribución de "extremos". Esto significa que en estos sistemas cuánticos, lo que importa son los casos extremos, no el promedio.

2. Cuando los bailarines son "extraños" (Diferentes Macroestados)

Ahora, imagina que tomas un bailarín de una fiesta muy fría (cercano al cero absoluto) y otro de una fiesta muy caliente (temperatura infinita).

• El hallazgo: Aquí la desconexión es aún más brutal. Los números de conexión caen mucho más rápido, como si la distancia fuera cuadrática ($e^{-L^2}$).
• La analogía: Es como intentar conectar dos personas que están en galaxias diferentes. No solo es difícil, es casi imposible. La probabilidad de que interactúen es tan pequeña que es como si no existieran en el mismo universo.

🎲 La Sorpresa: No son "Normales"

En los sistemas caóticos normales (los que siguen la ETH clásica), se esperaba que estos números de conexión fueran aleatorios pero con una distribución "normal" (la famosa campana de Gauss).

Pero los autores descubrieron que en estos sistemas integrables (los "perfectos"), la distribución es Gumbel.

• Traducción simple: Esto significa que el "caos" en estos sistemas perfectos es de un tipo muy especial. No es un ruido blanco y aburrido; es un ruido donde los eventos raros y extremos (los récords) dominan la estadística. Es como si el sistema cuántico dijera: "No me importa el promedio, solo me importan los casos más locos".

🚀 ¿Por qué es importante?

1. Validación: Confirmaron que incluso en sistemas "perfectos" y ordenados, la termalización (el proceso de llegar al equilibrio) ocurre, pero de una manera más sutil y con estadísticas diferentes a las que pensábamos.
2. Herramientas: Usaron matemáticas avanzadas (Bethe Ansatz) para ver cosas que los ordenadores normales no pueden ver (sistemas de 240 átomos). Esto es como usar un telescopio para ver estrellas que antes eran invisibles.
3. Futuro: Esto ayuda a entender cómo se comportan los materiales cuánticos reales y cómo podrían funcionar las futuras computadoras cuánticas, que a menudo se basan en sistemas que intentan ser "integrables" para mantener la información.

En resumen

El artículo nos dice que, incluso en un universo cuántico que parece tener reglas de baile perfectas y ordenadas, si miras de cerca cómo interactúan las partículas, descubres que la desconexión es brutalmente rápida y que la estadística de estas interacciones no es la que esperábamos (no es una campana normal, sino una distribución de récords extremos). Es un nuevo capítulo en la historia de cómo entendemos el equilibrio en el mundo cuántico.

Resumen técnico

1. Planteamiento del Problema

El estudio de la física cuántica de muchos cuerpos fuera del equilibrio busca entender cómo la evolución unitaria de sistemas aislados conduce a la termalización local. La Hipótesis de Termalización de Autoestados (ETH) es el marco estándar para sistemas caóticos (no integrables), donde los elementos de matriz de operadores locales entre autoestados de energía siguen una forma específica: los elementos diagonales son funciones suaves de la energía, mientras que los elementos fuera de la diagonal están suprimidos exponencialmente con el tamaño del sistema ($L$) y siguen una distribución gaussiana.

Sin embargo, en sistemas integrables (como la cadena de Heisenberg XXX), la presencia de un número extensivo de cargas conservadas viola la ETH estándar. Aunque se ha estudiado el comportamiento de los elementos diagonales en sistemas integrables, el comportamiento de los elementos de matriz fuera de la diagonal ($\langle E_i | O | E_j \rangle$) en redes de espín (modelos de red) con estados ligados (conocidos como "cuerdas" o strings en el ansatz de Bethe) permanecía poco explorado. Un trabajo previo (Ref. [26]) había estudiado esto en el modelo de Lieb-Liniger (bosones en continuo), pero no en cadenas de espín con interacciones y estados ligados, lo cual es crucial para describir macroestados de equilibrio realistas.

2. Metodología

Los autores emplean técnicas de Ansatz de Bethe Algebraico (ABA) de última generación para calcular elementos de matriz de operadores locales en la cadena de Heisenberg isotrópica de espín-1/2 (modelo XXX).

• Herramientas Computacionales: Utilizan fórmulas exactas para los elementos de matriz de operadores con soporte en uno y dos sitios. A diferencia de la diagonalización exacta (limitada a $L \sim 20-30$), el ABA permite alcanzar tamaños de sistema de hasta $L \sim 240$.
• Manejo de Singularidades: Abordan el desafío técnico de las singularidades ficticias que surgen al tomar el límite homogéneo y al tratar con soluciones complejas (cuerdas/Bethe strings). Implementan estrategias de regularización para eliminar estas singularidades, permitiendo el cálculo numérico de elementos de matriz entre autoestados genéricos.
• Muestreo de Macroestados: Utilizan un esquema de Monte Carlo (Metropolis) para muestrear autoestados dentro de macroestados termodinámicos (ensembles de Gibbs a diferentes temperaturas, incluyendo temperatura infinita y cero).
• Análisis Estadístico: Definen la variable $M_{ij} = \ln |\langle E_i | O | E_j \rangle|^2$ y analizan su distribución de probabilidad (PDF), su dependencia con la diferencia de energía y el tamaño del sistema.

3. Contribuciones Clave

1. Extensión a Modelos de Red: Se investiga por primera vez la estructura de los elementos de matriz fuera de la diagonal en un modelo de red integrable con estados ligados (cadenas de espín), más allá de los modelos de campo continuo.
2. Distinción de Macroestados: Se analiza sistemáticamente la diferencia entre elementos de matriz entre autoestados del mismo macroestado (termal) y entre macroestados diferentes (ej. estado fundamental vs. temperatura infinita).
3. Caracterización de la Distribución: Se identifica que la distribución estadística de los elementos de matriz fuera de la diagonal en sistemas integrables no es gaussiana, sino que sigue una distribución de Gumbel, contradiciendo la predicción estándar de ETH para sistemas caóticos.
4. Escalado de Tamaño Finito: Se cuantifica la dependencia exponencial de los elementos de matriz con el tamaño del sistema $L$, revelando correcciones logarítmicas específicas para el mismo macroestado y un decaimiento mucho más rápido ($\sim e^{-L^2}$) para macroestados diferentes.

4. Resultados Principales

A. Elementos de Matriz en el Mismo Macroestado

• Decaimiento Exponencial: Los elementos de matriz fuera de la diagonal decaen exponencialmente con el tamaño del sistema:
  $$|\langle E_i | O | E_j \rangle| \propto \exp\left(-c_O L \ln L - \tilde{M}_{ij} L\right)$$
  Esto confirma que la presencia de estados ligados (cuerdas) no altera cualitativamente el escenario de decaimiento observado en el modelo de Lieb-Liniger, aunque introduce correcciones logarítmicas ($L \ln L$).
• Distribución Estadística (Gumbel): La distribución de probabilidad de la variable $\tilde{M}_{ij}$ (tras un reescalado adecuado) se ajusta perfectamente a una distribución de Gumbel.
  • Esto es fundamentalmente diferente a la predicción de ETH para sistemas caóticos, donde los elementos de matriz seguirían una distribución gaussiana.
  • Se demuestra que la no-gaussianidad es robusta y detectable mediante el análisis de ratios de escala invariante ($\Gamma$), los cuales divergen con $L$, confirmando la naturaleza no gaussiana.
• Independencia del Operador: Este comportamiento se observa tanto para operadores de un espín ($E^{(11)}_i$) como para operadores de dos espines ($E^{(11)}_i E^{(22)}_{i+1}$), aunque para estos últimos el tamaño de sistema accesible es menor ($L \sim 60$) debido al costo computacional.

B. Elementos de Matriz entre Diferentes Macroestados

• Decaimiento Más Rápido: Cuando los autoestados pertenecen a macroestados termodinámicos diferentes (ej. estado fundamental vs. estado de temperatura infinita), el decaimiento es mucho más rápido:
  $$|\langle E_i | O | E_j \rangle| \propto \exp(-d_O L^2)$$
• Estadística: A pesar del cambio en la tasa de decaimiento, la distribución de probabilidad de los elementos de matriz sigue siendo descrita por una distribución de Gumbel.

C. Implicaciones sobre la Integrabilidad

Los resultados sugieren que la integrabilidad se manifiesta no solo en la violación de la termalización estándar (debido a las cargas conservadas), sino también en la estadística de las fluctuaciones de los elementos de matriz. Mientras que en sistemas caóticos las fluctuaciones son gaussianas, en sistemas integrables siguen una ley de valores extremos (Gumbel), lo que implica fluctuaciones más grandes y una estructura estadística distinta.

5. Significado e Impacto

• Validación Teórica: El trabajo proporciona evidencia numérica sólida de que el escenario de ETH modificado para sistemas integrables (propuesto teóricamente y verificado en modelos de campo continuo) es universal y se aplica a modelos de red con estados ligados.
• Conexión con Teoría de Probabilidad Libre: Los resultados abren la puerta a investigar la relación entre la estadística de los elementos de matriz en sistemas integrables y la Teoría de Probabilidad Libre (Free Probability Theory), sugiriendo que esta teoría podría ser aplicable más allá de los sistemas caóticos.
• Dinámica Cuántica: Comprender la estadística de estos elementos de matriz es crucial para reconstruir la dinámica temporal completa de observables después de un quench cuántico en sistemas integrables, permitiendo predecir la relajación y el estado estacionario con mayor precisión.
• Avance Computacional: Demuestra la viabilidad de usar el Ansatz de Bethe Algebraico para estudiar propiedades estadísticas de elementos de matriz en sistemas de muchos cuerpos de tamaño mesoscópico, superando las limitaciones de la diagonalización exacta.

En resumen, el artículo establece que, aunque los sistemas integrables no termalizan en el sentido estándar de ETH, sus elementos de matriz fuera de la diagonal exhiben un patrón de decaimiento exponencial sistemático y una estadística de Gumbel universal, diferenciándose claramente de los sistemas caóticos gaussianos.