Solution of Canonical Differential Equations for Integrals on Arbitrary Geometries

Technische Samenvatting: Oplossing van Canonieke Differentiaalvergelijkingen voor Integralen op Willekeurige Geometrieën

Probleemstelling
De evaluatie van multi-loop scattering amplitudes in kwantumveldentheorie berust vaak op het reduceren van Feynman-integralen tot een set master-integralen (MI's) via integratie-door-partiële identiteiten (IBP). Een krachtige methode voor het bepalen van deze MI's is het oplossen van systemen van eerste-orde lineaire differentiaalvergelijkingen met betrekking tot kinematische invarianten en massa's. De introductie van $\epsilon$-gefactoriseerde canonieke vormen heeft dit proces voor polylogarithmische gevallen gestroomlijnd, waarbij de connectiematrix bestaat uit $d\log$-vormen met rationale coëfficiënten.

Echter, voor problemen die betrokken zijn bij massieve deeltjes en meerdere schalen, strekt de onderliggende geometrie zich vaak uit voorbij eenvoudige polylogaritmen naar elliptische curven, Calabi-Yau-variëteiten en structuren met een hogere genus. In deze gevallen vereist het bereiken van een canonieke vorm een transformatie die gebruikmaakt van transcendente functies (perioden van de geometrie). Hoewel dergelijke vormen bestaan, ontstaat er een praktische flessenhals: de numerieke evaluatie van deze systemen wordt gehinderd door het gebrek aan gesloten vorm-expressies voor de vereiste speciale functies en de complexiteit van hun analytische continuatie. Bestaande numerieke benaderingen vertrouwen ofwel op reeksexpansies (beperkt in multiscale regimes) of vereisen de expliciete implementatie van complexe speciale functies, wat niet universeel haalbaar is.

Methodologie
De auteurs stellen een strategie voor om de noodzaak voor de expliciete evaluatie van transcendente functies te omzeilen. De kernobservatie is dat de speciale functies die verschijnen in de transformatie naar een canonieke basis niet arbitrair zijn; ze zijn afkomstig van het oorspronkelijke systeem van master-integralen, die voldoen aan lineaire differentiaalvergelijkingen met rationale coëfficiënten.

Consequent demonstreren de auteurs dat deze transcendente bouwstenen (bijv. elliptische perioden, vierkantswortels) zelf voldoen aan systemen van differentiaalvergelijkingen met rationale coëfficiënten. De voorgestelde methode omvat:

1. Augmentatie: Het construeren van een vergroot systeem van differentiaalvergelijkingen dat de canonieke differentiaalvergelijking voor de master-integralen koppelt aan de hulp-differentiaalvergelijkingen die de speciale functies beheersen.
2. Rationalisering: Dit vergrote systeem bevat uitsluitend rationale coëfficiënten. De afgeleiden van de speciale functies worden direct berekend uit hun eigen rationale differentiaalvergelijkingen in plaats van uit gesloten vorm-expressies.
3. Numerieke Integratie: Het gehele systeem wordt numeriek opgelost langs een gecomplexeerde pad in de kinematische ruimte om singulariteiten te vermijden. De analytische continuatie van de speciale functies wordt automatisch afgehandeld door de integratie van hun differentiaalvergelijkingen, waardoor ad-hoc algoritmen voor analytische continuatie overbodig worden.
4. Randvoorwaarden: De methode vereist randwaarden voor de speciale functies op een niet-singulier punt. Deze waarden kunnen willekeurig gekozen worden (bijv. ingesteld op 1), mits ze consistent zijn met eventuele algebraïsche relaties die gebruikt worden om de connectiematrix te vereenvoudigen. De canonieke master-integralen worden vervolgens geroteerd naar de canonieke basis met behulp van deze randwaarden.

Belangrijkste Bijdragen

• Theoretisch Inzicht: Het artikel stelt vast dat de obstructie voor het numeriek oplossen van canonieke systemen met transcendente entries slechts schijnbaar is. Door de speciale functies te behandelen als dynamische variabelen die hun eigen rationale differentiaalvergelijkingen bevredigen, wordt het probleem gereduceerd tot het oplossen van een groter, puur rationaal systeem.
• Software Implementatie: De auteurs bieden een standalone C++ pakket (diffeqs) aan die is ontworpen voor hoogwaardige numerieke evaluatie. De software:
  • Is geschreven in C++ en vertrouwt enkel op de standaard Boost-bibliotheek, waardoor afhankelijkheden van commerciële software (zoals Mathematica) tijdens runtime worden vermeden.
  • Ondersteunt willekeurige floating-point precisie via de QD-bibliotheek en getypeerde templates.
  • Implementeert een modulaire interface waarmee gebruikers de connectiematrix, de vectorveld voor speciale functies en het integratiepad kunnen definiëren.
  • Bevat functies voor complexe paddeformatie om numerieke stabiliteit en analytische continuatie te waarborgen.
• Realistische Toepassing: De methode wordt toegepast op een niet-triviaal fysiek geval: de twee-loop master-integralen voor di-foton, di-jet en $\gamma$+jet hadroproductie via een zware-quark lus. Dit systeem omvat 220 master-integralen, 16 vierkantswortels en transcendente functies gerelateerd aan elliptische perioden.

Resultaten
De auteurs hebben de implementatie getest op drie specifieke fase-ruimte punten.

• Nauwkeurigheid: In dubbele precisie bereikte de methode absolute fouten van $O(10^{-6})$ en $O(10^{-9})$ voor de geteste punten. In quadrupel precisie werden de fouten gereduceerd tot $O(10^{-8})$ en $O(10^{-19})$ respectievelijk, wat de correctheid van de analytische continuatie en de numerieke stabiliteit van de aanpak valideert.
• Prestaties: De evaluatietijd voor de connectiematrix bleek efficiënt te zijn (sub-milliseconde tot enkele milliseconden, afhankelijk van precisie en hardware). De auteurs merken op dat de quadrupel-precisie versie ongeveer 10 keer trager is dan de dubbele-precisie versie met betrekking tot de evaluatie van de connectiematrix.
• Vergelijking: De resultaten zijn gecontroleerd tegen waarden verkregen met de AMFlow-package, wat de consistentie bevestigt. De prestaties waren superieur op Apple Silicon chips vergeleken met AMD-processoren.

Betekenis en Claims
Het artikel claimt dat deze strategie canonieke differentiaalvergelijkingen een universeel toepasbare tool maakt voor willekeurige Feynman-integralen, inclus_ing die met elliptische en Calabi-Yau geometrieën. De betekenis ligt in het vermogen om numerieke fenomenologie uit te voeren zonder dat er closed-form implementaties van complexe speciale functies nodig zijn. Door de geometrische bouwstenen te propageren via hun eigen rationale differentiaalvergelijkingen, vermijdt de methode de "black box" implementatie van speciale functies en de daarmee gepaard gaande moeilijkheden van analytische continuatie.

De auteurs positioneren dit werk als een praktische stap naar de automatisering van hogere-orde perturbatieve berekeningen voor complexe scattering-processen. Zij benadrukken dat hoewel de huidige implementatie functioneel is zonder zware optimalisatie, verdere winst kan worden behaald door expressie-optimalisatie en partial-fraction decomposition strategieën. De methode wordt gepresenteerd als klaar voor grootschalige toepassingen in collider-fenomenologie, met name voor processen die meerdere schalen en massieve deeltjes bevatten waar traditionele reeksexpansie-methoden moeite hebben.