Computation and properties of the Epstein zeta function with applications to quantum systems

Résumé technique : Calcul et propriétés de la fonction zêta d'Epstein

Énoncé du problème
La fonction zêta d'Epstein, $Z_{\Lambda, \nu}(x, y)$, généralise la fonction zêta de Riemann aux sommes de réseaux oscillatoires dans des dimensions supérieures. C'est un outil critique en physique théorique et en chimie pour modéliser les systèmes à corps multiples classiques et quantiques à interactions à longue portée, incluant les constantes de Madelung, les énergies de Casimir et les relations de dispersion des ondes de spin. Malgré son importance, un traitement rigoureux et complet de son calcul et de ses propriétés analytiques faisait défaut. Les méthodes existantes, telles que la formule de Chowla–Selberg ou la représentation de Crandall, se heurtent souvent à des limitations : elles peuvent être restreintes à des cas particuliers de vecteurs de réseau ou de paramètres, souffrir d'instabilités numériques près des singularités (par exemple, des arguments négatifs dans les fonctions gamma incomplètes), ou manquer d'une implémentation de bibliothèque unifiée. De plus, la continuation holomorphe conjointe et la structure précise des singularités complexes de la fonction dans tous les paramètres ($\nu, x, y$) n'avaient pas été rigoureusement établies, entravant le développement d'algorithmes numériques robustes capables d'éviter la annulation catastrophique.

Méthodologie
Les auteurs comblent ces lacunes par une combinaison de reformulation analytique rigoureuse et de conception d'algorithmes numériques de haute performance :

1. Reformulation analytique :

   • Le papier introduit une représentation de type Crandall basée sur la décomposition du noyau $|z|^{-\nu}$ en « fonctions de Crandall » : une fonction supérieure $G_\nu$ (à décroissance superexponentielle) et une fonction inférieure $g_\nu$ (lisse).
   • En appliquant la formule de sommation de Poisson à la partie inférieure, les auteurs dérivent une représentation de la fonction zêta d'Epstein sous la forme d'une somme de deux sommes de réseaux à convergence superexponentielle. Cette formulation sert de base tant pour l'algorithme numérique que pour les preuves analytiques.
   • Les auteurs prouvent rigoureusement la continuation holomorphe conjointe de la fonction zêta d'Epstein dans l'exposant $\nu$ et les arguments vectoriels $x, y$ au sein de cônes complexes spécifiques. Ils caractérisent la structure des singularités, identifiant des pôles simples à $\nu = d$ (lorsque $y \in \Lambda^*$) et des singularités effaçables ailleurs.
   • Une fonction zêta d'Epstein régularisée est définie en soustrayant la singularité de puissance (dérivée de la transformée de Fourier distributionnelle de $| \cdot |^{-\nu}$). Cette régularisation élimine les singularités sans annulation catastrophique, permettant une évaluation efficace et une pré-computation par interpolation.

2. Algorithme numérique (EpsteinLib) :

   • Les auteurs implémentent une bibliothèque C, EpsteinLib, avec des liaisons Python, Julia et Mathematica.
   • Prétraitement : L'algorithme redimensionne les réseaux pour obtenir un volume unitaire et projette les vecteurs dans des cellules de réseau élémentaires pour exploiter les symétries de translation et d'échelle.
   • Cas particuliers : Un traitement explicite est fourni pour les singularités effaçables (par exemple, $\nu \in -2\mathbb{N}_0$) et les pôles ($\nu = d$).
   • Troncation et Sommation : Des bornes d'erreur rigoureuses sont dérivées pour la troncation des sommes à décroissance superexponentielle. Pour atténuer les erreurs d'arrondi de virgule flottante dans les dimensions élevées, les auteurs emploient une sommation compensée de Kahan.
   • Fonctions Gamma Incomplètes : Une implémentation personnalisée et performante de la fonction gamma incomplète est développée pour gérer les arguments négatifs et la plage complète de paramètres avec une précision double complète, surmontant les limitations des bibliothèques existantes.

Principales contributions

• Premier traitement analytique complet : Le papier fournit la première preuve de la continuation holomorphe conjointe de la fonction zêta d'Epstein dans tous les paramètres et offre une caractérisation complète de la structure de ses singularités.
• Algorithme à convergence superexponentielle : Un nouvel algorithme est introduit qui calcule la fonction pour des paramètres réels arbitraires dans n'importe quelle dimension avec des bornes d'erreur rigoureuses.
• Élimination des singularités : Une méthode de décomposition de la fonction en une partie singulière et une partie analytique régularisée est présentée, permettant l'évaluation des intégrales et le traitement des anomalies de type Wood sans annulation catastrophique.
• Implémentation EpsteinLib : La première bibliothèque open-source de haute performance pour calculer la fonction zêta d'Epstein pour des arguments réels arbitraires, présentant une évaluation en précision double complète et des benchmarks étendus.

Résultats

• Précision et Performance : L'implémentation atteint une précision de précision double complète à travers les dimensions $d=1, 2, 3, 4, 6, 8$ et toute la plage de paramètres $\nu \in [-12.5, 12.5]$. Les tests comparatifs par rapport aux formules analytiques connues et aux implémentations en précision arbitraire confirment la précision de l'algorithme.
• Efficacité : L'algorithme démontre une convergence superexponentielle. Par exemple, les sommes 2D sont calculées en moins de 20 microsecondes, et les sommes 8D en moins de 0,5 seconde sur un ordinateur portable standard. L'utilisation de la sommation de Kahan assure la stabilité même dans les hautes dimensions où les erreurs d'arrondi s'accumulent typiquement.
• Applications démontrées :
  • Systèmes de spins quantiques : La méthode a été appliquée pour calculer les relations de dispersion quantiques dans des systèmes de spins 3D à interactions à longue portée, capturant avec succès les comportements de mise à l'échelle anormaux ($\omega(k) \sim |k|^{\nu-d}$) pour $d < \nu < d+2$.
  • Effet Casimir : La bibliothèque a été utilisée pour calculer les énergies de Casimir dans des géométries 3D, reproduisant avec précision les lois de force asymptotiques et analysant la minimisation de l'énergie dans des boîtes déformées.

Signification
Le papier affirme que ce travail établit la fonction zêta d'Epstein comme une fonction spéciale standard pour la physique computationnelle et les mathématiques. En fournissant un fondement analytique rigoureux et une implémentation robuste et performante, les auteurs permettent la simulation efficace de systèmes à interactions à longue portée qui étaient auparavant trop coûteux en calcul ou numériquement instables. La capacité de gérer des réseaux généraux, des paramètres réels arbitraires et des singularités sans erreurs d'annulation facilite de nouvelles recherches en chimie théorique (stabilité de la matière), en physique de la matière condensée (diagrammes de phase quantiques, escaliers du diable) et dans les problèmes de diffusion d'ondes. Les auteurs notent que la bibliothèque est déjà utilisée par la communauté de la matière condensée quantique pour étudier des phénomènes exotiques, tels que les plateaux de magnétisation fractionnaire et la fusion des escaliers du diable.