Deficiency Identity for Zeta Values and Optimal Approximation

Resumen Técnico: Una Identidad de Déficit para Valores Zeta y Aproximación Óptima

1. Planteamiento del Problema

La función Zeta de Riemann, definida como $\zeta(s) = \sum_{n=1}^{\infty} n^{-s}$ para $\Re(s) > 1$, es un objeto central en la teoría de números analítica. Mientras que los valores en enteros pares admiten fórmulas cerradas explícitas (en términos de números de Bernoulli), no existen expresiones cerradas comparables para los valores en enteros impares ($\zeta(3), \zeta(5), \dots$).

Desde una perspectiva computacional, la aproximación estándar mediante truncamiento de la serie de Dirichlet ($\sum_{k=1}^n k^{-q}$) tiene una tasa de convergencia lenta, del orden de $O(n^{-(q-1)})$. Los métodos clásicos como la suma de Euler-Maclaurin mejoran esta tasa mediante correcciones explícitas, pero requieren un mayor complejo simbólico y el cálculo de derivadas o números de Bernoulli. El objetivo de este trabajo es proponer un nuevo marco de aproximación basado en una identidad estructural exacta que evite la expansión explícita de la cola de la serie.

2. Metodología: El Marco de Déficit No Lineal

El autor introduce una representación basada en la no linealidad y la convexidad. La idea central es que la acumulación global de una masa antes de aplicar una transformación convexa es más eficiente que aplicar la transformación localmente y luego sumar.

Definiciones Clave:

• Secuencia Base: Se define una suma parcial $S_n^p = \sum_{k=1}^n k^{-p}$.
• Transformación No Lineal: Se considera la transformación convexa $x \mapsto x^{q/p}$ con $q > p$.
• Funcional de Déficit ($D_{p,q}^n$): Se define la diferencia entre la transformación global de la suma parcial y la suma de las transformaciones locales:
  $$D_{p,q}^n := (S_n^p)^{q/p} - \sum_{k=1}^n k^{-q}$$
  Debido a la convexidad, $(S_n^p)^{q/p} \ge \sum_{k=1}^n k^{-q}$, por lo que el déficit es no negativo.

Identidad Exacta:
Al tomar el límite $n \to \infty$, se establece la siguiente identidad exacta para los valores zeta:
$$\zeta(q) = (\zeta(p))^{q/p} - D_{p,q}^{\infty}$$
Donde $D_{p,q}^{\infty}$ es el déficit acumulado total. Esta ecuación interpreta $\zeta(q)$ como la "masa base agregada no linealmente" menos la "pérdida debida a la fragmentación".

3. Esquema de Aproximación y Corrección de Sesgo

El método propone aproximar $\zeta(q)$ utilizando el déficit truncado a nivel $n$:
$$A_{p,q}^n = (\zeta(p))^{q/p} - D_{p,q}^n$$
Sin embargo, el estimador básico $A_{p,q}^n$ presenta un sesgo de primer orden de $O(n^{-1})$.

Estimador Corregido:
Para eliminar este sesgo, el autor deriva un estimador corregido $B_{p,q}^n$ que resta explícitamente el término de error dominante basado en la cola de la secuencia base $t_n = \zeta(p) - S_n^p$:
$$B_{p,q}^n = A_{p,q}^n - \frac{q}{p} (\zeta(p))^{q/p - 1} t_n$$
Este estimador elimina el término $O(n^{-1})$, mejorando la tasa de convergencia.

4. Contribuciones Clave y Resultados Teóricos

• Descomposición Estructural Exacta: Se proporciona una identidad rigurosa que descompone $\zeta(q)$ en una parte no lineal y un déficit acumulativo, sin necesidad de expansiones asintóticas previas.
• Tasas de Convergencia Óptimas:
  • Para el estimador corregido general, la tasa de error es:
    $$B_{p,q}^n - \zeta(q) = O(n^{-\min(2p-2, \, q-1)})$$
  • Optimización: Para valores impares $q = 2m+1$, la elección óptima del parámetro base es $p = q-1$. Esto resulta en una tasa de convergencia:
    $$B_{q-1,q}^n - \zeta(q) = O(n^{-(q-1)})$$
    Esta tasa coincide con la tasa de decaimiento natural del error de truncamiento, pero se logra manteniendo la estructura del déficit.
• Relación con Euler-Maclaurin: El artículo demuestra que el método actúa como una expansión de Euler-Maclaurin implícita. Las correcciones asintóticas surgen automáticamente de la expansión de Taylor de la transformación no lineal $(S_n^p)^{q/p}$, evitando el uso explícito de números de Bernoulli o derivadas de orden superior.
• Extensión a Funciones Zeta Espectrales: El marco se generaliza a funciones zeta asociadas a secuencias de autovalores $\{\lambda_k\}$ de operadores autoadjuntos positivos. Se demuestra que el comportamiento de convergencia depende de la interacción entre el exponente base $p$ y la tasa de crecimiento espectral.

5. Validación Numérica

Los experimentos numéricos realizados para $\zeta(3), \zeta(5)$ y $\zeta(7)$ confirman las predicciones teóricas:

• El estimador corregido con $p=2$ logra consistentemente una convergencia de $O(n^{-2})$.
• La elección optimizada $p = q-1$ acelera drásticamente la convergencia, alineándose con la tasa predicha $O(n^{-(q-1)})$.
• Aunque el método de Euler-Maclaurin sigue siendo superior asintóticamente cuando se permiten correcciones explícitas de alto orden, el método de déficit ofrece una ventaja conceptual y una implementación unificada más simple.
• Las pruebas en contextos espectrales (con autovalores $\lambda_k = k^\alpha$) validan la adaptabilidad del método a diferentes tasas de crecimiento.

6. Significado y Conclusión

La principal contribución de este trabajo es estructural y conceptual, más que puramente numérica en términos de velocidad de convergencia extrema.

1. Nueva Perspectiva: Ofrece una interpretación geométrica y física (eficiencia de volumen bajo fragmentación) para los valores de la función Zeta.
2. Unificación: Proporciona un mecanismo de corrección unificado para todos $q > p$ que genera automáticamente términos de corrección asintótica.
3. Limitaciones: El método no produce nuevas fórmulas cerradas para valores impares de Zeta ni supera el orden asintótico máximo alcanzable por las expansiones de Euler-Maclaurin truncadas adecuadamente.

En resumen, Mohammady presenta un marco elegante que transforma el problema de aproximar series de Dirichlet en un problema de gestión de "déficit" mediante transformaciones no lineales, logrando convergencias de alto orden de manera implícita y extendible a contextos espectrales generales.