Deferred Cyclotomic Representation for Stable and Exact Evaluation of q-Hypergeometric Series

1. El Problema: Evaluación de Series q-Hipergeométricas

Las series q-hipergeométricas finitas son fundamentales en teoría de representación, topología cuántica y física matemática (ej. símbolos 6j cuánticos, invariantes de nudos, modelos de gravedad cuántica). Sin embargo, su evaluación computacional enfrenta dos obstáculos críticos que limitan su escalabilidad y precisión:

• Inestabilidad Numérica (Cancelación Catastrófica): En aritmética de punto flotante, estas series a menudo involucran sumas alternas de términos exponencialmente grandes que se cancelan para producir un resultado pequeño. Esto genera un número de condición ($\kappa$) extremadamente alto, provocando una pérdida masiva de precisión (cancelación catastrófica) incluso en aritmética de doble precisión.
• Hinchazón de Expresiones (Expression Swell): En la evaluación simbólica exacta (álgebra computacional), las series se expresan tradicionalmente como funciones racionales o polinomios densos en el parámetro de deformación $q$. Al expandir factoriales cuánticos, el tamaño de las expresiones intermedias crece rápidamente, superando la capacidad de memoria y tiempo de cómputo debido a la necesidad de cálculos repetidos de MCD (máximo común divisor) y simplificación de polinomios.

El artículo argumenta que estos problemas no son meramente algorítmicos, sino que surgen de una mismatch de representación: las estructuras multiplicativas intrínsecas de los enteros y factoriales cuánticos se pierden al representarlas como polinomios densos aditivos.

2. Metodología: Representación Ciclotómica Diferida (DCR)

El autor introduce la Representación Ciclotómica Diferida (DCR), un marco que separa la estructura algebraica de la evaluación numérica. La metodología se basa en los siguientes pilares:

A. Factorización Ciclotómica

En lugar de tratar los enteros cuánticos $[n]_q$ y factoriales $[n]_q!$ como polinomios en $q$, se descomponen en sus factores irreducibles sobre los racionales: los polinomios ciclotómicos $\Phi_d(q^2)$.
La factorización clave es:
$$[n]_q = q^{1-n} \prod_{d|n, d>1} \Phi_d(q^2)$$
$$[n]_q! = q^{\frac{n(1-n)}{2}} \prod_{d=2}^{n} \Phi_d(q^2)^{\lfloor n/d \rfloor}$$

B. Base de Exponentes Enteros Esparsos

Cada término de la serie se representa como un monomio ciclotómico definido por un vector de exponentes enteros esparsos sobre la base $\{q, \Phi_2(q^2), \Phi_3(q^2), \dots\}$.

• Las operaciones multiplicativas (producto, cociente) se convierten en suma y resta de vectores de enteros.
• Las cancelaciones entre numerador y denominador se resuelven exactamente mediante la resta de exponentes, antes de cualquier evaluación numérica.

C. Compilación Diferida y Proyección

El proceso se divide en dos etapas:

1. Etapa de Compilación (Algebraica): Se construye un objeto combinatorio fijo ($S_{DCR}$) que codifica toda la estructura multiplicativa de la serie como vectores de enteros. Esta etapa es independiente del valor de $q$ y no genera expresiones densas.
2. Etapa de Proyección (Evaluación): La evaluación en un campo objetivo (enteros, punto flotante, raíces de la unidad) se realiza aplicando un homomorfismo de anillo $\Pi_T$ que mapea los generadores ciclotómicos a sus valores numéricos específicos.

3. Contribuciones Clave

1. Representación Ciclotómica: Demostración de que las amplitudes q-hipergeométricas admiten una representación en una base de exponentes ciclotómicos esparsos, revelando su estructura multiplicativa oculta.
2. Evaluación Diferida: Introducción del marco DCR, que separa la compilación combinatoria (una sola vez) de la evaluación posterior mediante proyección.
3. Reducción de Complejidad: El DCR escala casi linealmente con la longitud de la suma, evitando el crecimiento exponencial de expresiones típico de las representaciones polinómicas.
4. Estabilidad Numérica: Al realizar cancelaciones algebraicas exactas en el nivel de exponentes antes de la evaluación numérica, el DCR actúa como un precondicionador algebraico, reduciendo drásticamente el rango dinámico de los términos intermedios y mitigando la pérdida de precisión.
5. Unificación Estructural: Proporciona una perspectiva unificada donde la admisibilidad en raíces de la unidad, la deformación $q$ y el límite clásico ($q \to 1$) emergen como propiedades intrínsecas de un único objeto combinatorio.

4. Resultados Empíricos y de Rendimiento

El autor valida el marco utilizando el símbolo cuántico 6j como caso de prueba (implementado en el paquete Julia QRecoupling.jl):

• Uso de Memoria:
  • La construcción del DCR escala linealmente ($O(j)$) con el espín $j$.
  • La proyección exacta escala cúbicamente ($O(j^3)$).
  • En contraste, la evaluación "eager" (ansiosa) con CAS (Sistemas de Álgebra Computacional) sufre de hinchazón de expresiones, excediendo 50 GB de memoria para $j=120$, mientras que el DCR se mantiene en el orden de kilobytes.
• Estabilidad Numérica:
  • En aritmética de doble precisión (Float64), el método tradicional (Log-Sum-Exp) pierde precisión y signo correcto alrededor de $j \approx 90$.
  • El DCR mantiene la precisión y el signo correcto hasta valores mucho mayores, extendiendo el rango de evaluación fiable.
  • El factor de amplificación de error ($\gamma_{rep}$) se reduce en varios órdenes de magnitud (ej. de ~9500 a ~1000 en $j=500$) gracias a la compresión del rango dinámico.
• Latencia y Amortización:
  • Aunque la construcción inicial tiene un costo, la evaluación posterior es extremadamente rápida (microsegundos) porque solo requiere la proyección.
  • Esto permite un barrido eficiente de parámetros continuos sin recalcular la estructura combinatoria.

5. Significado e Impacto

El trabajo tiene implicaciones profundas tanto computacionales como teóricas:

• Computacional: Hace viables cálculos de topología cuántica y gravedad cuántica de bucles que antes eran prohibitivos debido a la inestabilidad numérica o el agotamiento de memoria. Permite evaluar invariantes de estados (como la función de partición de Turaev-Viro) en redes trianguladas grandes mediante la compilación de bloques locales reutilizables.
• Teórico: Cambia la visión de las amplitudes cuánticas. En lugar de ser funciones que deben re-calcularse para cada $q$, se ven como evaluaciones de un objeto combinatorio fijo.
  • La admisibilidad en raíces de la unidad se interpreta como la anulación de modos ciclotómicos específicos ($\Phi_h(q^2) = 0$).
  • El límite clásico ($q \to 1$) se reduce a la factorización prima de factoriales clásicos.
• Generalidad: El principio de alinear la representación con la estructura multiplicativa intrínseca sugiere que enfoques similares podrían aplicarse a otras funciones especiales y amplitudes cuánticas, mejorando simultáneamente la computabilidad exacta, la estabilidad numérica y la eficiencia algorítmica.

En resumen, el DCR transforma un problema de evaluación numérica inestable y costoso en un problema de manipulación de enteros esparsos y proyección algebraica, resolviendo los cuellos de botella fundamentales en la computación de series q-hipergeométricas.