Zeta expansion for long-range interactions under periodic boundary conditions with applications to micromagnetics

Resumen Técnico: Expansión Zeta para Interacciones de Largo Alcance bajo Condiciones de Contorno Periódicas

Planteamiento del Problema
El artículo aborda el desafío computacional de evaluar potenciales de interacción basados en leyes de potencia (específicamente potenciales de Riesz de la forma $|r|^{-\nu}$) para cuerpos homogéneos de $d$ dimensiones que interactúan con un arreglo infinito de copias de $n$ dimensiones bajo condiciones de contorno periódicas (PBC). Este problema es central en la micromagnética, donde la simulación precisa de texturas magnéticas en cables ($n=1$) o películas ($n=2$) requiere el manejo de sumas de red infinitas de interacciones dipolares ($\nu=1$) y sus derivadas.

Las prácticas estándar actuales en el software micromagnético popular (por ejemplo, MuMax3, OOMMF) dependen de la truncación de estas sumas de red infinitas a un número finito de imágenes. Este enfoque introduce errores incontrolados porque los núcleos de interacción decaen algebraicamente en lugar de exponencialmente. Como consecuencia, la periodicidad de las interacciones locales es exacta, pero las interacciones no locales son inexactas. Estas discrepancias no pueden eliminarse mediante el refinamiento de la malla. Aunque existe la sumatoria de Ewald para cargas puntuales, esta no es directamente aplicable a los dominios de volumen finito y a las derivadas de orden superior requeridas en micromagnética.

Metodología: La Expansión Zeta
Los autores proponen un método de "Expansión Zeta" que descompone la suma de red infinita en dos partes:

1. Una Suma Directa Pequeña: Una suma truncada sobre una red de "campo cercano" ($L_{near}$), similar a las implementaciones estándar.
2. Un Término de Corrección: Un término derivado analíticamente que elimina el error de truncación.

La idea teórica central es que el término de corrección puede expresarse como una serie de Taylor convergente que involucra derivadas de una función zeta de conjunto generalizada (una continuación meromorfa de las sumas de red). Para un dominio compacto $\Omega$, el potencial $U^{(m)}(r)$ viene dado por:
$$U^{(m)}(r) = \sum_{z \in L_{near}} S^{(m)}(r+z) + D_\Omega Z^{(m)}_{L_{far}, \nu}(r)$$
donde $S^{(m)}$ es la integral de interacción directa, $Z$ es la función zeta del conjunto para la red de campo lejano, y $D_\Omega$ es un operador diferencial dependiente de la forma definido por los momentos del dominio $\Omega$. Para cuboides alineados con los ejes, este operador toma una forma explícita que involucra derivadas de orden par.

Para evaluar la función zeta y sus derivadas de manera eficiente, los autores utilizan una representación de Crandall. Esta representación divide la función zeta en:

• Una suma en el espacio real que involucra "funciones de Crandall superiores" ($G_\nu$).
• Una suma en el espacio recíproco que involucra "funciones de Crandall inferiores" ($g_\nu$) y funciones de Bessel incompletas ($G_\nu(k, r)$).

Ambas sumas en la representación de Crandall decaen de forma superexponencial, lo que permite su truncación con un error despreciable. Los autores proporcionan algoritmos para calcular estas funciones especiales (incluidas las funciones de Bessel incompletas) con precisión de máquina completa, manejando singularidades y derivadas analíticamente mediante relaciones de recurrencia y fórmulas de reflexión.

Contribuciones Clave

1. Representación Exacta: El artículo deriva una fórmula exacta para las sumas de red infinitas de potenciales de Riesz y sus derivadas, complementando una suma directa con un término de corrección basado en la función zeta.
2. Convergencia Superexponencial: El algoritmo propuesto alcanza la precisión de máquina con un costo computacional comparable a los esquemas de sumatoria truncada estándar, ya que el término de corrección converge superexponentialmente.
3. Aplicabilidad General: El marco soporta exponentes de interacción $\nu$ arbitrarios, dimensiones $d$ y dimensiones periódicas $n$ arbitrarias ($n < d$), y geometrías de origen/objetivo arbitrarias (incluyendo fronteras curvas mediante integrales de momento).
4. Evaluación de Funciones Especiales: Los autores proporcionan un algoritmo robusto y de convergencia superexponencial para evaluar las funciones de Bessel incompletas y sus derivadas, las cuales son críticas para el componente de espacio recíproco del método.
5. Nuevas Correcciones Asintóticas: El método genera nuevas correcciones al límite asintótico del campo de desmagnetización para cuboides de tamaño finito en arreglos periódicos, identificando específicamente un término de corrección que escala como $N^{-7}$ para $N$ grande.

Resultos y Pruebas de Rendimiento (Benchmarks)

• Precisión: El método alcanza consistentemente la precisión de máquina completa (errores relativos $\approx 10^{-16}$) a través de diversas geometrías y órdenes de derivada.
• Comparación con Sumatoria Directa: Para exponentes de interacción grandes ($\nu \ge 10$), la expansión zeta coincide con la sumatoria directa sobre mallas grandes. Para el caso micromagnético ($\nu=1$), los autores demuestran que la truncación estándar (por ejemplo, $N_{cut}=10$) resulta en errores relativos que superan el 15%, mientras que la expansión zeta permanece exacta.
• Factores de Desmagnetización: Los autores tabulan valores de alta precisión para el factor de desmagnetización $D_z$ de cuboides 3D en arreglos periódicos 2D, sirviendo como referencia para futuros solvers.
• Comportamiento Asintótico: El método reproduce con éxito los límites asintóticos conocidos (por ejemplo, $D_z \to 1/3$ para dipolos puntuales, $D_z = 1$ para capas continuas) y proporciona el término de corrección de siguiente orden.

Significado y Reivindicaciones
Los autores afirman que este trabajo resuelve el problema de larga data de los errores de truncación incontrolados en las simulaciones micromagnéticas periódicas. Al aprovechar las funciones zeta generalizadas, el método permite la evaluación de las energías e interacciones de campo con precisión de máquina a un costo computacional no mayor que los esquemas truncados actuales.

El artículo posiciona el método como una "pre-rutina" para los solvers micromagnéticos: los núcleos de interacción pueden precomputarse con alta precisión y luego utilizarse en algoritmos iterativos estándar (por ejemplo, convolución basada en FFT). Los autores enfatizan que la técnica es ampliamente aplicable más allá de la micromagnética, citando aplicaciones potenciales en:

• Ferroeléctricos: Debido a la dualidad electromagnética, el método se aplica a texturas de polarización.
• Dinámica de Espín Atómica: Para modelar interacciones de largo alcance como RKKY en multicapas.
• Dinámica Molecular: Para el manejo de condiciones de contorno periódicas en dimensiones no acotadas.

El trabajo no propone nuevos fenómenos físicos, sino que proporciona una herramienta numérica rigurosa para eliminar artefactos numéricos en modelos físicos existentes, permitiendo así estudios más fiables de sistemas donde las diferencias de energía mínimas (por ejemplo, en sistemas frustrados) son críticas.