FORTVSH: A Fortran Module For Vector Spherical Harmonics Computations

Resumen Técnico: FORTVSH – Un módulo de Fortran para computaciones de Armónicos Esféricos Vectoriales

Planteamiento del Problema
Los armónicos esféricos vectoriales (VSH) son fundamentales para representar campos vectoriales arbitrarios en sistemas de coordenadas curvilíneas dentro de la física matemática, particularmente para el análisis espectral en magnetohidrodinámica (MHD). Un desafío computacional primario en las simulaciones de MHD, como aquellas que modelan sistemas geofísicos y astrofísicos, es la descomposición eficiente de los campos magnéticos en componentes poloidales y toroidales. Si bien los paquetes de código abierto existentes como SHTOOLS y SPHEREPACK proporcionan una infraestructura robusta para armónicos esféricos escalares y transformaciones vectoriales específicas, carecen de soporte directo para la representación completa del campo magnético poloidal-toroidal requerida para las ecuaciones de inducción 3D de la MHD. SHTOOLS está limitado a campos libres de rotación (carece de componentes toroidales), y SPHEREPACK está restringido a campos tangenciales 2D en una única capa, fallando en capturar el acoplamiento radial-tangencial inherente a los campos poloidales. Existe la necesidad de una biblioteca ligera y de alto rendimiento que preserve las convenciones algebraicas de la física aplicada y facilite la expansión espectral de campos vectoriales 3D sin la sobrecarga de abstracciones de alto nivel o dependencias externas.

Metodología
El artículo presenta FORTVSH, un módulo de Fortran de código abierto diseñado para evaluar diversas formas de VSH y funciones matemáticas relacionadas. La metodología se fundamenta en el marco teórico de los armónicos esféricos tensoriales (TSH) y el acoplamiento específico de los momentos angulares orbital ($L$) y de espín ($S$).

1. Base Teórica: El código implementa formas de VSH derivadas de la expansión de TSH donde el rango tensorial $S=1$. Utiliza dos bases primarias:
   • La base de funciones propias $L^2$ ($Y^L_{JM}$).
   • La base alternativa ($Y^{(\lambda)}_{JM}$), que no es una función propia de $L^2$ pero forma una base ortonormal completa adecuada para la descomposición poloidal-toroidal. Esta base está linealmente relacionada con las funciones propias de $L^2$ mediante coeficientes de Clebsch-Gordan.
2. Implementación Algebraica: El módulo incorpora identidades para las operaciones de rotacional (curl) sobre modos VSH acoplados con funciones radiales, permitiendo la transformación de operadores diferenciales (crucial para la ecuación de inducción magnética $\partial_t \mathbf{B} = -c\nabla \times \mathbf{E}$) en relaciones algebraicas entre modos espectrales.
3. Enfoque Numérico:
   • Batch vs. Modo Único: La biblioteca proporciona tanto funciones de modo único como rutinas de procesamiento por lotes ("batch") (con el sufijo _ALL) que computan rejillas completas de modos simultáneamente.
   • Estabilidad: La implementación utiliza relaciones de recurrencia normalizadas (Holmes y Featherstone, 2002) para los polinomios de Legendre asociados, asegurando la estabilidad numérica hasta grados elevados ($\ell_{max} = 2000$).
   • Validación: El código se valida contra formas analíticas independientes, valores de referencia de mpmath de Python y soluciones de forma cerrada derivadas de los coeficientes de Geppert-Wiebicke.

Contribuciones Clave

• Paquete FORTVSH: Una biblioteca de Fortran ligera y sin dependencias (sin capa de Python, backend de FFT o requisitos de BLAS/LAPACK) que proporciona un catálogo de rutinas numéricas para evaluar formas de VSH, incluyendo componentes poloidales y toroidales.
• Descomposición Poloidal-Toroidal: El paquete soporta explícitamente la representación de campos magnéticos descompuestos en potenciales poloidales y toroidales, una capacidad ausente en las principales bibliotecas existentes.
• Identidades Completas: El manuscrito establece e implementa un conjunto completo de identidades analíticas, incluyendo productos internos de modos VSH, operaciones de rotacional y relaciones de rotación entre diferentes bases VSH, utilizando símbolos de Clebsch-Gordan y Wigner.
• Optimización del Rendimiento: La introducción de rutinas de procesamiento por lotes que computan campos vectoriales sobre una rejilla completa simultáneamente, ofreciendo ventajas computacionales significativas sobre los bucles de modo único.

Resultados
El artículo presenta una extensa validación numérica y resultados de benchmarking:

• Precisión: Las rutinas de FORTVSH reproducen las formas analíticas y los valores de referencia con una precisión de máquina (errores absolutos típicamente por debajo de $10^{-15}$). Esto incluye la verificación de la ortonormalidad, la ortogonalidad poloidal-toroidal y la inversión de las relaciones de rotación de VSH.
• Estabilidad: Las rutinas de Legendre no normalizadas son estables hasta $\ell_{max} = 200$, mientras que las rutinas normalizadas permanecen estables hasta $\ell_{max} = 2000$.
• Rendimiento: Las pruebas de tiempo en un AMD Ryzen 7 5700G demuestran que las rutinas de procesamiento por lotes superan a los bucles de modo único por un orden de magnitud dentro del rango estable superpuesto. Además, las rutinas de procesamiento por lotes extienden el rango de parámetros estables por un orden de magnitud adicional en comparación con las evaluaciones de modo único.
• Casos de Aplicación: El artículo demuestra la utilidad de la biblioteca a través de tres casos de uso:
  1. Síntesis de Dipolo Puro: Reconstrucción fiel de un campo dipolar de vacío con contaminación espectral insignificante.
  2. Acoplamiento de Frontera Magnetostática: Solución precisa de las condiciones de contorno en la interfaz de una esfera uniformemente magnetizada, satisfaciendo la continuidad del campo radial y las condiciones de salto para el campo tangencial.
  3. Descomposición Espectral: Descomposición exitosa de campos vectoriales no axisimétricos arbitrarios, incluyendo un caso de prueba con un bulto gaussiano toroidal, mostrando una convergencia espectral característica.

Significancia
El artículo posiciona a FORTVSH como una base numéricamente fiable para la representación espectral de campos vectoriales en sistemas geofísicos y astrofísicos. Su importancia radica en su capacidad para facilitar simulaciones numéricas de evolución temporal de sistemas MHD en geometría esférica, al proporcionar las herramientas algebraicas necesarias para descomponer los campos magnéticos en componentes poloidales y toroidales. Al ofrecer un enfoque de "primacía de la física" que preserva las convenciones algebraicas de la literatura de física aplicada, FORTVSH está diseñado para ser embebido directamente en esquemas de avance temporal creados específicamente para sistemas MHD radialmente estratificados, tales como aquellos que modelan interacciones núcleo-corteza de estrellas de neutrones, dinamos planetarios y magnetosferas astrofísicas. El autor señala que la versión actual (v1.0.0) se centra en las operaciones de transformación directa, con la posibilidad de capacidades de transformación inversa previstas para futuras versiones.