CHESS: CHEbyshev pSeudo-Spectral transport for Feynman integral differential equations

Sintesi Tecnica: Chess (CHEbyshev pSeudo-Spectral transport per equazioni differenziali di integrali di Feynman)

Definizione del Problema
La valutazione degli integrali master di Feynman multi-loop rappresenta un collo di bottiglia computazionale nella fisica dei collider di precisione. Man mano che i calcoli si spostano verso ordini di loop e molteplicità più elevati, le famiglie di integrali coinvolgono numerose scale cinematiche, grandi insiemi di integrali master e strutture di singolarità complesse. Sebbene esistano varie strategie (ad esempio, decomposizione dei settori, flusso di massa ausiliaria, espansioni di serie locali), un problema numerico sottostante specifico rimane impegnativo: il trasporto ad alta precisione unidimensionale di equazioni differenziali $\epsilon$-fattorizzate lungo un percorso prestabilito nello spazio cinematico.

I metodi di propagazione locale standard richiedono spesso molti segmenti di espansione o passi piccoli per controllare l'alta precisione e il comportamento agli estremi, in particolare quando la matrice di coefficienti richiamata (pull-back) coinvolge combinazioni sparse di matrici costanti e pull-back scalari di una forma singola, o quando il sistema include ingredienti algebrici o ellittici. La sfida consiste nel risolvere il sistema di equazioni differenziali ordinarie (ODE) risultante in modo efficiente, mantenendo l'accuratezza spettrale e gestendo potenziali singolarità regolari spurie agli estremi di integrazione.

Metodologia
Il documento introduce Chess (CHEbyshev pSeudo-Spectral transport), un pacchetto in Wolfram Language progettato per risolvere il problema del trasporto unidimensionale:
$$\frac{d}{dt}J(t, \epsilon) = \epsilon B(t)J(t, \epsilon)$$
dove $J$ è un vettore di integrali master (o i loro coefficienti dell'espansione in $\epsilon$) e $B(t)$ è la matrice di coefficienti richiamata.

La metodologia centrale si basa sulla collocazione spettrale di Chebyshev–Lobatto. Il metodo sfrutta l'analiticità della connessione richiamata nel parametro del percorso complesso. Se la matrice di coefficienti $B(t)$ è analitica in un intorno complesso dell'intervallo di trasporto, la soluzione converge geometricamente con il numero di nodi di collocazione.

Caratteristiche algoritmiche chiave:

1. Gerarchia $\epsilon$-Fattorizzata: Il sistema è espanso in una finestra di Laurent finita. A causa della forma $\epsilon$-fattorizzata, la gerarchia dei coefficienti è triangolare. Ogni coefficiente $J^{(r)}$ è ottenuto dal precedente $J^{(r-1)}$ risolvendo un sistema del primo ordine disomogeneo. Ciò consente al solver di utilizzare una singola matrice di nodi scalare fattorizzata agendo su molti termini membri, invece di costruire un sistema monolitico per tutti i nodi e gli integrali simultaneamente.
2. Modalità di Input Doppia:
   • Modalità Dati Preparati (Prepared-Data Mode): L'utente fornisce matrici sparse costanti $A_\alpha$ e pull-back scalari precomputati $b_\alpha(t)$. Il pacchetto assembla $B(t) = \sum A_\alpha b_\alpha(t)$ nei nodi di collocazione. Questo separa le valutazioni scalari dipendenti dal percorso dai dati matriciali costanti, consentendo il caching e l'assemblaggio sparso.
   • Modalità Matrice Diretta (Direct-Matrix Mode): L'utente fornisce un valutatore numerico per la matrice completa $B(t)$. Ciò consente al solver di gestire sistemi con periodi ellittici o altre funzioni speciali senza richiedere una specifica decomposizione logaritmica.
3. Regolarizzazione degli Estremi: Per gestire singolarità regolari spurie agli estremi (dove $B(t)$ può divergere ma la soluzione rimane finita), Chess implementa una regolarizzazione basata sul residuo.
   • Per un estremo destro, la riga dell'equazione differenziale è sostituita da una condizione di sorgente finita derivata dal residuo e dall'estrazione della parte finita della connessione.
   • Per un estremo sinistro, il valore limite noto viene rimosso e un "sollevamento" (lifted) dell'incognita (il differenziale unidirezionale regolare) viene risolto, sostituendo l'equazione differenziale singolare con una condizione di derivata finita.
4. Parallelizzazione: La parallelizzazione è applicata solo a compiti indipendenti, come la valutazione dei pull-back scalari ai nodi o l'estrazione dei dati degli estremi. La ricorsione $\epsilon$-sequenziale è preservata, poiché lo strato $r$ dipende dallo strato $r-1$.

Contributi Chiave

• Efficienza Algoritmica: Utilizzando la struttura triangolare dei sistemi $\epsilon$-fattorizzati e l'assemblaggio di matrici sparse, Chess evita l'overhead di memoria e computazionale dei solver monolitici. Raggiunge una rapida convergenza dei nodi, spesso toccando i limiti di precisione della macchina con relativamente pochi nodi di collocazione.
• Flessibilità: Il pacchetto non è limitato ai sistemi logaritmici. L'interfaccia a matrice diretta permette di trasportare sistemi contenenti ingredienti ellittici o modulari, purché sia disponibile un affidabile valutatore numerico per $B(t)$.
• Robustezza: L'implementazione include la gestione specifica per gli estremi con singolarità regolari spurie, un problema comune nei trasporti cinematici che può destabilizzare i metodi di collocazione standard.
• Implementazione: Il pacchetto è scritto in Wolfram Language, offrendo un'interfaccia di alto livello per utenti che hanno già completato la riduzione degli integrali, la costruzione della base e la determinazione delle costanti di contorno.

Risultati e Benchmark
Il documento valida Chess attraverso tre gruppi di benchmark:

1. Doppio Pentagono a Sei Punti a Due Loop: Trasportato nella regione fisica e confrontato con dati AMFlow indipendenti. Con 24 nodi, la discrepanza era di $1.37 \times 10^{-31}$, diminuendo a $10^{-101}$ con 96 nodi, saturando alla precisione dei dati di riferimento.
2. Famiglie a Cinque Punti a Tre Loop: Quattro famiglie planari prive di massa (PBB, BPB, BHB, PBP) sono state testate con regolarizzazione dell'estremo sinistro. Il metodo ha dimostrato una raffinazione stabile dei nodi con differenze intorno a $10^{-70}$ (circa 70 cifre). Un confronto a protocollo fisso con il codice di serie locale DiffExp ha mostrato che Chess è da 8.8 a 17.9 volte più veloce in termini di tempo di esecuzione (wall time). Per il sistema più grande (734 integrali master), Chess ha anche dimostrato un minor utilizzo di memoria dell'albero di processo rispetto a DiffExp.
3. Triangolo Non-Planare (Caso Ellittico): Una famiglia di triangoli non-planari a due loop con periodi ellittici è stata trasportata utilizzando la modalità a matrice diretta. Il metodo ha raggiunto la convergenza geometrica fino al limite di precisione a 200 cifre del riferimento AMflow indipendente, raggiungendo errori inferiori a $10^{-200}$ con 512 nodi.

Significato e Ambito
Il documento posiziona Chess come uno strumento specializzato ad alte prestazioni all'interno dell'ecosistema più ampio del calcolo degli integrali di Feynman. Afferma esplicitamente che il pacchetto non esegue la riduzione degli integrali, non costruisce basi canoniche, non determina le costanti di contorno né seleziona percorsi analitici. Questi passaggi preparatori si assumono siano stati completati dall'utente.

Il significato di Chess risiede nella sua capacità di eseguire lo step finale di trasporto, computazionalmente intensivo, con alta precisione ed efficienza. Sfrutta le proprietà analitiche del parametro del percorso per ottenere la convergenza spettrale, rendendolo particolarmente efficace per sistemi multi-scala di grandi dimensioni dove i metodi tradizionali di serie locali potrebbero richiedere una segmentazione eccessiva. Gli autori sottolineano che il successo del metodo dipende dall'analiticità delle funzioni di coefficiente fornite; non è un algoritmo automatico di pianificazione del percorso per attraversare singolarità fisiche. Il pacchetto funge da robusto strato di trasporto per utenti che hanno già stabilito il necessario framework analitico.