Solution of Canonical Differential Equations for Integrals on Arbitrary Geometries

Sintesi Tecnica: Risoluzione di Equazioni Differenziali Canoniche per Integrali su Geometrie Arbitrarie

Enunciato del Problema
La valutazione di ampiezze di scattering multi-loop nella teoria quantistica dei campi si basa spesso sulla riduzione degli integrali di Feynman a un insieme di integrali master (MI) tramite identità di integrazione per parti. Un metodo potente per determinare questi MI consiste nel risolvere sistemi di equazioni differenziali lineari del primo ordine rispetto alle invarianti cinematiche e alle masse. L'introduzione di forme canoniche $\epsilon$-fattorizzate ha snellito questo processo per i casi polilogaritmici, dove la matrice di connessione consiste in forme $d\log$ con coefficienti razionali.

Tuttavia, per problemi che coinvolgono particelle massive e molteplici scale, la geometria sottostante si estende spesso oltre i semplici polilogaritmi per includere curve ellittiche, varietà di Calabi-Yau e strutture di genere superiore. In questi casi, il raggiungimento di una forma canonica richiede una trasformazione che coinvolge funzioni trascendenti (periodi della geometria). Sebbene tali forme esistano, un collo di bottiglia pratico sorge dal fatto che la valutazione numerica di questi sistemi è ostacolata dalla mancanza di espressioni in forma chiusa per le funzioni speciali richieste e dalla complessità della loro continuazione analitica. Gli approcci numerici esistenti si basano o su espansioni in serie (limitate nei regimi multi-scala) o richiedono l'implementazione esplicita di complesse funzioni speciali, il che non è universalmente fattibile.

Metodologia
Gli autori propongono una strategia per aggirare la necessità di valutare esplicitamente le funzioni trascendenti. L'osservazione centrale è che le funzioni speciali che compaiono nella trasformazione verso una base canonica non sono arbitrarie; esse originano dal sistema originale di integrali master, che soddisfano equazioni differenziali lineari con coefficienti razionali.

Di conseguenza, gli autori dimostrano che questi mattoni trascendenti (ad esempio, periodi ellittici, radici quadrate) soddisfano essi stessi sistemi di equazioni differenziali con coefficienti razionali. Il metodo proposto prevede:

1. Augmentazione: Costruzione di un sistema ampliato di equazioni differenziali che accoppia l'equazione differenziale canonica per gli integrali master con le equazioni differenziali ausiliarie che governano le funzioni speciali.
2. Razionalizzazione: Questo sistema ampliato contiene solo coefficienti razionali. Le derivate delle funzioni speciali sono calcolate direttamente dalle loro equazioni differenziali razionali definitorie piuttosto che da espressioni in forma chiusa.
3. Integrazione Numerica: L'intero sistema viene risolto numericamente lungo un percorso complesso nello spazio cinematico per evitare singolarità. La continuazione analitica delle funzioni speciali è gestita automaticamente dall'integrazione delle loro equazioni differenziali, eliminando la necessità di algoritmi di continuazione analitica ad hoc.
4. Condizioni al Contorno: Il metodo richiede valori al contorno per le funzioni speciali in un punto non singolare. Tali valori possono essere scelti arbitrariamente (ad esempio, impostati a 1) a patto che siano coerenti con qualsiasi relazione algebrica utilizzata per semplificare la matrice di connessione. I master integrali canonici sono poi ruotati verso la base canonica utilizzando questi valori al contorno.

Contributi Chiave

• Intuizione Teorica: Il lavoro stabilisce che l'ostacolo alla risoluzione numerica di sistemi canonici con entrate trascendenti è solo apparente. Trattando le funzioni speciali come variabili dinamiche che soddisfano le proprie equazioni differenziali razionali, il problema si riduce alla risoluzione di un sistema più grande e puramente razionale.
• Implementazione Software: Gli autori forniscono un pacchetto C++ standalone (diffeqs) progettato per la valutazione numerica ad alte prestazioni. Il software:
  • È scritto in C++ e si affida solo alla libreria standard Boost, evitando dipendenze da software commerciali (ad esempio, Mathematica) a runtime.
  • Supporta l'arbitraria precisione del punto mobile tramite la libreria QD e tipi templated.
  • Implementa un'interfaccia modulare che consente agli utenti di definire la matrice di connessione, il campo vettoriale per le funzioni speciali e il percorso di integrazione.
  • Include funzionalità per la deformazione del percorso complesso per garantire la stabilità numerica e la continuazione analitica.
• Applicazione Realistica: Il metodo è applicato a un caso fisico non banale: gli integrali master a due loop per la produzione di fotone-fotone, jet-fotone e $\gamma$+jet tramite un loop di quark pesante. Questo sistema coinvolge 220 integrali master, 16 radici quadrate e funzioni trascendenti legate ai periodi ellittici.

Risultati
Gli autori hanno testato l'implementazione su tre punti specifici dello spazio delle fasi.

• Accuratezza: In doppia precisione, il metodo ha raggiunto errori assoluti di $O(10^{-6})$ e $O(10^{-9})$ per i punti testati. In quadrupla precisione, gli errori sono stati ridotti rispettivamente a $O(10^{-8})$ e $O(10^{-19})$, validando la correttezza della continuazione analitica e la stabilità numerica dell'approccio.
• Prestazioni: La valutazione del tempo per la matrice di connessione è risultata efficiente (da meno di un millisecondo a pochi millisecondi a seconda della precisione e dell'hardware). Gli autori notano che la versione a quadrupla precisione è circa 10 volte più lenta della versione a doppia precisione per quanto riguarda la valutazione della matrice di connessione.
• Confronto: I risultati sono stati verificati rispetto ai valori ottenuti con il pacchetto AMFlow, confermando la coerenza. Le prestazioni sono risultate superiori sui chip Apple Silicon rispetto ai processori AMD.

Significato e Rivendicazioni
Il documento sostiene che questa strategia rende le equazioni differenziali canoniche uno strumento universalmente applicabile per arbitrari integrali di Feynman, inclusi quelli che coinvolgono geometrie ellittiche e di Calabi-Yau. La significatività risiede nella capacità di eseguire la fenomenologia numerica senza richiedere implementazioni in forma chiusa di complesse funzioni speciali. Propagando i mattoni geometrici attraverso le proprie equazioni differenziali razionali, il metodo evita l'implementazione "black box" delle funzioni speciali e le relative difficoltà di continuazione analitica.

Gli autori posizionano questo lavoro come un passo pratico verso l'automazione dei calcoli perturbativi di ordine superiore per processi di scattering complessi. Essi sottolineano che, sebbene l'attuale implementazione sia funzionale senza ottimizzazioni pesanti, ulteriori guadagni possono essere ottenuti tramite l'ottimizzazione delle espressioni e strategie di decomposizione in frazioni parziali. Il metodo è presentato come pronto per applicazioni su larga scala nella fenomenologia dei collider, in particolare per processi che coinvolgono molteplici scale e particelle massive dove i metodi tradizionali di espansione in serie faticano.