Resumen Técnico: Aproximaciones Simples y Precisas a la Función Zeta de Riemann

Planteamiento del Problema

El cálculo de la función zeta de Riemann, $\zeta(s)$, para valores grandes de la parte imaginaria, $t = \text{Im}(s)$, presenta desafíos significativos. Si bien el método de suma de Euler-Maclaurin ofrece cotas de error rigurosas, su complejidad computacional escala como $O(|s|)$, lo que lo hace poco práctico para $t$ grandes. El enfoque estándar para $t$ grandes es la fórmula de Riemann-Siegel, que reduce la complejidad a $O(t^{1/2})$ mediante el uso de una suma principal y una expansión asintótica para el término de resto. Sin embargo, las correcciones de orden superior en la fórmula de Riemann-Siegel involucran derivadas complejas de funciones especiales (específicamente la función $F(z)$), las cuales son difíciles de evaluar numéricamente debido a singularidades y fórmulas complicadas. Además, los algoritmos existentes de alta precisión, como los de Galway, a menudo requieren el ajuste de múltiples parámetros para lograr una precisión óptima.

El artículo aborda la necesidad de aproximaciones que sean tanto computacionalmente eficientes (requiriendo pocos términos) como altamente precisas en franjas verticales arbitrarias, incluida la franja crítica, sin las dificultades de implementación asociadas con las correcciones de orden superior de Riemann-Siegel.

Metodología

El autor desarrolla una familia de aproximaciones, denotadas $\zeta_p(s)$, que combinan la suma principal de la fórmula de Riemann-Siegel con una nueva aproximación del término de resto. La metodología central se basa en los siguientes pasos:

1. Representación Integral: El término de resto de la fórmula exacta de Riemann-Siegel se expresa como una integral $I_M(s)$ que involucra una medida de valor complejo $\eta(dx) = \frac{e^{-\pi x^2}}{\cosh(\pi \theta x)} dx$, donde $\theta = e^{-i\pi/4}$.
2. Cuantificación Gaussiana: En lugar de utilizar expansiones en series asintóticas, el autor aproxima esta integral mediante una regla de cuantificación gaussiana. El objetivo es encontrar un conjunto de nodos $\{x_{p,j}\}$ y pesos $\{\omega_{p,j}\}$ tales que la cuantificación sea exacta para un conjunto específico de $4p+2$ funciones de prueba.
3. Derivación de Coeficientes: Los coeficientes $\{\omega_{p,j}\}$ y $\{\lambda_{p,j}\}$ (donde $\lambda_{p,j} = \theta x_{p,j}$) se derivan resolviendo un sistema de ecuaciones lineales. Este sistema se construye exigiendo que la cuantificación sea exacta para funciones de la forma $f_k(x) = e^{-\pi x^2 + 2\pi \theta y_{p,k} x}$, donde $y_{p,k}$ son puntos equidistantes en $[-1, 1]$.
   • El problema se transforma en encontrar los nodos y pesos de una cuantificación gaussiana para un funcional lineal $L$ definido por los momentos de la integral de Mordell $H(y)$.
   • Estos nodos y pesos se calculan utilizando polinomios ortogonales generados mediante una relación de recurrencia de tres términos.
4. Fórmula de Aproximación: La aproximación final se define como:
   $$\zeta_p(s) = \sum_{n=1}^{N_t} n^{-s} + \chi(s) \sum_{n=1}^{N_t} n^{s-1} - \frac{(-1)^{N_t}}{2} \left[ I_{N_t+1/2, p}(s) + \chi(s) I_{N_t+1/2, p}(1-s) \right]$$
   donde $I_{M,p}(s)$ es la aproximación por cuantificación del término integral, que involucra únicamente funciones elementales (exponenciales, logaritmos) y coeficientes complejos precalculados.

Contribuciones Clave

• Implementación Simplificada: Las aproximaciones propuestas dependen de un único parámetro entero $p$ y de un conjunto de coeficientes precalculados $\{\omega_{p,j}\}$ y $\{\lambda_{p,j}\}$. A diferencia de la fórmula de Riemann-Siegel, que requiere evaluar derivadas complejas de funciones especiales para términos de orden superior, $\zeta_p(s)$ requiere únicamente funciones elementales y la función Gamma.
• Coeficientes Precalculados: El autor proporciona coeficientes precalculados para $1 \le p \le 30$ y valores seleccionados hasta $p=150$. Se demuestra que estos coeficientes tienen magnitud moderada y exhiben una rápida decadencia, facilitando una evaluación numérica estable.
• Aproximación de Derivadas: El artículo extiende la metodología para aproximar la derivada $\zeta'(s)$ diferenciando la fórmula de aproximación con respecto a $s$ (evaluada en $N_t$ fijo), evitando las discontinuidades asociadas con la diferenciación de la función parte entera en $N_t$.
• Algoritmo Híbrido: Se propone una estrategia para calcular $\zeta(s)$ en todo el plano complejo cambiando entre suma directa, suma de Euler-Maclaurin y la aproximación $\zeta_p(s)$ según la región de $s$.

Resultados

El artículo presenta evidencia numérica extensa sobre la precisión de $\zeta_p(s)$:

• Precisión en la Franja Crítica: Los experimentos numéricos en la franja crítica ($0 \le \sigma \le 1$) demuestran que el error $|\zeta_p(s) - \zeta(s)|$ disminuye significativamente a medida que aumenta $p$.
  • Para $p=10$, los errores son $< 10^{-15}$ para $t > 250$.
  • Para $p=50$, los errores son $< 10^{-100}$ para $t > 4000$.
  • Para $p=150$, los errores son $< 10^{-300}$ para $t > 6900$.
• Patrones de Error: El error exhibe un patrón distintivo donde se anula en los extremos de los intervalos $[t_n, t_{n+1}]$ (donde $t_n = 2\pi n^2$) y en $4p$ puntos equidistantes dentro de estos intervalos. Este comportamiento se explica analíticamente por el hecho de que el término de error dominante es proporcional a $H_p(B(t)) - H(B(t))$, que es cero en los nodos de cuantificación.
• Comparación con Riemann-Siegel: Para $t$ moderado (por ejemplo, $10^2 \le t \le 10^4$), se demuestra que la aproximación $\zeta_7(s)$ es más precisa que la aproximación de Riemann-Siegel con $K=5$ ($\zeta^{RS}_5$). Aunque $\zeta^{RS}_5$ tiene una tasa de decaimiento asintótico más rápida ($O(t^{-13/4})$) en comparación con $\zeta_7$ ($O(t^{-1/4})$), $\zeta_7$ requiere significativamente menos coeficientes precalculados (15 frente a ~100–170 para alta precisión) y es más fácil de implementar.
• Límites de Precisión: En implementaciones de doble precisión, los errores de redondeo dominan el error de aproximación para $t > 200$. Sin embargo, en precisión cuádruple, las aproximaciones mantienen alta precisión hasta valores de $t$ mucho mayores (por ejemplo, $t > 10^4$ con errores $< 10^{-29}$).
• Rendimiento: Se informa que la implementación en MATLAB del algoritmo híbrido es cientos a miles de veces más rápida que la función zeta integrada de MATLAB.

Significado y Afirmaciones

El artículo afirma que estas aproximaciones permiten el cálculo de alta precisión de $\zeta(s)$ y $\zeta'(s)$ en la franja crítica y otras franjas verticales para $t$ grandes. El significado radica en la compensación entre la complejidad de implementación y la precisión: el método ofrece una alternativa "simple" a la fórmula de Riemann-Siegel que evita la carga computacional de los términos de corrección de alto orden mientras mantiene una precisión competitiva o superior para un amplio rango de $t$.

El autor señala que, aunque no se proporcionan cotas de error rigurosas para $\zeta_p(s)$ en este trabajo, la evidencia numérica es extensa. El método se presenta como una herramienta práctica para cálculos de alta precisión, los cuales son cruciales para descubrimientos matemáticos como la refutación de la conjetura de Mertens. El artículo concluye proporcionando implementaciones de código abierto en MATLAB, Python y Fortran, facilitando el uso inmediato por parte de la comunidad investigadora.