Automated NRQCD and NRQED simulations of quarkonium and leptonium production with P-wave states and physical-mass effects

Sintesi Tecnica: Simulazioni Automatizzate di NRQCD e NRQED della Produzione di Quarkonium e Leptonium con Stati a Onda P ed Effetti di Massa Fisica

Definizione del Problema
La produzione e il decadimento di quarkonia pesanti (stati legati di quark pesanti) e leptonia (stati legati di leptoni) fungono da laboratori critici per testare la Cromodinamica Quantistica (QCD) e l'Elettrodinamica Quantistica (QED) nel regime non relativistico. Sebbene i framework della QCD non relativistica (NRQCD) e della QED non relativistica (NRQED) forniscano una rigorosa fattorizzazione dei coefficienti a breve distanza e degli elementi di matrice a lunga distanza (LDME), l'automazione della generazione di eventi per questi processi è stata storicamente limitata. Strumenti precedenti, come MadOnia e HELAC-Onia, si sono concentrati principalmente su stati a onda S o hanno carezato di una piena automazione per un numero arbitrario di stati a onda P. Inoltre, i calcoli perturbativi convenzionali spesso approssimano la massa dello stato legato come la somma delle masse dei suoi costituenti ($M_B \approx m_Q + m_{\bar{Q}}$). Questa approssimazione introduce incertezze sistematiche, in particolare nelle regioni soglia cinematiche e nei processi che coinvolgono molteplici stati legati con costanti identiche ma masse fisiche differenti (ad esempio, $J/\psi + \psi(2S)$), dove la differenza di massa fisica è significativa.

Metodologia
Gli autori presentano MadSONS (MadGraph Simulations Of NRQFT States), un nuovo modulo per il framework MadGraph5_aMC@NLO (MG5_aMC). Questo modulo estende l'implementazione precedentemente stabilita per gli stati a onda S includendo processi di produzione a livello di albero (tree-level) arbitrari che coinvolgono stati a onda P ed effetti di massa fisica.

1. Implementazione P-wave tramite Numeri Duali:
   A differenza degli stati a onda S, gli stati a onda P hanno una funzione d'onda nulla all'origine, il che significa che i loro tassi di produzione dipendono dalla derivata prima della funzione d'onda. Nello spazio dei momenti, ciò corrisponde a proiettori di momento angolare orbitale che coinvolgono derivate rispetto al momento relativo.

   • Gli autori implementano queste derivate utilizzando numeri duali multivariati, una forma di differenziazione automatica. Questo approccio fornisce derivate esatte e numericamente stabili senza fare affidamento su differenziazione simbolica o approssimazioni a differenze finite.
   • Il formalismo dei numeri duali consente la propagazione delle derivate attraverso la catena di generazione dell'ampiezza, permettendo la valutazione diretta dei proiettori di momento angolare orbitale ($P_L$) e di momento angolare totale ($P_J$) all'interno del calcolo dell'elemento di matrice.
   • Questo metodo è progettato per essere scalabile per future estensioni a momenti angolari orbitali più elevati (ad esempio, stati a onda D) e requisiti di derivata di ordine superiore per calcoli al Prossimo Ordine (NLO).

2. Rimescolamento del Momento per Masse Fisiche (Momentum Reshuffling):
   Per affrontare la discrepanza tra la massa fisica dello stato legato e la somma delle masse dei costituenti, gli autori introducono una procedura di rimescolamento del momento.

   • Generazione dello Spazio delle Fasi: Gli eventi vengono generati utilizzando le masse fisiche degli stati legati, garantendo il corretto volume dello spazio delle fasi disponibile e le soglie cinematiche.
   • Valutazione dell'Elemento di Matrice: L'elemento di matrice viene valutato utilizzando le condizioni on-shell delle particelle costituenti (somma delle masse dei costituenti), come richiesto dal formalismo NRQFT per i coefficienti a breve distanza.
   • Algoritmi di Rimescolamento: La mappatura tra i momenti fisici e i momenti a livello di costituenti è ottenuta tramite due distinti schemi di rinculo:
     • Rimescolamento dello stato iniziale (Initial-state reshuffling): Modifica i momenti dei partoni in entrata per conservare il quadrimomento mantenendo fissi i momenti dello stato finale.
     • Rimescolamento dello stato finale (Final-state reshuffling): Modifica i momenti in uscita (e potenzialmente l'energia del centro di massa partonico) per soddisfare le condizioni on-shell per le masse fisiche.
   • Questa strategia approssima l'inclusione di correzioni relativistiche di ordine superiore nella misura dello spazio delle fasi, migliorando la descrizione degli osservabili vicino alle soglie.

3. Modello UFO e Automazione:
   L'implementazione è racchiusa in un modello UFO (Universal Feynman Output) aggiornato, sm_onia. Questo modello definisce gli stati di Fock (inclusi le configurazioni a onda P come $^3P_J$ e $^1P_1$) e permette agli utenti di specificare le masse fisiche e gli LDME tramite schede di input (onia_card.dat). Il framework supporta sia la quarkonia (NRQCD) che la leptonia (NRQED), incluse le configurazioni color-singlet e color-octet.

Contributi Chiave e Risultati

• Validazione: L'implementazione P-wave è stata validata rispetto al generatore HELAC-Onia.

  • Ampiezze al Quadrato: I confronti a singoli punti dello spazio delle fasi mostrano un accordo al livello dell'aritmetica a doppia precisione ($10^{-14}$ a $10^{-16}$ di deviazione relativa).
  • Sezioni d'Urto Integrate: I benchmark per la produzione singola e doppia di quarkonia in collisioni $pp$ ed $e^+e^-$ mostrano un accordo entro le incertezze statistiche di Monte Carlo.
  • Leptonium: Per quanto riguarda la NRQED, il codice riproduce i calcoli analitici delle sezioni d'urto per la produzione di positronium e ditauonium a onda P a un livello sub-permille su un ampio intervallo di energie di collisione.

• Applicazioni Fenomenologiche:

  • Produzione Singola e Doppia di Quarkonium: Gli autori presentano le sezioni d'urto per vari processi all'LHC ($\sqrt{s}=13$ TeV) e ai B-factory ($\sqrt{s}=10.58$ GeV). Quantificano l'impatto del rimescolamento del momento, notando che, sebbene gli effetti siano spesso piccoli per la produzione singola inclusiva, possono essere significativi (fino al ~20-50%) per i canali di doppia produzione che coinvolgono stati con grandi differenze di massa (ad esempio, $J/\psi + \psi(2S)$).
  • Caso di Studio $J/\psi + \psi(2S)$: Uno studio dettagliato della produzione associata di $J/\psi + \psi(2S)$ all'LHC dimostra l'utilità del framework. Incorporando le masse fisiche e i contributi di feed-down P-wave (da $\psi(2S)$ e $\chi_{cJ}$), la sezione d'urto predetta viene ridotta rispetto ai trattamenti convenzionali. Questa riduzione porta a una migliore descrizione dei dati di LHCb, lasciando più spazio per potenziali contributi di Scattering di Doppio Partone (DPS). Lo studio evidenzia che i contributi P-wave rappresentano circa il 17% della sezione d'urto diretta e il 25% della sezione d'urto fiduciale totale tramite feed-down.
  • Produzione di Leptonium: Il framework consente lo studio di atomi esotici come il positronium e il ditauonium. Gli autori derivano previsioni NRQED in forma chiusa e confermano che la simulazione automatizzata corrisponde ai risultati analitici, validando i proiettori P-wave nel settore QED.

Significato e Rivendicazioni
L'articolo afferma che MadSONS costituisce un passo significativo verso la simulazione completamente automatizzata della produzione di stati legati. Completando l'automazione dei processi a onda S e P a Leading Order (LO), lo strumento permette l'inclusione sistematica di tutti gli stati di Fock fenomenologicamente rilevanti fino all'ordine $v^7$ in NRQCD.

Gli autori sottolineano che l'uso dei numeri duali fornisce una base robusta per future estensioni NLO, dove la cancellazione delle divergenze infrared richiede un'elevata precisione numerica. Inoltre, la procedura di rimescolamento del momento offre un metodo pratico per incorporare gli effetti di massa fisica senza richiedere correzioni relativistiche complete all'elemento di matrice, migliorando così l'affidabilità delle simulazioni vicino alle soglie cinematiche.

Questo lavoro posiziona MadSONS come uno strumento versatile per investigare i meccanismi di produzione di quarkonia, testare la fattorizzazione NRQCD ed esplorare sistemi esotici di leptonium, rimanendo al contempo compatibile con gli standard di strumenti di parton shower e hadronizzazione come Pythia8. Gli autori notano con modestia che, sebbene le previsioni LO migliorino la descrizione dei dati, potrebbero rimanere discrepanze dovute a correzioni QCD di ordine superiore mancanti, effetti relativistici e potenziali contributi DPS, che sono soggetti di indagine futura.