Computation and properties of the Epstein zeta function with applications to quantum systems

Resumen Técnico: Computación y Propiedades de la Función Zeta de Epstein

Planteamiento del Problema
La función zeta de Epstein, $Z_{\Lambda, \nu}(x, y)$, generaliza la función zeta de Riemann a sumas de redes oscilatorias en dimensiones superiores. Es una herramienta crítica en la física teórica y la química para modelar sistemas de muchos cuerpos clásicos y cuánticos con interacciones de largo alcance, incluyendo constantes de Madelung, energías de Casimir y relaciones de dispersión de ondas de espín. A pesar de su importancia, la falta de un tratamiento riguroso y exhaustivo de su computación y propiedades analíticas ha sido una carencia constante. Los métodos existentes, como la fórmula de Chowla–Selberg o la representación de Crandall, a menudo enfrentan limitaciones: pueden estar restringidos a casos especiales de vectores de red o parámetros, sufrir de inestabilidades numéricas cerca de las singularidades (por ejemplo, argumentos negativos en funciones gamma incompletas) o carecer de una implementación de librería unificada. Además, la continuación holomorfa conjunta y la estructura precisa de las singularidades complejas de la función en todos los parámetros ($\nu, x, y$) no habían sido establecidas rigurosamente, lo que dificultaba el desarrollo de algoritmos numéricos robustos que evitaran la cancelación catastrófica.

Metodología
Los autores abordan estas brechas mediante una combinación de reformulación analítica rigurosa y diseño de algoritmos numéricos de alto rendimiento:

1. Reformulación Analítica:

   • El artículo introduce una representación de tipo Crandall basada en la descomposición del núcleo $|z|^{-\nu}$ en "funciones de Crandall": una función superior $G_\nu$ (de decaimiento superexponencial) y una función inferior $g_\nu$ (suave).
   • Al aplicar la fórmula de sumación de Poisson a la parte inferior, los autores derivan una representación de la función zeta de Epstein como la suma de dos sumas de red de convergencia superexponencial. Esta formulación sirve como base tanto para el algoritmo numérico como para las demostraciones analíticas.
   • Los autores demuestran rigurosamente la continuación holomorfa conjunta de la función zeta de Epstein en el exponente $\nu$ y en los argumentos vectoriales $x, y$ dentro de conos complejos específicos. Caracterizan la estructura de singularidad, identificando polos simples en $\nu = d$ (cuando $y \in \Lambda^*$) y singularidades eliminables en otros casos.
   • Se define una función zeta de Epstein regularizada mediante la sustracción de la singularidad de ley de potencia (derivada de la transformada de Fourier distributiva de $| \cdot |^{-\nu}$). Esta regularización elimina las singularidades sin provocar una cancelación catastrófica, permitiendo la evaluación estable y la precomputación eficiente mediante interpolación.

2. Algoritmo Numérico (EpsteinLib):

   • Los autores implementan una biblioteca en C, EpsteinLib, con enlaces (bindings) para Python, Julia y Mathematica.
   • Preprocesamiento: El algoritmo reescala las redes a volumen unitario y proyecta los vectores en celdas de red elementales para explotar las simetrías de traslación y escala.
   • Casos Especiales: Se proporciona un manejo explícito para las singularidades eliminables (por ejemplo, $\nu \in -2\mathbb{N}_0$) y los polos ($\nu = d$).
   • Truncamiento y Sumación: Se derivan cotas de error rigurosas para el truncamiento de las sumas de decaimiento superexponencial. Para mitigar los errores de redondeo de punto flotante en dimensiones altas, los autores emplean una sumación compensada de Kahan.
   • Funciones Gamma Incompletas: Se desarrolla una implementación personalizada y de alto rendimiento de la función gamma incompleta para manejar argumentos negativos y el rango completo de parámetros con precisión doble completa, superando las limitaciones de las librerías existentes.

Contribuciones Clave

• Primer Tratamiento Analítico Integral: El artículo proporciona la primera prueba de la continuación holomorfa conjunta de la función zeta de Epstein en todos los parámetros y ofrece una caracterización completa de su estructura de singularidades.
• Algoritmo de Convergencia Superexponencial: Se introduce un nuevo algoritmo que computa la función para parámetros reales arbitrarios en cualquier dimensión con cotas de error rigurosas.
• Eliminación de Singularidades: Se presenta un método para descomponer la función en una parte singular y una parte analítica regularizada, lo que permite la evaluación de integrales y el tratamiento de anomalías de tipo Wood sin cancelación catastrófica.
• Implementación de EpsteinLib: La primera librería de alto rendimiento y código abierto para computar la función zeta de Epstein para argumentos reales arbitrarios, que cuenta con evaluación de precisión completa y un extenso benchmarking.

Resultados

• Precisión y Rendimiento: La implementación logra una precisión de doble precisión completa en dimensiones $d=1, 2, 3, 4, 6, 8$ y en todo el rango de parámetros $\nu \in [-12.5, 12.5]$. Las pruebas de rendimiento contra fórmulas analíticas conocidas e implementaciones de precisión arbitraria confirman la precisión del algoritmo.
• Eficiencia: El algoritmo demuestra una convergencia superexponencial. Por ejemplo, las sumas en 2D se computan en menos de 20 microsegundos, y las sumas en 8D en menos de 0.5 segundos en una computadora portátil estándar. El uso de la sumación de Kahan garantiza la estabilidad incluso en dimensiones altas donde los errores de redondeo suelen acumularse.
• Aplicaciones Demostradas:
  • Sistemas de Espín Cuántico: El método se aplicó para computar relaciones de dispersión cuántica en sistemas de espín 3D con interacciones de largo alcance, capturando con éxito comportamientos de escala anómalos ($\omega(k) \sim |k|^{\nu-d}$) para $d < \nu < d+2$.
  • Efecto Casimir: La librería se utilizó para computar energías de Casimir en geometrías 3D, reproduciendo con precisión las leyes de fuerza asintóticas y analizando la minimización de energía en cajas deformadas.

Significado
El artículo sostiene que este trabajo establece la función zeta de Epstein como una función especial estándar para la física computacional y las matemáticas. Al proporcionar una base analítica rigurosa y una implementación robusta y de alto rendimiento, los autores permiten la simulación eficiente de sistemas con interacciones de largo alcance que anteriormente eran computacionalmente prohibitivos o numéricamente inestables. La capacidad de manejar redes generales, parámetros reales arbitrarios y singularidades sin errores de cancelación facilita nuevas investigaciones en química teórica (estabilidad de la materia), física de la materia condensada (diagramas de fase cuánticos, escaleras de diablo) y problemas de dispersión de ondas. Los autores señalan que la librería ya está siendo utilizada por la comunidad de materia condensada cuántica para estudiar fenómenos exóticos, tales como las plataformas de magnetización fraccionaria y la fusión de las escaleras de diablo.