Diameter truncated operator evolution

Resumen Técnico: Evolución de Operadores con Truncamiento de Diámetro (DTOE)

Planteamiento del Problema
Simular la dinámica de sistemas cuánticos de muchos cuerpos fuera del equilibrio es fundamentalmente desafiante debido al rápido crecimiento de la complejidad. A medida que estos sistemas evolucionan, el entrelazamiento, la no-estabilicidad y la no-gaussianidad aumentan, causando que los operadores inicialmente locales se vuelvan altamente no locales y "difusos" (scrambled). Este crecimiento hace que las simulaciones tradicionales de redes de tensores sean intratables más allá de escalas de tiempo muy cortas, ya que no pueden capturar eficientemente la evolución dinámica completa de la función de onda o de todos los observables posibles. Si bien la termalización tardía permite descripciones hidrodinámicas, el acercamiento al equilibrio en tiempos cortos e intermedios sigue siendo difícil de computar con precisión.

Metodología
Los autores introducen la Evolución de Operadores con Truncamiento de Diámetro (DTOE), un método numérico diseñado para simular la dinámica de operadores restringiendo la simulación a un subespacio de operadores con una extensión espacial acotada.

• Concepto Central: A diferencia de métodos anteriores como la Evolución de Operadores Asistida por Disipación (DAOE) o la Propagación de Pauli, que truncan los operadores basándose en su peso (el número de términos de Pauli no triviales en una expansión de cadenas), DTOE trunca basándose en el diámetro. El diámetro se define como el tamaño de la región en el espacio real en la red donde el operador tiene un soporte no trivial.
• Implementación: El método se formula dentro del marco de circuitos cuánticos de tipo brickwork en 1+1D. Utiliza una notación "plegada" donde los operadores se vectorizan en un espacio de Hilbert doble. La evolución temporal se computa aplicando matrices de transferencia que mapean intervalos de operadores de diámetro-$d$.
• Esquema de Truncamiento: Se emplea un truncamiento estricto: todos los operadores con un diámetro que exceda un corte fijo $d$ son omitidos estrictamente de la simulación. Las funciones de correlación se reconstruyen mediante la concatenación de estas matrices de transferencia para generar caminos entre operadores locales.
• Motivación Teórica: El enfoque se basa en la observación de que, si bien la mayor parte del flujo de la norma del operador se dirige hacia trayectorias altamente no locales, las funciones de correlación de dos puntos locales suelen estar dominadas por trayectorias "atípicas" que permanecen dentro de una región de diámetro constante y acotado. Para circuitos genéricos, existe una fase "fácil" en el espacio-tiempo (típicamente por encima de una velocidad crítica $v_c$) donde las correlaciones están dominadas por operadores de diámetro constante $O(1)$. En este régimen, DTOE ofrece una complejidad de tiempo polinomial ($O(t^2)$) con un escalamiento de memoria constante $O(t)$.

Contribuciones Clave y Resultados
El artículo demuestra la eficacia de DTOE mediante extensas simulaciones numéricas tanto en sistemas integrables como no integrables, sin el uso de recursos de computación de alto rendimiento.

1. Circuitos Dual-Unitary Perturbados (Modelo de Ising con Impulso):

   • Los autores simularon el modelo de Ising con impulsos con parámetros $J=B=\pi/4+\eta$ y $h=1.2$, variando la fuerza de la perturbación $\eta \in [0, 0.5]$.
   • Computaron funciones de correlación de dos puntos a temperatura infinita ($C_{XX}$) hasta $t=100$ capas de circuito.
   • Resultos: Para $\eta = 0.24588$, las aproximaciones de bajo diámetro ($d=1, 2$) divergieron de los datos exactos, pero diámetros mayores ($d \ge 9$) lograron errores relativos del orden de $10^{-2}$ a $10^{-3}$. El método se mantuvo efectivo incluso en un régimen altamente no perturbativo, superando las simulaciones actuales de computadores cuánticos digitales que típicamente gestionan $\sim 25$ ciclos de Floquet.

2. Modelo Heisenberg XXZ:

   • El método se aplicó a circuitos XXZ para extraer los exponentes de transporte hidrodinámico ($\alpha$) para el escalamiento del ancho de correlación $\sigma(t) \sim t^\alpha$.
   • Régimen Balístico ($|\Delta| < 1$): DTOE con $d=10$ extrajo con precisión $\alpha = 1.00 \pm 0.01$. Los autores observaron que la escala de tiempo $t^*$ sobre la cual el método permanece preciso antes de derivar hacia la difusión escala exponencialmente con el diámetro $d$ ($t^* \sim O(\exp(d))$).
   • Régimen Superdifusivo ($|\Delta| = 1$): En el punto isotrópico, el método capturó el exponente predicho teóricamente $\alpha = 2/3$ (específicamente $\alpha \in [0.6663, 0.6684]$ para $d=10$), a pesar de la dificultad conocida de capturar largas cadenas de cuasipartículas en este régimen.
   • Régimen Difusivo ($|\Delta| > 1$): El método reprodujo con precisión el exponente difusivo $\alpha = 1/2$ incluso con diámetros bajos, ya que el transporte está dominado por cadenas de cuasipartículas cortas en este régimen.

Significancia y Reivindicaciones
El artículo afirma que DTOE proporciona un protocolo de truncamiento más "esbelto" en comparación con los métodos basados en el peso. El análisis teórico sugiere que para una red de tamaño $L$, el número de operadores con diámetro $d$ escala como $O(L)$, mientras que el número de operadores con peso $d$ escala como $O(L^d)$. Por consiguiente, para $L$ grandes, la razón de operadores retenidos en DTOE frente al truncamiento basado en el peso se desvanece, lo que sugiere que DTOE es intrínsecamente más eficiente.

Los autores sostienen que DTOE está físicamente bien motivado para sistemas donde las correlaciones locales están dominadas por trayectorias de diámetro acotado. Destacan que el método extrae con éxito propiedades de transporte y funciones de correlación en escalas de tiempo que superan las capacidades del hardware cuántico actual. Además, señalan que aunque sus resultados se presentan en 1+1D, la literatura existente sugiere que en circuitos genéricos en 2D, toda la red espacio-temporal puede caer en la fase "fácil", lo que implica que DTOE podría ser particularmente poderoso en dimensiones superiores, aunque esta generalización no se prueba explícitamente en este trabajo. Los autores también observan que mientras que DAOE se enmarca como un método MPO con un superoperador de disipación específico, enmarcar DTOE de manera similar sigue siendo un desafío para trabajos futuros.