Circular Rosenzweig-Porter random matrix ensemble

Este artículo propone y valida numéricamente un análogo unitario (circular) del ensamble de matrices aleatorias de Rosenzweig-Porter, definido mediante el movimiento browniano de Dyson, para modelar las estadísticas de niveles y la fractalidad de los autoestados de la localización de muchos cuerpos en sistemas periódicamente impulsados (Floquet).

Lo esencial

Imagina que tienes una pista de baile gigante y caótica llena de miles de bailarines. En el mundo de la física cuántica, estos bailarines representan los estados posibles de un sistema complejo (como un grupo de átomos interactuando).

Este artículo trata sobre comprender cómo se mueven e interactúan estos bailarines cuando la música cambia de dos maneras diferentes:

1. Sistemas Estáticos: La música es un zumbido constante e inmutable (como una habitación normal y tranquila).
2. Sistemas Floquet (Impulsados Periódicamente): La música es un ritmo repetitivo que sigue cambiando las reglas cada pocos segundos (como una luz estroboscópica o un láser pulsante).

Durante mucho tiempo, los físicos tuvieron un excelente "libro de reglas" (llamado el conjunto de Rosenzweig-Porter) para el primer escenario (la habitación estática). Este libro de reglas les ayuda a predecir si los bailarines se mezclarán libremente (caos) o se quedarán atrapados en sus propios rincones (localización).

Sin embargo, nadie tenía un buen libro de reglas para el segundo escenario (la habitación pulsante y rítmica). Dado que las reglas de la mecánica cuántica cambian cuando las cosas son impulsadas por un ritmo, las matemáticas antiguas no encajaban del todo.

La Nueva Idea: Una Pista de Baile Circular

Los autores de este artículo se preguntaron: "¿Podemos construir una versión de ese viejo libro de reglas que funcione para los sistemas rítmicos y pulsantes?"

Crearon un nuevo modelo al que llaman el conjunto circular de Rosenzweig-Porter.

Así es como lo construyeron, utilizando una analogía simple:

• La Vieja Forma (Movimiento Browniano): Imagina a los bailarines moviéndose aleatoriamente sobre una línea recta y plana. Si los empujas aleatoriamente con el tiempo, se dispersan de una manera predecible. Así es como funcionaba el modelo antiguo.
• La Nueva Forma (Movimiento Circular): Para los sistemas rítmicos, los autores se dieron cuenta de que los bailarines no se mueven sobre una línea recta; se mueven sobre un círculo. Imagina a los bailarines corriendo alrededor de una pista circular. Sus posiciones se miden mediante ángulos (como la esfera de un reloj) en lugar de distancias rectas.

Definieron su nuevo modelo como el resultado de una "caminata aleatoria" que ocurre específicamente en este círculo. No solo adivinaron las matemáticas; simularon este proceso en una computadora para ver qué sucedía.

Lo Que Encontraron

Los autores ejecutaron simulaciones masivas por computadora (involucrando hasta 1.000 bailarines) para ver si su nuevo modelo "circular" se comportaba como el antiguo modelo de "línea recta". Verificaron dos cosas principales:

1. El Espaciamiento de los Bailarines (Niveles de Energía)
Observaron los espacios entre los bailarines.

• En la zona "Caótica": Los bailarines están distribuidos uniformemente y los espacios entre ellos siguen un patrón específico y complejo (como una fiesta concurrida donde todos se empujan).
• En la zona "Localizada": Los bailarines se agrupan o se mantienen muy separados de una manera muy predecible y simple (como personas paradas en una fila).
• El Resultado: Su nuevo modelo circular mostró exactamente el mismo cambio de "caótico" a "agrupado" que el modelo antiguo. El "punto de inflexión" donde cambia el comportamiento ocurrió en el mismo lugar.

2. La Forma de los Bailarines (Estados Propios)
Observaron qué tan "extendida" está la influencia de un solo bailarín.

• Extendido: La energía de un bailarín se comparte entre muchos otros.
• Fractal: Un bailarín está en un medio terreno extraño: extendido, pero no completamente. Es como una nube que tiene una forma difusa y autosimilar.
• Localizado: Un bailarín está atrapado en un solo lugar.
• El Resultado: El modelo circular reprodujo exactamente estas mismas formas. Ya sea que los bailarines estuvieran completamente mezclados, parcialmente mezclados (fractales) o atrapados, el nuevo modelo coincidió perfectamente con el antiguo.

La Conclusión

El artículo afirma que han construido exitosamente una versión unitaria (circular) del famoso modelo de Rosenzweig-Porter.

Al tratar el sistema como un círculo en lugar de una línea recta, crearon una herramienta que describe con precisión el comportamiento de los sistemas cuánticos impulsados periódicamente (pulsantes). Al igual que el modelo antiguo era una herramienta "fenomenológica" (descriptiva) para sistemas estáticos, este nuevo modelo circular sirve como una herramienta descriptiva para sistemas que son sacudidos o impulsados rítmicamente.

Lo demostraron mostrando que las "huellas dactilares" estadísticas de su nuevo modelo (cómo se espacian los niveles y cómo se forman los estados) son indistinguibles de las huellas dactilares del modelo original, bien comprendido. Esto ofrece a los físicos una nueva y confiable manera de estudiar sistemas cuánticos complejos y rítmicos sin tener que resolver ecuaciones increíblemente difíciles desde cero.

Resumen técnico

Resumen Técnico: Ensemble Circular de Matrices Aleatorias Rosenzweig-Porter

Planteamiento del Problema
El ensemble de matrices aleatorias Rosenzweig-Porter (RP) es un modelo fenomenológico bien establecido utilizado para describir las estadísticas de niveles y la fractalidad de los autoestados a través de la transición de localización de muchos cuerpos (MBL) en sistemas cuánticos interactivos, estáticos y desordenados. En sistemas estáticos, la dinámica está gobernada por Hamiltonianos hermitianos. Sin embargo, los sistemas impulsados periódicamente (Floquet) se describen mediante operadores unitarios de Floquet, donde los autovalores se encuentran sobre el círculo unitario en el plano complejo. Aunque la MBL se ha extendido a sistemas de Floquet, no se había construido un análogo unitario (circular) directo del ensemble RP capaz de capturar la misma fenomenología para sistemas impulsados periódicamente. Este artículo aborda la cuestión de si tal análogo unitario puede definirse y si comparte las propiedades estadísticas clave del ensemble RP original.

Metodología
Los autores proponen una construcción del ensemble Circular Rosenzweig-Porter (CRP) basada en la observación de que el ensemble RP estándar puede considerarse como el resultado de un proceso de movimiento browniano de Dyson (DBM) de tiempo finito.

1. Construcción: Los autores definen el ensemble CRP como el resultado de una evolución estocástica de una matriz unitaria simétrica dependiente del tiempo $S(t)$ de dimensión $N$. La evolución sigue un proceso de movimiento browniano de Dyson adaptado para matrices unitarias:
   $$dS(t) = i\sqrt{dt} U^T(t) X U(t)$$
   donde $X$ es una matriz del Ensemble Ortogonal Gaussiano (GOE) muestreada en cada paso de tiempo infinitesimal, y $U(t)$ se deriva de la descomposición $S(t) = U^T(t)U(t)$.
2. Condiciones Iniciales: El proceso comienza con una matriz diagonal $S(0) = \text{diag}(e^{i\theta_1(0)}, \dots, e^{i\theta_N(0)})$, donde las fases iniciales $\theta_n(0)$ se muestrean independientemente de una distribución uniforme sobre $[-\pi, \pi)$. Esto asegura una densidad de estados uniforme, consistente con operadores de Floquet genéricos.
3. Definición del Ensemble: El ensemble CRP se define como el estado de la matriz $S(t)$ en un tiempo específico $t = \epsilon^2 / N^\gamma$, donde $\epsilon$ es la fuerza de perturbación y $\gamma$ controla la fuerza relativa de los términos diagonales y no diagonales.
4. Implementación Numérica: Los autores evalúan numéricamente el proceso estocástico utilizando un algoritmo que actualiza la matriz mediante $S(t+dt) = S(t)e^{i\sqrt{dt}X}$. Este enfoque preserva la unitariedad para pasos de tiempo finitos y evita descomposiciones matriciales explícitas. Las simulaciones se realizaron para dimensiones de matriz $N \leq 1000$, promediando sobre al menos $10^6$ autovalores o autoestados.

Contribuciones Clave

• Definición del CRP: El artículo introduce el ensemble Circular Rosenzweig-Porter como un análogo unitario del ensemble RP estándar, definido explícitamente mediante un proceso de movimiento browniano de Dyson circular.
• Mapeo Teórico: Los autores proporcionan un argumento heurístico de que las estadísticas espectrales del ensemble CRP en el medio del espectro deben coincidir con las del ensemble RP, siempre que la densidad de estados se escale apropiadamente (específicamente, mapeando el parámetro RP $\epsilon$ a $\epsilon/\sqrt{2\pi}$ en el contexto CRP).
• Verificación Numérica: El estudio proporciona evidencia numérica exhaustiva que caracteriza las propiedades estadísticas de los autovalores y autoestados del ensemble CRP.

Resultados
La investigación numérica confirma que el ensemble CRP exhibe comportamientos estadísticos análogos al ensemble RP a través de diferentes regímenes del parámetro $\gamma$:

1. Estadísticas de Espaciado de Niveles: La relación promedio de espaciados de niveles consecutivos, $r$, sirve como diagnóstico para el caos cuántico.
   • Para $\gamma < 2$, el sistema exhibe estadísticas de Wigner-Dyson ($r \approx 0.530$), características de sistemas caóticos.
   • Para $\gamma \geq 2$, el sistema transita a estadísticas de Poisson ($r \approx 0.386$), características de sistemas integrables.
   • Ocurre una transición de fase en el punto crítico $\gamma_c = 2$ con un exponente crítico $\nu = 1$. Los datos muestran un colapso de escala de $r$ en función de $(\gamma - \gamma_c) \ln(N)^{1/\nu}$, consistente con el ensemble RP.
2. Fractalidad de Autoestados: Se analizó la relación de participación inversa (IPR) de los autoestados en la base donde $S(0)$ es diagonal.
   • Para $\gamma \leq 1$, los autoestados están extendidos ($IPR_2 \sim N^{-1}$).
   • Para $\gamma \geq 2$, los autoestados están localizados ($IPR_2 \sim O(1)$).
   • En el régimen intermedio $1 < \gamma < 2$, los autoestados son fractales, escalando como $IPR_2 \sim N^{-(2-\gamma)}$.
3. Forma de Autoestados: Se encontró que la distribución de las amplitudes de los autoestados $|\langle m | \psi_n \rangle|^2$ sigue estadísticas de Breit-Wigner. Los datos se alinean con la predicción teórica:
   $$|\langle m | \psi_n \rangle|^2 \sim \frac{1}{(\theta_n - \theta_m^{(0)})^2 + \Gamma^2}$$
   donde $\Gamma$ es el ancho de dispersión. Aunque los efectos de tamaño finito y las aproximaciones de continuo conducen a discrepancias cuantitativas leves, el acuerdo cualitativo y las diferencias estructurales distintivas entre los regímenes $\gamma = 0.5$ y $\gamma = 1.5$ se observan claramente.

Significado y Afirmaciones
Los autores afirman que el ensemble Circular Rosenzweig-Porter sirve como un modelo fenomenológico válido para las estadísticas de niveles y la fractalidad de autoestados observadas a través de la transición de localización de muchos cuerpos en sistemas impulsados periódicamente (Floquet). Al establecer que el ensemble CRP comparte propiedades estadísticas clave con el ensemble RP, el trabajo proporciona una herramienta teórica manejable para estudiar fases no ergódicas en sistemas cuánticos impulsados.

El artículo enmarca modestamente este trabajo como un "primer paso" hacia la construcción de generalizaciones más amplias, como modelos que incorporan multifractalidad u otros modelos de Floquet con autoestados multifractales. Los autores sugieren que futuras generalizaciones podrían lograrse considerando procesos estocásticos con incrementos correlacionados, similares a propuestas recientes para el ensemble RP estándar. Además, señalan la utilidad potencial del ensemble CRP como un bloque de construcción no máximamente aleatorio en estudios de circuitos cuánticos aleatorios.