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

Technische Samenvatting: Berekening en Eigenschappen van de Epstein-zetafunctie

Probleemstelling
De Epstein-zetafunctie, $Z_{\Lambda, \nu}(x, y)$, generaliseert de Riemann-zetafunctie naar oscillerende roosterzommen in hogere dimensies. Het is een cruciaal instrument in de theoretische natuurkunde en chemie voor het modelleren van lang reikende interagerende klassieke en kwantum veel-deeltjes-systemen, waaronder Madelung-constanten, Casimir-energieën en spingolf-dispersierelaties. Ondanks het belang ervan ontbrak het tot nu toe aan een rigoureuze en uitgebreide behandeling van de berekening en de analytische eigenschappen ervan. Bestaande methoden, zoals de Chowla–Selberg-formule of de representatie van Crandall, kampen vaak met beperkingen: ze kunnen beperkt zijn tot speciale gevallen van roostervectoren of parameters, lijden onder numerieke instabiliteit nabij singulariteiten (bijv. negatieve argumenten in incomplete gammafuncties), of missen een uniforme bibliotheekimplementatie. Bovendien was de gezamenlijke holomorfe voortzetting en de precieze complexe singulariteitsstructuur van de functie in alle parameters ($\nu, x, y$) nog niet rigoureus vastgesteld, wat de ontwikkeling van robuuste numerieke algoritmen die catastrofale annulering voorkomen, in de weg stond.

Methodologie
De auteurs pakken deze hiaten aan door een combinatie van een rigoureuze analytische herformulering en hoogwaardig numeriek algoritmeontwerp:

1. Analytische Herformulering:

   • Het artikel introduceert een Crandall-achtige representatie gebaseerd op de decompositie van de kernel $|z|^{-\nu}$ in "Crandall-functies": een bovenste functie $G_\nu$ (superexponentieel vervallend) en een onderste functie $g_\nu$ (gladde functie).
   • Door de Poisson-sommatieformule toe te passen op het onderste deel, leiden de auteurs een representatie van de Epstein-zetafunctie af als de som van twee superexponentieel convergerende roosterzommen. Deze formulering dient als basis voor zowel het numerieke algoritme als de analytische bewijzen.
   • De auteurs bewijzen rigoureus de gezamenlijke holomorfe voortzetting van de Epstein-zetafunctie in de exponent $\nu$ en de vectorargumenten $x, y$ binnen specifieke complexe kegels. Ze karakteriseren de singulariteitsstructuur en identificeren eenvoudige polen bij $\nu = d$ (wanneer $y \in \Lambda^*$) en verwijderbare singulariteiten elders.
   • Een geregulariseerde Epstein-zetafunctie wordt gedefinieerd door de machtswet-singulariteit af te trekken (afgeleid van de distributionele Fourier-transformatie van $| \cdot |^{-\nu}$). Deze regularisatie verwijdert singulariteiten zonder catastrofale annulering, wat stabiele evaluatie en efficiënte precomputatie via interpolatie mogelijk maakt.

2. Numeriek Algoritme (EpsteinLib):

   • De auteurs implementeren een C-bibliotheek, EpsteinLib, met Python-, Julia- en Mathematica-bindings.
   • Preprocessing: Het algoritme schaalt roosters naar een eenheidsvolume en projecteert vectoren in elementaire roostercellen om gebruik te maken van translatie- en schaal-symmetrieën.
   • Speciale Geval: Expliciete afhandeling wordt geboden voor verwijderbare singulariteiten (bijv. $\nu \in -2\mathbb{N}_0$) en polen ($\nu = d$).
   • Truncatie en Sommatie: Rigoureuze foutmargen worden afgeleid voor het afkappen van de superexponentieel vervallende sommen. Om floating-point afrondingsfouten in hoge dimensies te beperken, gebruiken de auteurs Kahan-gecompenseerde sommatie.
   • Incomplete Gammafuncties: Een aangepaste, hoogwaardige implementatie van de incomplete gammafunctie is ontwikkeld om negatieve argumenten en het volledige parametrange aan te kunnen met volledige double-precisie, waarmee de beperkingen van bestaande bibliotheken worden overwonnen.

Belangrijkste Bijdragen

• Eerste Uitgebreide Analytische Behandeling: Het artikel biedt het eerste bewijs van de gezamenlijke holomorfe voortzetting van de Epstein-zetafunctie in alle parameters en biedt een volledige karakterisering van de singulariteitsstructuur.
• Superexponentieel Convergerend Algoritme: Een nieuw algoritme wordt geïntroduceerd dat de functie berekent voor willekeurige reële parameters in elke dimensie met rigoureuze foutmargen.
• Verwijdering van Singulariteiten: Een methode voor het decomponeren van de functie in een singuliere component en een geregulariseerd analytisch deel wordt gepresenteerd, wat de evaluatie van integralen en de behandeling van Wood-type anomalieën mogelijk maakt zonder catastrofale annulering.
• EpsteinLib Implementatie: De eerste hoogwaardige, open-source bibliotheek voor het berekenen van de Epstein-zetafunctie voor willekeurige reële argumenten, met volledige-precisie evaluatie en uitgebreide benchmarking.

Resultaten

• Nauwkeurigheid en Prestaties: De implementatie bereikt volledige double-precision nauwkeurigheid over dimensies $d=1, 2, 3, 4, 6, 8$ en het gehele parametrange $\nu \in [-12.5, 12.5]$. Benchmarks tegen bekende analytische formules en precisie-implementaties met willekeurige precisie bevestigen de nauwkeurigheid van het algoritme.
• Efficiëntie: Het algoritme vertoont superexponentiële convergentie. Zo worden 2D-zommen berekend in minder dan 20 microseconden, en 8D-zommen in minder dan 0,5 seconden op een standaard laptop. Het gebruik van Kahan-sommatie zorgt voor stabiliteit, zelfs in hoge dimensies waar afrondingsfouten gewoonlijk accumuleren.
• Toegepaste Demonstraties:
  • Kwantum Spinsystemen: De methode werd toegepast om kwantum dispersierelaties in 3D spin-systemen met lang reikende interacties te berekenen, waarbij de anomalieën in schaling ($\omega(k) \sim |k|^{\nu-d}$) voor $d < \nu < d+2$ succesvol werden gevangen.
  • Casimir-effect: De bibliotheek werd gebruikt om Casimir-energieën in 3D-geometrieën te berekenen, waarbij de asymptotische krachtwetten nauwkeurig werden gereproduceerd en de energie-minimalisatie in gedeformeerde dozen werd geanalyseerd.

Betekenis
Het artikel beweert dat dit werk de Epstein-zetafunctie vestigt als een standaard speciale functie voor computationele fysica en wiskunde. Door een rigoureuze analytische fundering en een robuuste, hoogwaardige implementatie te bieden, stellen de auteurs de efficiënte simulatie van lang reikende interagerende systemen mogelijk die voorheen computationeel onhaalbaar of numeriek instabiel waren. Het vermogen om algemene roosters, willekeurige reële parameters en singulariteiten zonder annulering-fouten te behandelen, vergemakkelijkt nieuw onderzoek in theoretische chemie (stabiliteit van materie), gecondenseerde materie fysica (kwantum fase-diagrammen, devil's staircases) en golfverstrooiingsproblemen. De auteurs merken op dat de bibliotheek al wordt gebruikt door de kwantum-condensed-matter gemeenschap voor het bestuderen van exotische fenomenen, zoals fractionele magnetisatieplateaus en het smelten van devil's staircases.