Abriendo el espacio de Krylov para acceder a la dinámica de todos los tiempos mediante simetrías dinámicas

Información Básica

• ID del Artículo: 2503.07403
• Título: Opening Krylov space to access all-time dynamics via dynamical symmetries
• Autores: Nicolas Loizeau, Berislav Buča, Dries Sels
• Instituciones: Niels Bohr Institute (Copenhague), Université Paris-Saclay, University of Oxford, New York University, Flatiron Institute
• Clasificación: quant-ph (Física Cuántica)
• Fecha de Publicación: 22 de agosto de 2025 (versión v2)
• Enlace del Artículo: https://arxiv.org/abs/2503.07403v2

Resumen

Este artículo aborda los problemas fundamentales de la dinámica a corto y largo plazo en sistemas cuánticos de muchos cuerpos cerrados. Para sistemas cerrados con interacciones locales, la dinámica de observables locales puede expandirse en modos propios (pseudo-locales) del Liouvilliano, denominados simetrías dinámicas. Estos operadores se dividen en dos categorías: operadores transitorios que decaen con el tiempo y operadores eternos que oscilan permanentemente o permanecen invariantes (leyes de conservación). El artículo propone un método para derivar estas simetrías dinámicas numérica y analíticamente mediante la introducción de condiciones de frontera abiertas que emergen naturalmente en la cadena de Krylov, dividiendo el espacio de Krylov en grados de libertad del sistema y del ambiente, donde los operadores no-locales constituyen un baño efectivo para los operadores locales.

Antecedentes de Investigación y Motivación

Problemas Centrales

1. Desafíos Computacionales en Dinámica Cuántica de Muchos Cuerpos: Resolver la dinámica a corto y largo plazo en sistemas cuánticos de muchos cuerpos cerrados es uno de los principales desafíos en física atómica y física de materia condensada
2. Dificultad en la Extracción de Simetrías Dinámicas: Aunque teóricamente la dinámica de observables locales puede caracterizarse completamente mediante modos propios del Liouvilliano (simetrías dinámicas), la derivación de estos operadores ha sido imposible excepto para muy pocos modelos

Importancia del Problema

• Caracterización Dinámica Completa: Las simetrías dinámicas proporcionan una caracterización completa de la dinámica del sistema, para cualquier operador O y estado inicial ρ, su evolución temporal puede expresarse como: Tr(O(t)ρ) = ∑_ω e^(iωt)μ_ω
• Comprensión de Fenómenos Fuera del Equilibrio: Es crucial para entender localización de muchos cuerpos, cristales de tiempo, cicatrices cuánticas, fragmentación y otros comportamientos dinámicos no triviales
• Mecanismos de Termalización: En sistemas caóticos, relaciona la hipótesis de crecimiento de operadores correlacionados con la termalización y el decaimiento exponencial de observables

Limitaciones de Métodos Existentes

1. Restricción de Sistemas Finitos: El espectro del Liouvilliano de un sistema cerrado finito es puramente real, lo que impide explorar propiedades del límite termodinámico
2. Métodos de Expansión Pseudo-Modal: Los métodos anteriores (Refs 96,97) se limitan a dinámicas cercanas al equilibrio
3. Pseudo-Localidad No Estudiada: Investigaciones previas no examinaron la (pseudo-)localidad de las simetrías dinámicas, que es la clave para afectar la física fuera del equilibrio

Motivación de la Investigación

Proponer un método más riguroso que implemente una descomposición sistema-ambiente mediante la introducción de condiciones de frontera abiertas en la cadena de Krylov, permitiendo:

• Que cantidades locales se soporten en el sistema
• Que cantidades no-locales constituyan el ambiente
• Identificar cuáles simetrías dinámicas son pseudo-locales y por lo tanto relevantes

Contribuciones Principales

1. Método Novedoso de Condiciones de Frontera Abiertas: Propone introducir condiciones de frontera disipativas en una posición específica L de la cadena de Krylov, evitando rebotes no-físicos de la función de onda, con la forma: ∂t φ_L = (b_L + b{L+1})φ_ - 2b_{L+1}φ_L
2. Marco de Descomposición Sistema-Ambiente: Establece una descomposición natural del espacio de Krylov donde:
   • Lado izquierdo (sistema): operadores locales simples con coeficientes de salto exactamente conocidos
   • Lado derecho (ambiente): operadores no-locales complejos que constituyen un baño efectivo
3. Soluciones Analíticas para Casos Ideales: Proporciona soluciones exactas para dos casos idealizados:
   • Crecimiento lineal b_n = n: vectores propios expresados en polinomios de Meixner
   • Crecimiento de raíz cuadrada b_n = √n: vectores propios expresados en polinomios de Hermite
4. Conexión Directa entre Crecimiento de Operadores y Termalización: Demuestra que en el caso de crecimiento lineal (sistemas caóticos), todos los observables decaen exponencialmente con velocidad -2α, donde α es la tasa de crecimiento del coeficiente de Lanczos
5. Verificación Numérica y Aplicaciones: Verifica la efectividad del método en sistemas de espín reales como cadenas XXZ e Ising con campo transversal, extrayendo exitosamente simetrías dinámicas

Explicación Detallada del Método

Definición de la Tarea

Entrada:

• Hamiltoniano H
• Operador inicial O_0 (o matriz de densidad inicial ρ)

Salida:

• Simetrías dinámicas A_ω que satisfacen H, A_ω = -ωA_ω
• Frecuencias ω (pueden ser complejas): la parte real corresponde a oscilaciones, la parte imaginaria a decaimiento

Restricciones:

• Enfoque solo en simetrías dinámicas pseudo-locales (que afectan la física local)
• Aplicable a sistemas con interacciones locales y estados iniciales con correlaciones de corto alcance

Arquitectura del Modelo

1. Construcción del Espacio de Krylov (Algoritmo de Lanczos)

Comenzando desde el operador semilla O_0, se construye recursivamente una base ortonormal {O_n} aplicando el Liouvilliano L = H, ·:

O_1 = LO_0/b_1 = [H,O_0]/b_1,  b_1 = ||LO_0||

Para n > 2:
O'_n = LO_{n-1} - b_{n-1}O_{n-2}
O_n = O'_n/b_n
b_n = ||O'_n||

donde la norma se define como ||O||² = 1/(2N) TrO²

2. Mapeo de Partícula Única Unidimensional

La evolución temporal en la base de Krylov se mapea a un problema de salto de partícula única 1D:

O(t) = 1/(2N) ∑_n i^n φ_n(t)O_n

La ecuación de evolución de los coeficientes: ∂t φ_n = b_n φ - b_{n+1} φ_{n+1}, φ_n(0) = δ_

3. Condiciones de Frontera Abiertas

Se trunca la cadena en la posición L, asumiendo que φ_n es suficientemente suave para aproximarse localmente como una función lineal: φ_{L+1} ≈ φ_L + (φ_L - φ_)

Obteniéndose la condición de frontera disipativa: ∂t φ_L = (b_L + b{L+1})φ_ - 2b_{L+1}φ_L

4. Matriz de Liouvilliano Modificada

El Liouvilliano truncado tiene forma tridiagonal más correcciones de frontera:

L̃ = i × [matriz tridiagonal, última fila de forma especial]
Última fila: [0, 0, ..., b_L + b_{L+1}, -2b_{L+1}]

Esta matriz no-hermítica permite valores propios complejos

Puntos de Innovación Técnica

1. Significado Físico de las Condiciones de Frontera

• Evita Introducción de Parámetros: A diferencia de Refs 96,97 con disipación externa + extrapolación, este método no requiere parámetros adicionales
• Hipótesis de Suavidad: Basada en la hipótesis natural de continuidad de derivadas de la función de onda de Krylov
• Separación Localidad-No-Localidad: Los operadores en el lado derecho de la cadena de Krylov están dominados principalmente por cadenas de Pauli k-locales con k>L

2. Diferencias con Métodos Tradicionales

• Frontera de Dirichlet: φ_{L+1} = 0 conduce a rebotes no-físicos
• Este Método: Permite flujo hacia el ambiente y retorno, reflejando más fielmente la propagación de operadores

3. Justificación del Diseño

En casos ideales (b_n = n o √n):

• La dinámica bajo condiciones de frontera abiertas es prácticamente indistinguible de la cadena infinita (Figura 2)
• El espectro real permanece casi invariante, solo la parte imaginaria introduce decaimiento físico (Figura 3)
• Caso de crecimiento lineal: todas las raíces no-triviales satisfacen Im(ω) = -2i (verificado numéricamente hasta n=2^12)

Configuración Experimental

Sistemas Modelo

1. Cadena XXZ (Modelo Integrable)

H = ∑_i (s^x_i s^x_{i+1} + s^y_i s^y_{i+1} + Δs^z_i s^z_{i+1} + hs^z_i)
Parámetros: Δ = -1/2, h = 2
Operador inicial: O_0 = Q_3 = ∑_i (s^+_i s^+_{i+1}s^+_{i+2} + s^-_i s^-_{i+1}s^-_{i+2})

Se conoce que existe una simetría dinámica con frecuencia ω = 12

2. Cadena Caótica

H = ∑_i (s^x_i s^x_{i+1} - 1.05s^z_i + 0.5s^x_i)
Operador inicial: O_0 = ∑_i (1.05s^x_i s^x_{i+1} + s^z_i)

Se elige un operador sin solapamiento con H para excluir cantidades conservadas triviales

Implementación Numérica

• Herramientas: Paquete Julia PauliStrings.jl (utiliza representación de cadenas de Pauli para cálculo eficiente de conmutadores anidados)
• Invariancia Traslacional: Solo se almacena la celda unitaria en lugar del operador extensivo
• Truncamiento de Pauli Mejorado: En cada paso del algoritmo de Lanczos se retienen M cadenas de Pauli con mayor peso (Algoritmo 1)

Métodos de Comparación

1. Simulación de Operador de Producto Matricial (MPO): Como referencia
2. Frontera de Dirichlet Estándar: Demuestra rebotes no-físicos
3. Método de Teretenkov et al.: Extrapolación lineal + disipación artificial (comparación en material suplementario)

Indicadores de Evaluación

• Precisión de Frecuencia de Simetría Dinámica: Comparación con resultados conocidos (XXZ) o predicciones analíticas (caótica)
• Fidelidad de Evolución Temporal: TrO_0(0)O_0(t) vs resultados exactos/MPO
• Medida de Localidad: Distribución del vector propio en la base de Krylov |ψ_0|² (más localizado más oscuro)

Resultados Experimentales

Resultados Principales

1. Verificación de Casos Ideales (Figuras 2-3)

Crecimiento lineal b_n = n:

• La dinámica bajo condiciones de frontera abiertas es indistinguible de la cadena infinita
• Todos los valores propios satisfacen Im(ω) = -2i (excepto dos raíces triviales ω = -i, -3i)
• La frontera de Dirichlet causa rebotes no-físicos significativos en t≈5

Crecimiento de raíz cuadrada b_n = √n:

• Las condiciones de frontera abiertas reproducen perfectamente la dinámica de la cadena infinita
• La parte real del espectro permanece casi invariante, la parte imaginaria introduce pequeño decaimiento

2. Resultados de la Cadena XXZ (Figuras 4a,c)

• Recuperación de Simetría Dinámica: Extracción exitosa de la simetría dinámica conocida cerca de ω = 12
  • L=26: Re(ω) ≈ 12.0, Im(ω) ≈ -0.5
  • Truncamiento de Pauli mejorado (L=50): mejora significativa en convergencia
• Evolución Temporal: TrO_0(0)O_0(t) muestra oscilaciones persistentes (coincide con MPO hasta t≈3)
• Localidad: Solo la simetría dinámica de valor real es local (|ψ_0|² concentrado en n pequeño)

3. Resultados de Cadena Caótica (Figuras 4b,d)

• Sin Simetrías Eternas: Todos los modos decaen, sin frecuencias puramente reales
• Agrupamiento de Tasa de Decaimiento: Los valores propios se concentran cerca de Im(ω) ≈ -0.72i
• Consistencia con Predicción Teórica: -0.72 ≈ -2λ, donde λ es la tasa de crecimiento del coeficiente de Lanczos
• Decaimiento Exponencial: TrO_0(0)O_0(t) decae rápidamente a 0 (coincide con MPO)
• No-Localidad: Todos los modos son extremadamente no-locales (|ψ_0|² distribuido uniformemente)

4. Templado desde Estado Estacionario (Figura 5)

Templado desde |ψ_0⟩ = 2^{-N/2}(|0⟩+|1⟩)^⊗N:

• Q_1 = ∑_i (s^+_i + s^-_i): decaimiento rápido
• Q_5 (expansión de simetría dinámica 3-local 5-orden): oscilaciones de largo tiempo
• Resultados para diferentes L (24-30) convergen bien
• El truncamiento estándar (línea gris) muestra comportamiento no-físico

Experimentos de Ablación

Impacto del Truncamiento de Pauli (Figura Suplementaria 4)

• Sin Truncamiento de Pauli: L_max = 26 (límite de cálculo exacto)
• Con Truncamiento de Pauli (retención de 2^22 cadenas): L_max = 50
  • Frecuencia de simetría dinámica: Re(ω) mejora de ≈12.1 a ≈12.05
  • Im(ω) mejora de ≈-0.8 a ≈-0.2

Análisis de Convergencia (Figuras Suplementarias 2-3)

• Modelo XXZ: Re(ω) se estabiliza después de L=10, Im(ω) mejora continuamente
• Cadena Caótica: El espectro converge básicamente después de L=30, los valores propios de L profundo se agrupan densamente cerca de la línea Im(ω)=-2λ

Comparación con Método de Teretenkov (Figuras Suplementarias 6-7)

Cadena Caótica:

• Los valores propios del método de Teretenkov (γ=0.2-0.4) no convergen al resultado analítico
• Este método reproduce precisamente el agrupamiento Im(ω)=-2λ

Cadena XXZ:

• Ambos métodos tienen precisión similar para simetrías dinámicas locales (ω=12)
• Las predicciones de modos no-locales son significativamente diferentes

Hallazgos Experimentales

1. Relación Cuantitativa entre Crecimiento de Operadores y Termalización: En el caso de crecimiento lineal (caótico), la tasa de decaimiento es exactamente -2α, verificando directamente la predicción teórica
2. Criterio de Localidad: |ψ_0|² en la representación de Krylov proporciona automáticamente una medida cuantitativa de la localidad de la simetría dinámica
3. Aplicabilidad del Método:
   • Sistemas caóticos (crecimiento lineal): excelente desempeño, la hipótesis de suavidad es válida
   • Sistemas integrables (crecimiento de raíz cuadrada): puede extraer simetrías dinámicas, pero convergencia más lenta
4. Cuello de Botella Computacional: La limitación principal es el cálculo de coeficientes de Lanczos de alto orden (los operadores se vuelven altamente no-locales), el truncamiento de Pauli proporciona alivio efectivo

Trabajo Relacionado

Métodos de Espacio de Krylov

• Complejidad de Operadores: Parker et al. 105 para explorar crecimiento de operadores en caos cuántico
• Hidrodinámica: Refs 97, 117-119 estudian propiedades de transporte
• Sistemas Floquet: Refs 120-124 aplicados a sistemas impulsados periódicamente

Expansión Pseudo-Modal

• Proyección de Mori: Método temprano para calcular funciones de correlación temporal 92
• Desarrollos Recientes: Refs 96, 97 sugieren agregar disipación + extrapolación o estudiar resonancias de Ruelle-Pollicott 98-104
• Limitaciones: Limitado a cerca del equilibrio, no examina pseudo-localidad

Teoría de Simetrías Dinámicas

• Equilibrio Eterno: Marco teórico de Buča 82, forma de estado ρ(t) = Z^{-1}exp(∑_u μ_u e^{iλ_u t}A_u)
• Álgebras Generadoras de Espectro: Refs 33, 34, 90, 91 relacionadas pero no necesariamente pseudo-locales
• Cristales de Tiempo y Cicatrices: Refs 12, 29-44 aplicaciones de simetrías dinámicas

Ventajas de Este Trabajo

1. Sin Parámetros: No introduce tasas de disipación artificial u otros parámetros adicionales
2. Sin Extrapolación: Utiliza directamente coeficientes de Lanczos conocidos
3. Pseudo-Localidad: Identifica explícitamente simetrías pseudo-locales relevantes
4. Lejos del Equilibrio: Aplicable a situaciones fuera del equilibrio como templado fuerte

Conclusiones y Discusión

Conclusiones Principales

1. Validez del Método: Las condiciones de frontera abiertas logran exitosamente el cálculo numérico del espectro del Liouvilliano de sistemas cerrados infinitos
2. Perspectivas Teóricas:
   • Sistemas caóticos: crecimiento lineal de Lanczos → todos los observables decaen con velocidad -2α
   • Sistemas integrables: crecimiento de raíz cuadrada → existen modos de oscilación eterna
3. Valor Práctico: Combinado con métodos de cadenas de Pauli, proporciona una herramienta poderosa para calcular dinámicas de muchos cuerpos

Limitaciones

1. Hipótesis de Suavidad: El método depende de la suavidad de φ_n en la representación de Krylov
   • Sistemas caóticos: el crecimiento lineal garantiza suavidad
   • Otros casos: requiere verificación
2. Complejidad Computacional:
   • Operadores de alto orden O_n se vuelven extremadamente no-locales
   • Modelo XXZ solo puede calcularse exactamente hasta L=26
   • El truncamiento de Pauli introduce aproximación pero es efectivo
3. Precisión de Largo Tiempo:
   • La parte imaginaria de la simetría dinámica no converge completamente causando error acumulativo a largo tiempo
   • Requiere L más grande o estrategia de truncamiento mejorada
4. Dificultad de Templado desde Estado Térmico: El templado desde estado térmico ρ = e^{-βO_0} requiere calcular la exponencial de sumas de operadores no-conmutativos, numéricamente difícil (solo puede manejarse templado desde estado estacionario)

Direcciones Futuras

1. Circuitos Cuánticos: Generalización a casos donde el Liouvilliano es matriz de Hessenberg superior 107
2. Otras Simetrías Pseudo-Locales:
   • Simetrías dinámicas semi-locales 137-139
   • Simetrías pseudo-locales de cicatrices cuánticas 30, 31, 36, 82
3. Simulación de Experimentos RMN: Cálculo directo de la transformada de Fourier de Tr(Z_tot(t)Z_tot(0)) 140
4. Mejora de Convergencia:
   • Estrategias de truncamiento más refinadas
   • Truncamiento de Pauli adaptativo
   • Optimización de algoritmo recursivo (material suplementario)

Evaluación Profunda

Fortalezas

1. Innovación del Método ⭐⭐⭐⭐⭐

• Condiciones de Frontera Originales: Basadas en la hipótesis física natural de suavidad de la función de onda, diferente de la adición artificial de disipación
• Descomposición Sistema-Ambiente: Marco ingenioso para implementar separación localidad-no-localidad a nivel de espacio de operadores
• Combinación Teoría-Numerología: Las soluciones analíticas en casos ideales (polinomios de Meixner/Hermite) se complementan con verificación numérica en sistemas reales

2. Profundidad Teórica ⭐⭐⭐⭐⭐

• Conexión Crecimiento de Operadores-Termalización: Primera demostración cuantitativa de la relación universal crecimiento lineal → decaimiento -2α
• Cuantificación de Pseudo-Localidad: Identificación automática de simetrías dinámicas relevantes mediante |ψ_0|²
• Integración con Marco de Equilibrio Eterno: Integración profunda con teoría de Buča 82

3. Suficiencia Experimental ⭐⭐⭐⭐

• Verificación Multi-Modelo: Dos casos típicos: integrable (XXZ) y caótica
• Múltiples Comparaciones: MPO, frontera de Dirichlet, método de Teretenkov
• Análisis de Convergencia: Estudio detallado de dependencia en L y truncamiento de Pauli
• Insuficiencia: Falta prueba en más sistemas físicos reales (como modelo de Hubbard)

4. Implementación Técnica ⭐⭐⭐⭐⭐

• Código Abierto: Implementación en Julia disponible (GitHub + PauliStrings.jl)
• Optimizaciones de Algoritmo: Invariancia traslacional, truncamiento de Pauli, algoritmo recursivo y otros trucos prácticos
• Reproducibilidad: Material suplementario detallado y configuración de parámetros

Insuficiencias

1. Rigor de la Hipótesis de Suavidad

• Falta Prueba Teórica: ¿Cuándo φ_n es suficientemente suave para aproximación lineal?
• Casos de Fallo: Puede fallar bajo fragmentación fuerte o condiciones de frontera especiales
• Recomendación: Requiere argumentación matemática más rigurosa o criterios de aplicabilidad

2. Escalabilidad Computacional

• Límite de L: XXZ solo alcanza 26 (50 requiere truncamiento), limitando aplicaciones de alta precisión
• Dependencia de Tamaño de Sistema: No se discute la dependencia en N (tamaño del sistema)
• Generalización 2D: Solo cadenas 1D, sistemas de dimensión superior son difíciles

3. Verificación Experimental Limitada

• Tipo de Templado: Solo templado desde estado estacionario, no se implementó templado desde estado térmico
• Escala de Tiempo: Precisión de largo tiempo limitada (t<5 para XXZ)
• Cantidades Físicas: Principalmente auto-correlación de operadores, falta de otros observables

4. Comparación Insuficiente con Método de Teretenkov

• Solo en Material Suplementario: No se discute detalladamente en texto principal
• Sensibilidad a Parámetros: No se estudia sistemáticamente el impacto de γ y esquema de extrapolación
• Análisis de Ventajas/Desventajas: ¿Cuándo este método es superior? ¿Cuándo son equivalentes?

Impacto

Contribución al Campo ⭐⭐⭐⭐⭐

• Cambio de Paradigma: Introduce ideas de sistemas abiertos al estudio de dinámica de sistemas cerrados
• Herramienta Universal: Aplicable a amplia gama de problemas de muchos cuerpos cuánticos
• Fundamento Teórico: Proporciona marco unificado para entender termalización, cristales de tiempo, cicatrices, etc.

Valor Práctico ⭐⭐⭐⭐

• Método Numérico: Puede aplicarse directamente al cálculo de dinámicas de sistemas reales
• Guía Experimental: Apoyo teórico para experimentos como RMN
• Limitación: El costo computacional aún limita aplicación a sistemas muy grandes

Reproducibilidad ⭐⭐⭐⭐⭐

• Código Abierto: Implementación completa en GitHub
• Documentación Detallada: Material suplementario contiene todos los detalles de algoritmos
• Transparencia de Parámetros: Todas las configuraciones experimentales claramente especificadas

Escenarios de Aplicabilidad

Escenarios Ideales

1. Sistemas 1D con Interacciones Locales: Cadenas de espín, bosones en red, etc.
2. Sistemas Caóticos: El crecimiento lineal de Lanczos garantiza efectividad del método
3. Dinámica de Corto-Medio Plazo: t<10τ (escala de tiempo característica)
4. Observables Locales: Medición de operadores k-locales

Escenarios que Requieren Precaución

1. Sistemas Fuertemente Fragmentados: La hipótesis de suavidad puede fallar
2. Interacciones de Largo Alcance: La estructura del espacio de Krylov puede ser diferente
3. Dinámica Extremadamente Larga: Convergencia incompleta de parte imaginaria causa error acumulativo
4. Sistemas de Dimensión Alta: El costo computacional aumenta dramáticamente

Escenarios No Aplicables

1. Sistemas Abiertos: Marcos de Lindbladian más apropiados
2. Sistemas con Medición Fuerte: Requiere herramientas teóricas diferentes
3. Límite Clásico: Métodos clásicos más simples

Suplemento de Detalles Técnicos

Derivación de Fórmulas Clave

Forma de Ecuación Principal de Condiciones de Frontera (Material Suplementario S20)

La evolución temporal exacta Ȯ = iH,O, la condición de frontera dura corresponde a:

Ȯ = iH,O - 1/(2N) Tr(O_L† iH,O)O_L

Sustituyendo la condición de frontera abierta:

Ȯ = iH,O - 1/(2N) Tr(O_L† iH,O)O_L + ib_{L+1}Tr(O_† O)O_L - 2b_{L+1}Tr(O_L† O)O_L

Demuestra cómo la condición de frontera equilibra el flujo de operadores hacia el ambiente

Solución de Polinomios de Meixner (Material Suplementario S3-S7)

Cuando el crecimiento es lineal b_n = n, el vector propio:

φ_n(ω) = i^n M_n(ω)/n!

donde M_n satisface la recursión: M_{n+1}(x) = xM_n(x) - n²M_(x)

La condición de frontera da la restricción de frecuencia: (ω + 2i(L+1))M_L(ω) = (2L+1)LM_(ω)

Las raíces triviales ω = -i, -3i corresponden a φ_n ∝ 1 y φ_n ∝ 2n-1

Trucos Numéricos

Algoritmo Recursivo (Material Suplementario Algoritmo)

Mejora en la extracción de simetrías dinámicas:

1. Inicialización aleatoria de O_0 k-local
2. Diagonalización de L̃, selección de modo propio con frecuencia ≈ω
3. Truncamiento de A_ω a k-local, actualización O_0 ← truncate(A_ω)
4. Repetición hasta convergencia

En modelo XXZ con k=3, converge en 4 pasos: Re(ω) de 12.14→12.02

Referencias (Seleccionadas)

Fundamentos Metodológicos

• 82 B. Buča, PRX 13, 031013 (2023) - Teoría de equilibrio eterno
• 105 D. Parker et al., PRX 9, 041017 (2019) - Complejidad de Krylov
• 96,97 A. Teretenkov, O. Lychkovskiy et al., PRB (2024) - Expansión pseudo-modal

Campos de Aplicación

• 12 B. Buča et al., Nat. Commun. 10, 1730 (2019) - Cristales de tiempo
• 29-44 Serie de literatura sobre cicatrices cuánticas
• 119 N. Loizeau et al., SciPost (2025) - PauliStrings.jl

Fundamento Teórico

• 1 L. D'Alessio et al., Adv. Phys. 65, 239 (2016) - Revisión de termalización
• 86-89 Teoría de pseudo-localidad (Prosen, Doyon, etc.)

Evaluación General: Este es un artículo excelente con innovación importante en el campo de dinámica cuántica de muchos cuerpos (calificación 4.5/5). El método central (cadena de Krylov con frontera abierta) es elegante y simple con imagen física clara, la derivación teórica es rigurosa (soluciones analíticas en casos ideales), y la verificación numérica es suficiente. Las principales contribuciones radican en establecer la conexión cuantitativa entre crecimiento de operadores y termalización, y proporcionar una herramienta numérica práctica. Las limitaciones principales están en que el rango de aplicabilidad de la hipótesis de suavidad requiere mayor clarificación, y la escalabilidad computacional es limitada. Para investigadores que estudian dinámicas fuera del equilibrio de sistemas cuánticos de muchos cuerpos, cristales de tiempo, cicatrices cuánticas y otros problemas de frontera, este es un artículo de lectura obligatoria.