First look at the evaluation of two-loop Feynman integrals for radiative return processes

Sintesi Tecnica: Prima analisi della valutazione degli integrali di Feynman a due loop per i processi di ritorno radiativo

Problema
I test di precisione del Modello Standard presso i collisionatori elettrone–positrone a bassa energia (ad es. BaBar, Belle II, BESIII, KLOE) si basano su previsioni teoriche accurate per i processi di ritorno radiativo, come $e^+e^- \to \pi^+\pi^-\gamma$ e $e^+e^- \to \mu^+\mu^-\gamma$. Questi processi sono critici per determinare il contributo della polarizzazione del vuoto adronico al momento magnetico anomalo del muone ($(g-2)_\mu$). Per risolvere le persistenti tensioni tra le misurazioni sperimentali e le previsioni teoriche, sono necessarie correzioni QED al Next-to-Next-to-Leading Order (NNLO).

Il principale collo di bottiglia nel raggiungere l'accuratezza NNLO per questi processi, mantenendo la piena dipendenza dalla massa dell'elettrone, è la valutazione degli integrali di Feynman a quattro punti a due loop. La presenza di propagatori massivi introduce una significativa complessità algebrica nelle riduzioni tramite Integrazione per Parti (IBP) e nelle Equazioni Differenziali (DE). Inoltre, questi integrali coinvolgono funzioni speciali associate a geometrie di genere superiore, specificamente curve ellittiche e, in alcuni casi, radici quadrate annidate. Costruire una base canonica per questi integrali — dove le DE sono $\epsilon$-fattorizzate ed espresse in termini di forme dlog — è matematicamente impegnativo e spesso non strettamente necessario per la fenomenologia numerica.

Metodologia
Gli autori si concentrano sul contributo della radiazione di stato iniziale (ISR) per il processo $e^+e^- \to \gamma\gamma^*$, che riduce il problema alla valutazione di ampiezze genuine a due loop a quattro punti. Considerano tre famiglie di integrali planari, denotate come PL1, PL2 e PL3, che contengono un totale di 212 integrali master (MI), inclusi quelli associati a cinque curve ellittiche non isomorfe e un settore con una radice quadrata annidata.

Invece di tentare di costruire una base canonica completa per tutti i settori (il che è difficile per le geometrie ellittiche e con radici annidate), gli autori adottano una strategia sviluppata in lavori precedenti [48, 49]. Il loro approccio prevede:

1. Costruzione della Base: Costruiscono una base di MI tale che il sistema di DE risultante sia polinomiale nel regolatore dimensionale $\epsilon$. Per i settori polilogaritmici, mirano a forme canoniche. Per i settori ellittici, identificano i MI associati a differenziali ellittici del primo, secondo e terzo tipo.
2. Gestione delle Geometrie Ellittiche: Per i settori in cui la natura ellittica non è manifesta nell'integrale d'angolo (ad es. settori kite a quattro punti), utilizzano super-settori riducibili e analizzano gli operatori di Picard-Fuchs per identificare la geometria e costruire le basi appropriate.
3. Gestione delle Radici Quadrate Annidate: In un settore specifico di PL1 che coinvolge una radice quadrata annidata, dimostrano che, sebbene esista una trasformazione per $\epsilon$-fattorizzare le DE (che comporta un dupletto di radici annidate), preferiscono lavorare con una base che mantenga le DE lineari in $\epsilon ed eviti l'introduzione di strutture trascendenti complesse, affidandosi invece a funzioni razionali e radici quadrate semplici.
4. Integrazione Numerica: Le DE vengono risolte numericamente. Le condizioni al contorno sono determinate utilizzando il pacchetto AMFlow in un punto dello spazio delle fasi scelto per simulare configurazioni cinematiche a bassa energia. L'integrazione è eseguita utilizzando un'implementazione in-house nel linguaggio di programmazione Julia, impiegando integratori non rigidi (Tsit5 e Vern).

Contributi Chiave e Risultati

• Framework Sistematico: Il documento presenta la prima costruzione sistematica di DE per le famiglie di integrali planari rilevanti per i processi di ritorno radiativo NNLO con piena dipendenza dalla massa dell'elettrone.
• Trattamento dei Settori Ellittici: Gli autori identificano e trattano con successo i settori associati a cinque distinte curve ellittiche. Forniscono basi esplicite per questi settori, inclusi i difficili topologie kite a quattro punti dove il comportamento ellittico deriva da integrali con punti (potenze dei propagatori $>1$) piuttosto che dagli integrali d'angolo.
• Analisi delle Radici Quadrate Annidate: Il lavoro dettaglia il trattamento di un settore in PL1 contenente una radice quadrata annidata. Costruisce una base che soddisfa DE polinomiali in $\epsilon$ senza richiedere l'introduzione esplicita della struttura della radice annidata nelle equazioni differenziali, facilitando la stabilità numerica.
• Validazione Numerica: Gli autori dimostrano valutazioni numeriche stabili attraverso la regione fisica rilevante per gli esperimenti a bassa energia ($s \in \{1.0404, 16, 100\} \text{ GeV}^2$).
  • Performance: I tempi di valutazione variano da 40 a 800 ms per punto di spazio delle fasi.
  • Precisione: I confronti con AMFlow mostrano un accordo di almeno otto cifre significative in specifici punti di test. I confronti con il pacchetto DiffExp (utilizzando la doppia precisione) mostrano circa 5 cifre di accordo per le famiglie PL1 e PL2 e 2 cifere per la famiglia PL3 nei casi peggiori.
  • Robustezza: Il metodo rimane stabile nonostante le grandi gerarchie di scale ($m_e^2/s \sim 10^{-7} - 10^{-8}$).

Significatività
Il documento sostiene di fornire "elementi costitutivi essenziali per le previsioni NNLO per i processi di ritorno radiativo". Avendo costruito con successo e integrato numericamente le DE per i settori planari, gli autori hanno superato l'ostacolo computazionale primario per la componente di radiazione di stato iniziale dell'ampiezza NNLO.

Il lavoro è modesto nella sua portata, dichiarando esplicitamente di coprire solo le famiglie planari. L'estensione alle famiglie non planari e l'assemblaggio dell'ampiezza completa sono indicati come lavori in corso. La significatività risiede nella dimostrazione che un approccio polinomiale in $\epsilon$, combinato con l'integrazione numerica, è una strategia valida ed efficiente per gestire le geometrie complesse (curve ellittiche e radici quadrate annidate) intrinseche nei calcoli a due loop massivi, evitando la necessità di forme canoniche completamente analitiche che sono attualmente difficili da ottenere per tali geometrie.