Génération automatisée d'événements pour la production de quarkonium et leptonium en onde S dans NRQCD et NRQED

Informations fondamentales

• ID de l'article: 2510.26773
• Titre: Automated event generation for S-wave quarkonium and leptonium production in NRQCD and NRQED
• Auteurs: Alice Colpani Serri, Chris A. Flett, Jean-Philippe Lansberg, Olivier Mattelaer, Hua-Sheng Shao, Lukas Simon
• Classification: hep-ph (physique des hautes énergies - phénoménologie), hep-ex (physique des hautes énergies - expérience), nucl-ex (physique nucléaire - expérience), nucl-th (physique nucléaire - théorie)
• Date de soumission: 30 octobre 2025
• Lien de l'article: https://arxiv.org/abs/2510.26773

Résumé

Cet article propose une extension du cadre MadGraph5_aMC@NLO capable de calculer automatiquement les sections efficaces à l'ordre dominant pour la production de quarkonium en onde S (quarkonium) et de leptonium en onde S dans le formalisme de factorisation NRQCD (QCD non-relativiste) et NRQED (QED non-relativiste). Le cadre a été validé par plusieurs processus de référence, démontrant sa robustesse et sa flexibilité pour les études phénoménologiques. L'avantage clé de cette implémentation réside dans son intégration transparente aux fonctionnalités existantes de MadGraph5_aMC@NLO, supportant non seulement les calculs du Modèle Standard, mais aussi divers scénarios au-delà du Modèle Standard ou de théorie effective des champs via une interface UFO (Universal Feynman Output) modifiée. De plus, ce cadre reste compatible avec les fonctionnalités de douche partonique et d'hadronisation des générateurs Monte Carlo standard. À travers de nombreux exemples, les auteurs soulignent que l'étude théorique des processus de quarkonium nécessite une attention particulière aux contributions sous-dominantes, dont les effets sont souvent difficiles à prédire uniquement par des arguments de comptage simples basés sur les règles de hiérarchie de couplage et d'échelle de vitesse.

Contexte et motivation de la recherche

Problème de recherche

Cette recherche vise à résoudre le problème de la génération automatisée d'événements pour les processus de production de quarkonium et de leptonium, en particulier dans le cadre théorique NRQCD et NRQED.

Importance du problème

1. Signification physique: Le quarkonium (comme J/ψ, Υ) constitue une sonde importante pour étudier les aspects perturbatifs et non-perturbatifs de la QCD. Depuis la découverte du J/ψ il y a 50 ans, qui a déclenché la "révolution de novembre", le mécanisme de production du quarkonium reste incompletement compris.
2. Besoins expérimentaux: Les collisionneurs comme le LHC et l'EIC produisent de grandes quantités de données sur le quarkonium, nécessitant des outils de prédiction théorique fiables pour la comparaison.
3. Exploration de nouvelle physique: Le leptonium (comme le positronium, le vrai muonium, le ditauonium) peut être utilisé pour tester précisément l'QED, chercher une nouvelle physique et mesurer les constantes fondamentales.

Limitations des méthodes existantes

1. MadOnia: Basé sur MadGraph/MadEvent v4, supporte uniquement un seul quarkonium en onde S ou P, n'a pas été porté au cadre moderne MG5_aMC.
2. HELAC-Onia: Supporte plusieurs quarkonium en onde S et P, mais nécessite une combinaison manuelle des coefficients de courte portée lors du traitement des états liés physiques.
3. Générateurs MC universels: Les fonctionnalités d'éléments de matrice et de douche partonique pour le quarkonium dans Pythia, Herwig et Sherpa ne sont pas encore complètement compatibles.
4. Absence d'outils pour le leptonium: Les outils théoriques disponibles publiquement pour la production de leptonium sont extrêmement limités.

Motivation de la recherche

Développer un outil complet, automatisé et convivial capable de:

• Traiter un nombre arbitraire de quarkonium et leptonium en onde S dans le cadre NRQCD/NRQED
• S'intégrer transparemment au cadre largement utilisé MG5_aMC
• Supporter les scénarios du Modèle Standard et au-delà
• Jeter les bases pour les calculs de précision NLO futurs

Contributions principales

1. Première implémentation: Première implémentation de la génération automatisée d'événements pour le quarkonium et leptonium en onde S dans le cadre MG5_aMC
2. Cadre unifié: Fournit une implémentation NRQCD/NRQED unifiée pour traiter le quarkonium et le leptonium
3. Nouveau modèle UFO: Développe le modèle sm_onia incluant un fichier boundstates.py spécialisé définissant le contenu des états de Fock
4. Validation étendue: Valide l'implémentation par comparaison détaillée avec HELAC-Onia et les expressions analytiques
5. Applications riches: Démontre diverses processus de production sur collisionneurs pp, e⁺e⁻, e⁻p, incluant:
   • Production de quarkonium unique et production associée avec bosons électrofaibles et jets
   • Production de double/triple quarkonium
   • Production de leptonium (positronium, vrai muonium, ditauonium)
6. Perspectives théoriques: Illustre par de nombreux exemples que l'importance des contributions sous-dominantes ne peut pas être simplement estimée par les règles d'échelle de vitesse

Détails méthodologiques

Cadre théorique

Factorisation NRQCD/NRQED

La section efficace différentielle d'un état lié B peut s'exprimer comme: $d\sigma(N_1 + N_2 \to B + X) = \sum_{I_1,I_2,n} \int dx_1 dx_2 f_{I_1/N_1}(x_1) f_{I_2/N_2}(x_2) d\hat{\sigma}(I_1I_2 \to (C_1C_2)[n] + X_p) \langle \mathcal{O}^B_n \rangle$

où:

• $(C_1C_2)$: paire de quarks lourds $(Q\bar{Q}')$ ou paire de leptons $(\ell^-\ell'^+)$
• $n = {}^{2S+1}L_J^{[C]}$: état quantique (spin S, moment angulaire orbital L, moment angulaire total J, état de couleur C)
• $d\hat{\sigma}$: section efficace de courte portée (calculable en QCD/QED perturbatif)
• $\langle \mathcal{O}^B_n \rangle$: élément de matrice de longue portée (LDME)

Méthode de projection

Projection des amplitudes sur des états de Fock spécifiques via une méthode de projection covariante:

1. Projection de couleur (uniquement pour le quarkonium):
   • Singulet de couleur: $P_{C=1} = \delta_{c_3c_4}/\sqrt{N_c}$
   • Octet de couleur: $P_{C=8} = \sqrt{2} t^{c_{34}}_{c_4c_3}$
2. Projection de spin: $P_S = \frac{\bar{v}_{\lambda_{C_2}}(k_4) \Gamma_S u_{\lambda_{C_1}}(k_3)}{2\sqrt{2m_{C_1}m_{C_2}}}$ où $\Gamma_{S=0} = \gamma_5$ (singulet), $\Gamma_{S=1} = \not{\varepsilon}^*_{\lambda_s}(K)$ (triplet)
3. Projection du moment angulaire orbital (onde S: L=0): $\mathcal{A}^{(n,0)}_{\{[C],S,L=0\}}(r) = \left[\mathcal{A}^{(n,0)}_{\{[C],S\}}(r)\right]_{q=0}$

Architecture d'implémentation

Extension du modèle UFO

1. Nouvelle classe Boundstate: Similaire à la classe Particle, définit les états de Fock de l'état lié
2. Convention de nommage: Jpsi(1|3S11) représente J/ψ avec nombre quantique principal N=1, triplet de spin ³S₁, singulet de couleur
3. Codage PDG:
   • Singulet de couleur: utilise le codage PDG standard (ex: J/ψ: 443)
   • Octet de couleur: codage non-standard 99nqnsnrnLnJ
   • Leptonium: nouveau codage ±990nℓ₁nℓ₂nrnLnJ

Flux de génération de code

1. Projection de couleur: Appliquée au niveau Python lors de la génération des diagrammes de Feynman
2. Projection de spin: Génère des sous-programmes dédiés via le module ALOHA, dépendant de la bibliothèque HELAS
3. Traitement LDME: Intégré directement dans le calcul de l'amplitude au carré, via entrée onia_card.dat
4. Intégration d'espace de phase: Adopte une stratégie de pseudo-multi-canal, utilisant l'intégrateur VEGAS avec parallélisation

Modes d'exécution

1. Mode ordre dominant: Simulation complète d'événements LO, génération d'événements au format LHEF
2. Mode autonome: Génère uniquement le code de calcul d'éléments de matrice (Fortran), pour intégration dans des cadres externes

Points d'innovation technique

1. Combinaison automatique d'états de Fock: L'utilisateur spécifie l'état lié physique (ex: Jpsi), le code inclut automatiquement toutes les contributions d'états de Fock pertinents
2. Réglage flexible de LDME: Les valeurs LDME peuvent être ajustées à l'exécution sans recompilation
3. Compatibilité BSM: Supporte via l'interface UFO tout modèle de nouvelle physique (exemple HEFT fourni)
4. Interface de douche partonique: Intégration transparente avec les générateurs MC comme Pythia
5. Support du leptonium: Première capacité de calcul automatisé pour le leptonium

Configuration expérimentale

Configuration de calcul

Paramètres fondamentaux

• Masses de quarks: $m_c = 1.55$ GeV (3FS), $m_b = 4.7$ GeV, $m_t = 173$ GeV
• Masses de leptons: $m_e = 0.511$ MeV, $m_\mu = 105.66$ MeV, $m_\tau = 1.777$ GeV
• Paramètres électrofaibles: $m_Z = 91.188$ GeV, $m_H = 125$ GeV, $\alpha^{-1}_{G_\mu} = 132.507$
• Échelles de renormalisation/factorisation: $\mu_R = \mu_F = H_T/2 = \frac{1}{2}\sum_i \sqrt{k^2_{i,T} + m^2_i}$

Configuration des collisionneurs

1. Collision proton-proton (LHC): $\sqrt{s} = 13$ TeV, PDF4LHC21_40 NNLO PDF
2. Collision électron-positron (Belle/BaBar): $\sqrt{s} = 10.58$ GeV (résonance Υ(4S))
3. Collision électron-proton (EIC): $\sqrt{s} = 140.7$ GeV ($E_{e^-} = 18$ GeV, $E_p = 275$ GeV)

Valeurs LDME

Utilise les valeurs par défaut du tableau 2:

• Singulet de couleur: Calculé via la fonction d'onde $R_Q(0)$, ex: $\langle \mathcal{O}^{J/\psi}_{{}^3S_1^{[1]}} \rangle = 1.16$ GeV³
• Octet de couleur: Ajusté à partir des données expérimentales, ex: $\langle \mathcal{O}^{J/\psi}_{{}^1S_0^{[8]}} \rangle = 0.0146$ GeV³

Stratégie de validation

Validation au niveau des éléments de matrice

• Comparaison avec HELAC-Onia pour 42 processus différents d'amplitudes au carré
• Écart relatif $\Delta_{\text{rel}} < 10^{-6}$ (niveau de précision en virgule flottante)

Validation au niveau des sections efficaces

• Comparaison des sections efficaces totales et différentielles avec HELAC-Onia
• Valeurs de pull $< 3$, cohérentes avec les attentes statistiques

Validation par expressions analytiques

• Comparaison des processus de leptonium avec les formules analytiques de la littérature
• Accord à la précision du millième dans les limites de découplage du boson Z et H

Résultats expérimentaux

Présentation des résultats principaux

1. Production de quarkonium unique (collision pp, √s = 13 TeV)

Sections efficaces totales du processus 2→1:

• $\sigma(pp \to \eta_c + X) = 2.9366(5)$ μb
• $\sigma(pp \to J/\psi + X) = 536.14(6)$ nb
• $\sigma(pp \to \eta_b + X) = 5.4935(7)$ μb
• $\sigma(pp \to \Upsilon + X) = 6.0655(4)$ nb

Découvertes clés:

• La section efficace de l'état pseudoscalaire est bien supérieure à celle de l'état vecteur, car le processus QCD 2→1 de l'état vecteur est interdit par la conservation de la parité C
• La section efficace de $\eta_c$ est extrêmement sensible à la région de faible x du PDF, différents ensembles PDF pouvant inverser la taille relative par rapport à $\eta_b$

2. Production associée (collision pp)

Association avec bosons électrofaibles:

• $\sigma(pp \to J/\psi + Z) = 1.3425(4)$ pb (dominée par octet de couleur)
• $\sigma(pp \to \Upsilon + Z) = 340.7(1)$ fb (dominée par singulet de couleur)

Association avec jets/photons (avec coupures cinématiques appliquées):

• $\sigma(pp \to J/\psi + j) = 329.8(2)$ nb ($p_{T,j} > 10$ GeV)
• $\sigma(pp \to J/\psi + \gamma) = 19.13(1)$ nb ($p_{T,\gamma} > 2$ GeV)

3. Production de double quarkonium

Paires de même quarkonium:

• $\sigma(pp \to J/\psi + J/\psi) = 10.756(3)$ nb (dominée par singulet de couleur)
• $\sigma(pp \to \Upsilon + \Upsilon) = 44.63(1)$ pb

Production mixte:

• $\sigma(pp \to J/\psi + \Upsilon) = 10.503(3)$ pb (le processus QCD de singulet de couleur disparaît à $\mathcal{O}(\alpha^4_s)$)

Production de paires $B_c$:

• $\sigma(pp \to B^+_c + B^-_c) = 472.5(1)$ pb
• $\sigma(pp \to B^{\ast+}_c + B^{\ast-}_c) = 940.9(2)$ pb

4. Production de triple J/ψ

$\sigma(pp \to J/\psi + J/\psi + J/\psi) = 1.038(3)$ pb

La contribution de triple singulet de couleur n'est que de 54.57 fb, car $gg \to 3J/\psi[^3S_1^{[1]}]$ est interdit par la conservation de la parité C, la contribution principale provenant de l'état initial $q\bar{q}$.

Résultats de production de leptonium

Positronium (Ps)

Collision pp ($\sqrt{s} = 13$ TeV):

• $\sigma(pp \to \text{Ps}_1 + j) = 9.483(2)$ fb

Collision e⁺e⁻ ($\sqrt{s} = 10.58$ GeV):

• $\sigma(e^+e^- \to \text{Ps}_0 + \gamma) = 3.300(1)$ ab
• $\sigma(e^+e^- \to \text{Ps}_1 + \gamma) = 8.420(2)$ ab

Vrai muonium (TM)

Collision pp:

• $\sigma(pp \to \text{TM}_1 + j) = 9.460(2)$ fb (pratiquement identique au positronium)

Collision e⁺e⁻:

• $\sigma(e^+e^- \to \text{TM}_0 + \gamma) = 438.16(1)$ yb (extrêmement petit, supprimé par $4m^2_\mu/s$)
• $\sigma(e^+e^- \to \text{TM}_1 + \gamma) = 8.422(2)$ ab

Ditauonium (T)

Collision pp:

• $\sigma(pp \to \text{T}_1) = 9.1906(4)$ fb (production directe, première prédiction fiable)
• $\sigma(pp \to \text{T}_1 + j) = 6.055(1)$ fb
• $\sigma(pp \xrightarrow{\gamma\gamma} p\text{T}_0p) = 108.055(7)$ ab (processus UPC)

Collision e⁺e⁻:

• $\sigma(e^+e^- \to \text{T}_0 + \gamma) = 110.001(3)$ zb
• $\sigma(e^+e^- \to \text{T}_1 + \gamma) = 9.576(2)$ ab

Cas particulier: Production J/ψ + H (cadre HEFT)

Dans la théorie effective du Higgs, via le sommet effectif $Hgg$:

• $\sigma(pp \to J/\psi + H) = 1.53^{+0.40}_{-0.29}$ fb (contribution d'octet de couleur uniquement)
• $\sigma(pp \to \Upsilon + H) = 91^{+23}_{-17}$ ab

Caractéristiques de la distribution de rapidité:

• L'état ${}^1S_0^{[8]}$ présente un creux dans la région de rapidité centrale (pic à $|y| \approx 4$)
• L'état ${}^3S_1^{[8]}$ présente un pic dans la région de rapidité centrale
• Le creux provient du diagramme de sommet effectif $Hgg$, partiellement annulé par le diagramme de couplage de Yukawa

Effets de douche partonique

Pour le processus $pp \to J/\psi + c\bar{c}$, comparaison entre ordre fixe (fLO) et appariement avec douche partonique (LO+PS):

• ${}^3S_1^{[1]}$ et ${}^1S_0^{[8]}$: PS élargit la distribution $p_T$
• ${}^3S_1^{[8]}$: PS rétrécit la distribution $p_T$
• Effet total: annulation partielle, conduisant à une distorsion modérée

Analyse d'ablation: Contributions d'états de Fock

Pour $pp \to \eta_c$ comme exemple, décomposition en 6 canaux de partons:

• $gg \to \eta_c[^1S_0^{[1]}]$: 2.2995 μb (dominant)
• $gg \to \eta_c[^1S_0^{[8]}]$: 33.43 nb
• $q\bar{q} \to \eta_c[^3S_1^{[8]}]$: 603.7 nb
• Autres canaux: $\sim 0$ (interdits par conservation des nombres quantiques)

Perspectives importantes: Les estimations simples d'échelle de vitesse échouent souvent, car:

1. La conservation des nombres quantiques fait disparaître certains canaux
2. Les effets d'amplification de propagateur (ex: $g^* \to (Q\bar{Q})[^3S_1^{[8]}]$)
3. Les suppressions ou amplifications cinématiques/dynamiques

Travaux connexes

Générateurs MC universels

• Pythia v8.3/Herwig v7.4: Introduisent la douche partonique spécifique au quarkonium, mais les fonctionnalités d'éléments de matrice ne sont pas encore complètement compatibles
• Limitations: La production de quarkonium passe principalement par l'hadronisation ou la désintégration, l'implémentation d'éléments de matrice étant limitée

Générateurs d'événements au niveau des partons indépendants du processus

• MadOnia: Basé sur MadGraph v4, supporte un seul quarkonium S/P, n'a pas été porté à MG5_aMC
• HELAC-Onia: Supporte plusieurs quarkonium S/P, mais nécessite une combinaison manuelle des états de Fock, les LDME ne sont pas ajustables à l'exécution
• Avantages du présent travail:
  1. Cadre MG5_aMC moderne avec maintenance active
  2. Combinaison automatique d'états de Fock
  3. LDME ajustables à l'exécution
  4. Premier support du leptonium

Générateurs spécifiques au processus

• SuperChic/STARlight/eSTARlight: Focalisés sur les processus exclusifs
• EPOS4: Quarkonium dans les collisions d'ions lourds, utilisant la forme de matrice de densité de Wigner
• BCVEGPY: Spécialisé pour la production hadronique du mésons $B_c$
• Avantages du présent travail: Généralité, support de toute combinaison d'états de Fock en onde S

Calculateurs de sections efficaces

• FDCHQHP: Calculs NLO pour la production hadronique de quarkonium S/P
• FONLL: Calculs NLO+NLL pour J/ψ non-prompt
• Positionnement du présent travail: Automatisation LO, fondation pour extension NLO

Kits de calcul symbolique

• FeynOnium/AmpRed: Calcul symbolique semi-automatique
• Distinction du présent travail: Calcul numérique et génération d'événements complètement automatisés

Conclusions et discussion

Conclusions principales

1. Réussite technique: Implémentation réussie de la génération automatisée d'événements pour le quarkonium et leptonium en onde S dans le cadre MG5_aMC
2. Validation suffisante: Comparaison détaillée avec HELAC-Onia et expressions analytiques confirmant la correction de l'implémentation
3. Applications larges: Couvrant divers processus de production sur collisionneurs pp, e⁺e⁻, e⁻p
4. Perspectives théoriques: Illustration par de nombreux exemples que l'importance des contributions sous-dominantes ne peut pas être simplement estimée par les règles d'échelle de vitesse, nécessitant une analyse spécifique au processus
5. Leptonium: Première prédiction systématique de production de leptonium, guidant les observations expérimentales futures

Limitations

1. Précision LO: Implémentation actuelle uniquement à l'ordre dominant, les corrections NLO pouvant être significatives (ex: production de triple J/ψ)
2. Restriction aux ondes S: Support uniquement des états en onde S, l'implémentation des états P laissée pour travaux futurs
3. Absence de DPS: N'inclut pas les contributions de double/triple diffusion de partons, pouvant dominer certains processus
4. Effets de feed-forward: Ne considère pas les contributions de feed-forward des désintégrations d'états excités
5. Intégration d'espace de phase: Adopte une stratégie de pseudo-multi-canal simplifiée, l'efficacité pouvant être insuffisante pour processus extrêmement complexes
6. Région basse énergie: Pour processus très basse énergie comme $pp \to \text{Ps}_1$, l'hypothèse de factorisation peut échouer

Directions futures

1. Implémentation NLO: Réaliser la précision NLO basée sur le schéma de soustraction FKS étendu (fondation théorique existante 67)
2. Support des ondes P: Implémenter les opérateurs de projection de moment angulaire orbital et total
3. Optimisation multi-canal: Développer une paramétrisation d'espace de phase efficace spécifique aux états liés
4. Ajustement global: Combiner avec données expérimentales pour ajustement global des LDME
5. Intégration DPS: Combiner avec cadre DPS pour prédictions complètes
6. Collisions d'ions lourds: Étendre aux PDF nucléaires et effets de matière nucléaire froide

Impact physique

1. Guidance expérimentale: Fournir des prédictions fiables pour recherches de quarkonium/leptonium au LHC, EIC, Belle II, etc.
2. Compréhension des mécanismes: Aider à distinguer les mécanismes de singulet/octet de couleur, contributions SPS/DPS
3. Exploration de nouvelle physique: Supporter l'étude de production d'états liés dans scénarios BSM
4. Tests de précision: Fournir outils théoriques pour tests de précision QCD/QED

Évaluation approfondie

Points forts

Innovativité technique

1. Conception architecturale excellente: Introduction élégante de la classe Boundstate et modification du format UFO pour intégrer les états liés au cadre MG5_aMC
2. Convivialité: Convention de nommage claire (ex: Jpsi(1|3S11)), LDME ajustables à l'exécution, courbe d'apprentissage douce
3. Implémentation modulaire: Séparation des projections de couleur/spin/moment angulaire orbital, facilitant extension et maintenance
4. Validation multi-niveaux: Validation à trois niveaux - éléments de matrice, sections efficaces, formules analytiques - assurant fiabilité

Contributions scientifiques

1. Caractère novateur: Première implémentation du calcul automatisé du leptonium dans cadre MC principal
2. Systématicité: Couvrant 3 types de collisionneurs, 40+ processus, démontrant applicabilité large
3. Perspectives profondes: Illustration par nombreux exemples des limitations des règles d'échelle de vitesse, corrigeant malentendus possibles
4. Valeur pratique: Fournissant aux physiciens expérimentaux outil facile d'utilisation, réduisant barrière d'entrée pour calculs théoriques

Qualité de l'article

1. Structure claire: Chaîne logique théorie-implémentation-validation-applications complète
2. Détails suffisants: Fournissant configuration paramétrique complète et exemples de commandes reproductibles
3. Richesse graphique: 68 pages incluant 5 figures et 10 tableaux, densité informationnelle élevée
4. Littérature complète: 297 références, synthèse suffisante des travaux connexes

Insuffisances et suggestions d'amélioration

Limitations méthodologiques

1. Efficacité d'espace de phase: Stratégie de pseudo-multi-canal limitée pour processus haute multiplicité, suggérant:
   • Développer échantillonnage d'importance spécifique aux états liés
   • Utiliser apprentissage machine pour optimiser paramétrisation d'espace de phase
2. Comportement basse énergie: Hypothèse de factorisation pour leptonium très basse énergie nécessite justification plus soignée
3. Incertitude LDME: Utilisation de valeurs LDME fixes sans évaluation systématique d'impact d'incertitude

Configuration expérimentale

1. Choix de coupures: Certaines coupures cinématiques (ex: $p_{T,j} > 2$ GeV) plutôt arbitraires, manquant étude systématique
2. Dépendance PDF: Utilisation d'ensemble PDF unique, sans évaluation d'incertitude PDF
3. Dépendance d'échelle: Certains processus manquant évaluation d'incertitude de variation d'échelle

Discussion physique

1. Distinction de mécanismes: Bien que montrant contributions différents états de Fock, manquant discussion approfondie sur distinction expérimentale
2. Observables: Focalisé principalement sur sections efficaces totales et distributions single-particule, manquant polarisation, corrélations
3. Estimation de bruit de fond: Pas de discussion sur bruit de fond expérimental et efficacité de détection

Évaluation d'impact

Impact académique

• Potentiel de citation: Élevé, comme première implémentation MG5_aMC d'états liés, deviendra outil standard du domaine
• Recherches ultérieures: Attendu catalyser nombreuses études d'application (analyses données LHC, prédictions physique EIC, etc.)
• Développement théorique: Fondation pour automatisation NLO, intégration DPS, travaux ultérieurs

Valeur pratique

• Application expérimentale: Expériences ATLAS, CMS, LHCb, Belle II, etc. peuvent utiliser directement
• Valeur éducative: Réduisant barrière d'entrée pour recherche physique états liés, favorisant formation nouveaux chercheurs
• Contribution open-source: Code sera publié via Launchpad et NLOAccess, promouvant développement communautaire

Reproductibilité

• Disponibilité code: Engagement de publication, matériaux supplémentaires fournis
• Documentation détaillée: Article incluant exemples utilisation complets et configuration paramétrique
• Validation suffisante: Vérifications croisées multiples assurant fiabilité résultats

Scénarios d'application

Applications idéales

1. Physique LHC: Prédictions section efficace production quarkonium, études mécanismes, recherche nouvelle physique
2. Physique EIC: Diffusion profondément inélastique avec production quarkonium, extraction TMD gluon
3. Usines B: Tests QED précision, recherche leptonium
4. Recherche théorique: Ajustement LDME, vérification mécanismes octet de couleur, études DPS

Scénarios non-applicables

1. Collisions d'ions lourds: Manquant effets matière nucléaire thermique/froide, nécessitant combinaison avec outils EPOS4
2. Prédictions ultra-précision: Précision LO insuffisante, attendant implémentation NLO
3. Processus dominés onde P: Ex: production χc, attendant implémentation onde P
4. Processus exclusifs: Ex: production centrale exclusive, nécessitant outils spécialisés SuperChic

Recommandations et perspectives

Pour auteurs

1. Implémenter rapidement précision NLO, besoin le plus urgent
2. Développer support onde P, complétant couverture états de Fock
3. Fournir cadre ajustement global LDME
4. Renforcer collaboration avec groupes expérimentaux, promouvant application outil

Pour utilisateurs

1. Comprendre limitations précision LO, interpréter prudemment prédictions absolues
2. Effectuer évaluation systématique incertitudes (échelle, PDF, LDME)
3. Combiner avec autres outils (ex: générateurs DPS) pour image complète
4. Fournir retours utilisation, promouvant amélioration outil

Pour domaine

1. Établir base de données LDME standardisée
2. Développer techniques automatisation NLO physique états liés
3. Promouvoir mesures expérimentales couvrant plus processus et observables
4. Renforcer collaboration théorie-expérience, réalisant tests QCD/QED précision

Références (sélection)

Fondations théoriques principales

• 9 Bodwin, Braaten, Lepage (1995): Travail fondateur formalisme factorisation NRQCD
• 64 Caswell, Lepage (1986): Lagrangienne effective NRQED
• 67 A A H, Shao, Simon (2024): Soustraction FKS production quarkonium NLO

Outils connexes

• 25 Artoisenet et al. (2008): MadOnia - outil automatisation NRQCD génération précédente
• 27,28 Shao (2013, 2016): HELAC-Onia - référence comparaison principale
• 6,7 Alwall et al. (2014, 2018): Cadre MadGraph5_aMC@NLO

Applications importantes

• 112 Lansberg, Shao (2013): Importance corrections α⁵s production paires J/ψ
• 194 Shao, Zhang (2019): Production triple J/ψ comme sonde diffusion multi-partons
• 4 d'Enterria, Shao (2023): Perspectives découverte ditauonium

Évaluation globale: Ceci est un article de haute qualité en méthodologie, avec implémentation technique solide, validation suffisante, et applications larges. Il comble un vide important dans les outils automatisés pour physique états liés, particulièrement pour domaine leptonium. Bien que présentant limitations de précision LO et restriction aux ondes S, comme première étape d'une série de travaux, il a jeté fondations solides pour développements futurs. Recommandé aux chercheurs expérimentaux et théoriques de surveiller étroitement développements ultérieurs, particulièrement implémentation NLO et extension onde P.