A resonance-aware MC@NLO QCD+EW-matched calculation of lepton-pair production

Sintesi Tecnica: Un Calcolo QCD+EW con Matching MC@NLO Sensibile alla Risonanza per la Produzione di Coppie di Leptoni

Definizione del Problema
Mentre il campo della fisica degli acceleratori transita verso l'era dell'HL-LHC (High-Luminosity LHC), le incertezze sperimentali si stanno riducendo, richiedendo miglioramenti commensurati nella precisione teorica. Sebbene i calcoli al Next-to-Leading Order (NLO) in QCD siano lo standard, gli effetti Elettrodeboli (EW) stanno diventando sempre più critici, particolarmente per osservabili che coinvolgono particelle prive di colore o grandi isospin deboli. Inoltre, la radiazione QED nei collisionatori di adroni può indurre grandi correzioni a specifici osservabili.

Una sfida significativa sorge quando l'accoppiamento (matching) di calcoli NLO QCD+EW a ordine fisso con gli sciami di partoni (parton showers, PS) in presenza di risonanze interne (ad esempio, i bosoni $Z$ o $W$). Gli algoritmi standard di parton shower, come quelli basati sul formalismo dei dipoli di Catani-Seymour, si affidano all'assegnazione del rinculo a partoni spettatori per mantenere la conservazione del momento. Nei processi risonanti, questa assegnazione del rinculo può alterare la virtualità del propagatore risonante, introducendo termini di ordine superiore spurii che sono numericamente grandi ma teoricamente non controllati. Tali artefatti distorcono la forma della risonanza e compromettono la precisione richiesta per processi come la produzione di coppie di leptoni Drell-Yan ($pp \to \ell^+\ell^-$).

Metodologia
Gli autori presentano il primo accoppiamento automatizzato di NLO QCD+EW a uno sciame di partoni QCD+QED intersecato (interleaved) utilizzando il metodo MC@NLO all'interno del formalismo dei dipoli di Catani-Seymour, implementato nel generatore di eventi SHERPA. La metodologia comprende tre componenti principali:

1. Sciame di Dipoli QCD+QED Intersecato:
   Gli autori estendono il framework CSSHOWER per evolvere simultaneamente le emissioni QCD e QED. L'evoluzione dello sciame è definita da un funzionale generatore dove il parametro di ordinamento $t$ è il momento trasverso per le emissioni e la virtualità per gli splitting. Il fattore di Sudakov è costruito come prodotto di componenti QCD e QED, garantendo un'evoluzione intersecata piuttosto che una semplice somma. Per gestire le funzioni di splitting non definite positive inerenti ai dipoli QED (a causa delle correlazioni di carica), viene impiegato un approccio di sciame pesato, con strategie specifiche per minimizzare i pesi negativi limitando gli spettatori dei dipoli.

2. Matching MC@NLO con Modifica di Sudakov:
   L'equazione maestra standard di MC@NLO è adattata per accoppiare lo sciame intersecato a elementi di matrice NLO QCD+EW. Ciò comporta la generazione di "S-events" (Born + virtuale + dipoli integrati) ed "H-events" (correzioni reali).

   • Modifica degli H-event: Gli autori identificano che negli sciami intersecati, la dinamica QCD può indurre grandi soppressioni di Sudakov ($\bar{\Delta} \ll 1$) che, combinate con le emissioni reali QED negli H-events, portano a contributi negativi non fisici nella sezione d'urto differenziale. Per mitigare questo, introducono un fattore di Sudakov aggiuntivo, $\Delta_n(t_n, t_{n+1})$, calcolato tramite lo standard sciame intersecato, applicato direttamente agli H-events. Questo approccio "MC@NLO$_{sud}$" sopprime i termini di ordine superiore non fisici ($O(\alpha_s \alpha)$) mantenendo al contempo l'accuratezza formale NLO.

3. Sottrazione e Evoluzione Sensibili alla Risonanza:
   Per affrontare il problema della risonanza, gli autori sviluppano una modifica della sottrazione e dello sciame sensibile alla risonanza. Lo spazio delle fasi è diviso in regimi "risonanti" e "non risonanti" basandosi su due criteri:

   • Virtualità ($\Delta_r$): La distanza della virtualità del propagatore risonante dalla sua massa nominale, normalizzata rispetto alla larghezza $\Gamma_r$.
   • Scala di Emissione ($t$): La durezza dell'emissione rispetto alla larghezza della risonanza al quadrato ($\Gamma_r^2$).

   Nel regime risonante (emissioni soft/collineari rispetto alla larghezza), i dipoli standard che attraversano la risonanza sono attivi. Nel regime non risonante (emissioni dure che risolvono la risonanza), il processo si fattorizza in sottoprocessi di produzione e decadimento. I dipoli che collegano i settori di produzione e decadimento sono soppressi o disattivati mediante funzioni $\Theta$. Ciò assicura che le assegnazioni del rinculo non distorcano artificialmente la distribuzione della virtualità del propagatore risonante. Per le risonanze neutre (come il bosone $Z$), ciò comporta la modifica dei dipoli final-initial e initial-final; per le risonanze cariche, sarebbero richiesti ulteriori dipoli risonanti (notato come lavoro futuro).

Risultati Chiave
Gli autori validano il loro framework utilizzando $pp \to \nu_e \bar{\nu}_e$ (non risonante) e $pp \to e^+ e^-$ (risonante) a $\sqrt{s} = 13.6$ TeV.

• Validazione contro l'Ordine Fisso: Per i processi non risonanti, i risultati MC@NLO QCD+EW corrispondono da vicino ai calcoli a ordine fisso NLO QCD+EW per le emissioni dure, con gli S-events che dominano le regioni soft ed esibiscono la corretta soppressione di Sudakov.
• Correzione degli Artefatti: Nel canale $pp \to e^+ e^-$, il calcolo standard MC@NLO QCD+EW (senza la modifica di Sudakov per gli H-events) produce sezioni d'urto negative non fisiche e "kink" negli spettri del momento trasverso del fotone ($p_{\perp, \gamma}$) e della massa dei dileptoni ($m_{\ell\ell}$). L'inclusione del fattore di Sudakov aggiuntivo in MC@NLO$_{sud}$ rimuove questi artefatti, ripristinando il comportamento fisico e l'accordo con le correzioni QED riassunte tramite YFS.
• Impatto della Sensibilità alla Risonanza: Confrontando i calcoli "resonance-aware" (sensibili alla risonanza) e quelli standard ("resonance-unaware") per $pp \to e^+ e^-$ emerge che:
  • Massa dei Dileptoni: Una distorsione fino allo 0.5% nella distribuzione della massa invariante sotto il picco di $Z$ ($m_{\ell\ell} \approx 60$ GeV) quando si utilizzano i dipoli standard. Il trattamento sensibile alla risonanza corregge questo aspetto, garantendo che la forma sia indipendente dai parametri dello schema di sottrazione.
  • Momento Trasverso del Fotone: Differenze significative nella coda del fotone duro. Il calcolo sensibile alla risonanza potenzia le emissioni di fotoni più duri fino al 15% a $p_{\perp, \gamma} \approx 1$ TeV rispetto all'approccio standard, grazie al corretto trattamento degli effetti di interferenza ISR-FSR tramite i dipoli che attraversano la risonanza.
• Potenziamento Statistico: Il paper investiga l'incremento artificiale degli splitting QED per migliorare la statistica. Sebbene un fattore $c=5$ migliori il campionamento dei fotoni senza artefatti, $c=10$ introduce deviazioni non fisiche sub-percentuali a causa delle fluttuazioni dei pesi, portando gli autori a raccomandare di non procedere con un potenziamento aggressivo nel primo emissione accoppiato.

Significatività e Rivendicazioni
Il paper rivendica di presentare il primo accoppiamento automatizzato di NLO QCD+EW a uno sciame di partoni QCD+QED intersecato utilizzando il metodo MC@NLO nel formalismo dei dipoli di Catani-Seymour. La sua importanza primaria risiede in:

1. Precisione per l'HL-LHC: Fornire uno strumento teorico capace di accoppiare la crescente precisione dei dati sperimentali includendo coerentemente sia le correzioni di ordine superiore QCD che EW.
2. Gestione della Risonanza: Risolvere il problema specifico dei termini di ordine superiore spurii introdotti dalle assegnazioni del rinculo nei processi risonanti. La modifica proposta, sensibile alla risonanza, assicura che la distribuzione della virtualità dei propagatori risonanti non sia distorta dallo sciame di partoni, un aspetto vitale per le misure di precisione del bosone $Z$ e l'estrazione delle PDF.
3. Rimozione degli Artefatti: Dimostrare che i pesi negativi non fisici e le distorsioni spettrali nell'accoppiamento intersecato QCD+EW possono essere curati con una modifica di Sudakov fisicamente motivata applicata agli H-events.

Gli autori notano che, sebbene l'algoritmo sensibile alla risonanza sia generale, la loro attuale implementazione è specifica per il processo (focalizzata sul Drell-Yan a corrente neutra). Essi riconoscono che l'estensione a risonanze con colore o carica (come i quark top o i bosoni $W$) richiede funzioni di splitting aggiuntive in cui la risonanza agisce sia come emittente che come spettatore, un compito riservato a lavori futuri.