Graph lattice sums and graph zeta functions for long-range interacting quantum lattice models

Resumen Técnico: Sumas de Redes de Grafos y Funciones Zeta de Grafos para Modelos de Red Cuánticos de Interacción de Largo Alcance

Planteamiento del Problema
La simulación de modelos de red cuánticos con brecha (gapped) y de interacciones de largo alcance (que decaen como $|x|^{-\nu}$) se ve obstaculizada por el crecimiento exponencial de la dimensión del espacio de Hilbert con el tamaño del sistema. Aunque las expansiones de clústeres vinculados de alto orden ofrecen una ruta rigurosa para resolver el problema de autovalores para sistemas infinitos, estas requieren el cálculo de sumas de redes de alta dimensión asociadas a estructuras de grafos. Actualmente, estas sumas se evalúan mediante métodos de Monte Carlo, los cuales sufren de errores estocásticos y tiempos de ejecución a escala de clúster (requiriendo a menudo miles de horas-núcleo para un solo punto de momento). Además, resolver la relación de dispersión completa requiere recomputar estas sumas para cada vector de onda $k$, lo que agrava el costo computacional. Los desafíos principales son:

1. Complejidad Exponencial: La discretización directa escala exponencialmente con el número de nodos del grafo.
2. Convergencia Lenta: Para núcleos de ley de potencia con exponentes $\nu$ cercanos a la dimensión de la red $d$, la convergencia es arbitrariamente lenta.
3. Resolución de la Rejilla de Momento: Los métodos estándar no pueden recuperar eficientemente toda la zona de Brillouin.

Metodología
Los autores proponen un marco determinista de alto rendimiento que reemplaza el muestreo estocástico por una combinación de teoría de grafos, funciones zeta generalizadas y redes de tensores. El núcleo del método consiste en factorizar las sumas de redes de grafos en componentes más pequeños y computables basados en la topología del grafo.

1. Factorización por Bloques:
   La suma de la red del grafo $Z_{\Lambda, G}(k)$ se descompone sobre los bloques del grafo (subgrafos máximos 2-conectados) y los puntos de corte (cutvertices). Utilizando la invarianza de traslación, la suma se factoriza en un producto de sumas sobre estos bloques. Solo los bloques de la "espina" (aquellos en el camino entre los terminales) dependen del momento $k$; los bloques de adherencia se evalúan en $k=0$. Esto reduce la complejidad de exponencial en el número total de nodos a exponencial en el tamaño del bloque más grande.

2. Evaluación Analítica de Bloques Básicos:

   • Puentes (Bridges): Para puentes (dos nodos conectados por aristas), la suma es la transformada de Fourier de la red del núcleo de interacción. Para núcleos de ley de potencia, esta es la función zeta de Epstein, la cual se computa analíticamente o mediante librerías de alto rendimiento (EpsteinLib).
   • Ciclos: Para ciclos sin dependencia de momento, la suma se representa como una integral sobre la zona de Brillouin del producto de las transformadas de Fourier, computable mediante transformadas de Duffy o cuadratura de Gauss.

3. Álgebra Semi-analítica para Bloques Serie-Paralelo:
   Para bloques con un ancho de árbol (treewidth) $tw \le 2$ (grafos serie-paralelo), los autores construyen la suma de la red del grafo mediante la multiplicación puntual y la convolución de objetos base (puentes).

   • Para manejar las colas de ley de potencia y las singularidades en $k=0$, introducen una representación semi-analítica que consiste en una serie de Fourier de decaimiento rápido (para la parte regular) y una suma finita de funciones zeta de Epstein (para las partes singulares).
   • Los Teoremas 6.2 y 6.3 demuestran que esta álgebra es cerrada bajo multiplicación y convolución. Las singularidades se rastrean analíticamente, mientras que las partes regulares se manejan mediante FFT.
   • Una estrategia de compresión descarta los términos de Epstein con exponentes altos (que se vuelven suaves) hacia la serie de Fourier, evitando la explosión de términos y evitando errores de cancelación numérica.

4. Eliminación de Cubetas por Redes de Tensores para Grafos Altamente Conectados:
   Para bloques con $tw \ge 3$ (donde la reducción serie-paralelo se detiene), el método emplea la eliminación de cubetas (bucket elimination) mediante redes de tensores.

   • La suma de la red se representa como una red de tensores donde los núcleos de las aristas actúan como operadores de convolución.
   • Explotando la invarianza de traslación, los autores reducen la escala numérica. El costo escala polinomialmente con el tamaño de la rejilla de momento $N$ como $O(N^{tw})$ (o $O(N^{tw-1} \log N)$ con FFT), en lugar de exponencialmente con el número de nodos.
   • Crucialmente, esto permite la evaluación de toda la rejilla de momento al costo de una sola evaluación de $k$.

Contribuciones Clave

• Resolución de Sumas de Redes de Grafos: El artículo proporciona el primer método para computar todas las sumas de redes de grafos requeridas (funciones zeta de grafos) para núcleos de ley de potencia de forma determinista, eliminando la necesidad de muestreo de Monte Carlo.
• Reducción de Complejidad: El método reduce el tiempo de evaluación de las expansiones de clústeres vinculados de vanguardia de decenas de miles de horas-núcleo a minutos en un solo núcleo.
• Resolución Completa del Momento: Mediante el uso de FFT y redes de tensores, se recupera la relación de dispersión completa al costo de una sola evaluación de un punto $k$.
• Librería de Zeta de Grafos (GZL): Los autores proporcionan una implementación de alto rendimiento y código abierto (GZL) que reproduce todos los resultados en hardware estándar.

Resultados

• Benchmarking: El método fue evaluado frente a soluciones analíticas para grafos simples y frente a datos de Monte Carlo publicados para el modelo de Ising de campo transversal de largo alcance (LRTFIM) en redes 1D, 2D y 3D.
• Acuerdo: Los resultados muestran un acuerdo total con los datos de Monte Carlo dentro de los errores estadísticos de las referencias. La desviación máxima a través de 998 coeficientes de serie fue de $3.1\sigma_{MC}$, con una mediana de $0.12\sigma_{MC}$.
• Rendimiento:
  • Una relación de dispersión 3D en una rejilla de 4096 puntos de momento se computó en aproximadamente 10 minutos en un solo núcleo.
  • En contraste, reproducir los mismos datos mediante Monte Carlo requirió aproximadamente $10^4$–$10^5$ horas-núcleo en un clúster.
• Precisión: El método logra una precisión numérica controlada incluso para exponentes de interacción $\nu$ cercanos a la dimensión del sistema $d$, donde la suma directa falla.

Significancia y Reivindicaciones
El artículo afirma resolver un cuello de botella fundamental en la simulación de sistemas cuánticos de interacción de largo alcance. Al hacer que las sumas de redes de alta dimensión sean computables con precisión controlada y escala polinomial, el método permite:

• Acceso Riguroso a Sistemas Infinitos: Permite la extracción de magnitudes físicas (energía del estado fundamental, relaciones de dispersión) para redes infinitas mediante la teoría de perturbaciones de alto orden sin errores estocásticos.
• Relevancia Experimental: El artículo complementario [22] utiliza este marco para reproducir observaciones experimentales en materiales cuánticos, específicamente la dispersión del imán de Ising de red triangular $KTmSe_2$.
• Dispersión Anómala: El método resuelve con éxito características no analíticas en la relación de dispersión (escalando como $|k|^\sigma$) que son huellas de interacciones de largo alcance, las cuales son inaccesibles para los enfoques estándar de Monte Carlo debido a la incapacidad de resolver eficientamente toda la rejilla de momento.

Los autores posicionan este trabajo como una base matemática que, combinada con la teoría de perturbaciones de clústeres vinculados, crea un marco integral para estudiar modelos de red cuánticos con brecha (especialmente en sistemas frustrados o de interacción de largo alcance) donde otros métodos como DMRG o Monte Carlo cuántico son ineficaces.