TreeLevelPhysique et méthodes
Ce que TreeLevel calcule exactement, comment, et contre quoi il a été vérifié.
Périmètre
TreeLevel calcule des amplitudes à l'arbre (aucune boucle), pour des états initiaux et finaux de particules élémentaires ou effectives du modèle : diffusions 2 → n et désintégrations 1 → n. Les résultats sont des sections efficaces en picobarns (partoniques quand il s'agit de quarks : pas de fonctions de distribution de partons), des largeurs partielles en GeV, des durées de vie en secondes (ħ = 6,582 × 10⁻²⁵ GeV·s), des rapports de branchement quand le modèle déclare la largeur totale. Les largeurs de résonance déclarées entrent dans les propagateurs (−i M Γ au dénominateur) ; les diagrammes d'un même processus sont sommés de façon cohérente, avec leurs signes relatifs.
Conventions
- Métrique (+, −, −, −) ; ε0123 = +1 ; γ⁵ = iγ⁰γ¹γ²γ³, PL,R = (1 ∓ γ⁵)/2 ; base de Weyl (chirale) pour les matrices γ, comme ALOHA.
- Règle de vertex = i × (terme du lagrangien, champs retirés), toutes les impulsions entrantes, ∂μ → −i pμ. Propagateurs : i/(p² − m² + iMΓ) pour les scalaires, i(p̸ + m)/(p² − m² + iMΓ) pour les fermions, −i gμν/(p² − M² + iMΓ) pour les vecteurs en jauge de Feynman, avec le terme −pμpν/M² en jauge unitaire.
- Modèle standard : conventions de FeynRules — Dμ = ∂μ − i gs TaGa − i (e/sW) (σi/2) Wi − i (e/cW) Y B, Φ = (−iG⁺, (v + h + iG⁰)/√2), jauge de Feynman–'t Hooft ; paramètres numériques de l'UFO
sm(α⁻¹ = 127,9, GF = 1,16637 × 10⁻⁵ GeV⁻², MZ = 91,1876 GeV, MW dérivée). - Couleur : Ta = λa/2, fabc standard, εijk pour les triplets. Les indices de couleur externes sont sommés, l'état initial moyenné (3 par quark, 8 par gluon).
Amplitudes numériques
Pour un point cinématique donné, chaque diagramme est un réseau de tenseurs : les règles des vertex (évaluées numériquement avec les impulsions des pattes), les propagateurs, et les deltas de couleur des lignes internes. Le réseau est contracté sur ses indices internes dans un ordre glouton (toujours le tenseur partageant le plus d'indices avec le résultat courant), ce qui évite de former de grands produits extérieurs — les vertex à quatre gluons y gagnent un facteur 25. Ce réseau ne dépend pas des hélicités : il est calculé une fois, puis contracté avec les fonctions d'onde externes (spineurs u, v, ū, v̄ d'hélicité ±½ ; vecteurs de polarisation εμ d'hélicité ±1, et 0 pour les bosons massifs) pour chaque configuration d'hélicités. Σ|ℳ|² est la somme sur ces configurations et sur les couleurs externes, divisée par la multiplicité de spin et de couleur de l'état initial.
Les diagrammes dont les fermions externes jouent les mêmes rôles (ligne/colonne) partagent leurs fonctions d'onde : leurs réseaux sont sommés avant la boucle sur les hélicités. Les intégrations sur l'angle ou l'énergie sont réparties sur tous les cœurs.
Fermions : chaînes, signes, Majorana
Chaque ligne fermionique continue d'un diagramme forme une chaîne lue du spineur ligne (ū sortant ou v̄ entrant) au spineur colonne (u entrant ou v sortant). Le signe relatif entre diagrammes est la parité de la permutation des fermions externes, et (−1) par boucle fermionique.
Avec des fermions de Majorana ou des couplages qui violent le nombre fermionique, une chaîne n'a plus de sens « naturel ». TreeLevel suit les règles de Denner, Eck, Hahn et Küblbeck (1992) : on oriente chaque chaîne arbitrairement ; un vertex traversé contre son orientation propre donne Γ′ = CΓᵀC⁻¹ (γμ change de signe, 1, γ⁵ et γμγ⁵ non), un propagateur remonté donne S(−p), et les spineurs externes dépendent de la position dans la chaîne et du caractère entrant/sortant — pas de particule/antiparticule. Les spineurs satisfont v = Cūᵀ exactement, condition sans laquelle les interférences entre diagrammes d'orientations différentes sont fausses. Le résultat est indépendant de l'orientation choisie, ce que la suite de tests vérifie en inversant systématiquement le point de départ des chaînes.
Espace des phases et observables
| Processus | Méthode | Résultat |
|---|---|---|
| 2 → 2 | dσ/dΩ = |ℳ|² |pf| / (64π² s |pi|) dans le centre de masse ; σ par quadrature de Gauss–Legendre en cos θ, coupure |cos θ| ≤ c optionnelle | ⟨|ℳ|²⟩(cos θ), dσ/dΩ, σ ; balayage en √s |
| 1 → 2 | formule fermée Γ = |p| |ℳ|² / (8π M²) | Γ, τ, T½, BR |
| 1 → 3 | intégration de Dalitz (deux masses invariantes), quadrature | Γ, τ, T½, BR |
| 2 → n, 1 → n | Monte Carlo RAMBO (points uniformes dans l'espace des phases massif), coupures en énergie, masse de paire et angle, erreur statistique | σ ± δσ, Γ ± δΓ |
Les particules identiques dans l'état final reçoivent le facteur 1/n! ; le facteur de symétrie d'un diagramme (1/|Aut|) est calculé sur son graphe. Conversion : 1 GeV⁻² = 0,389379 mb.
Moteur symbolique
Le même graphe est aussi traité symboliquement : chaque chaîne fermionique fermée par la conjugaison complexe devient une trace de Dirac, développée en produits scalaires pi·pj et masses (γ⁵ par le tenseur ε, contracté ensuite) ; les sommes de polarisation −gμν (+ pμpν/M² pour les bosons massifs) et les facteurs de couleur sont appliqués ; le résultat, une somme de fractions rationnelles regroupées par dénominateur, est écrit en s, t, u avec les substitutions demandées (masses nulles, u = Σm² − s − t, cW² = 1 − sW²) et rendu en LaTeX. Il coïncide avec le calcul numérique à 10⁻⁸ près sur tous les processus testés — c'est lui qui a permis de déceler, puis de corriger, une erreur d'interférence dans le traitement numérique des Majorana.
Validations
| Ce qui est vérifié | Résultat |
|---|---|
Règles du Modèle standard contre l'UFO sm de MadGraph | 129 vertex sur 129 identiques, CKM complexe comprise |
| QED : e⁺e⁻ → μ⁺μ⁻, Bhabha, Compton, e⁺e⁻ → γγ | formules des manuels (Peskin & Schroeder) ; Bhabha en trois termes |
| e⁺e⁻ → W⁺W⁻ | annulation de jauge γ + Z + ν ; σ(200 GeV) = 19,2 pb |
| Asymétrie avant-arrière au pic du Z | valeur attendue des couplages gV, gA |
| Largeurs Z → ff̄, W → ℓν, h → bb̄, ττ, WW*, ZZ* | formules fermées, rapports de branchement |
| Muon : τ = 2,19 µs ; π⁰ → γγ : 7,79 eV ; π → μν ; K → μν ; τ → πν, Kν : BR 10,6 %, 0,70 % ; Ds → μν 0,54 % ; Σ⁰ → Λγ 8,9 keV ; Λ → p e ν̄ à 8 % ; n → p e ν̄ : 946 s | mesures, à l'arbre |
| Monte Carlo RAMBO | volumes d'espace des phases exacts ; 2 → 2 et désintégration du muon retrouvés par MC |
| Majorana : Γ(Z → NN) = Γ(Z → νν̄) pour N sans masse ; Γ(N → ℓ∓W±) ; e⁻e⁻ → W⁻W⁻ | indépendance du flux, formules fermées |
| MSSM_SLHA2 importé : σ(e⁺e⁻ → χ̃⁰₁χ̃⁰₁, 500 GeV) | 269 fb à SPS1a (littérature : 250–300 fb) ; symbolique = numérique sur χ̃⁺χ̃⁻, χ̃⁰₁χ̃⁰₂, g̃g̃ |
| SU(5) : X → uu, e⁺d̄ ; Y → ud, e⁺ū, ν̄d̄ | Γ = g₅²M/(24π) ; BR ½ ½ et ½ ¼ ¼ ; ΓX = ΓY ; équivalence avec un champ conjugué indépendant (facteur 2 exact des chiralités) |
| Symbolique contre numérique | 10⁻⁸ sur une douzaine de processus, MS, MSSM et SU(5) |
La suite de tests du moteur (plus de cent cas) rejoue tout ceci à chaque modification.
Limites connues
- Arbre seulement : pas de boucles, pas de corrections radiatives, pas d'évolution des couplages. Les fantômes de Faddeev–Popov ne sont pas dérivés (inutiles à l'arbre).
- Pas de fonctions de distribution de partons ni d'hadronisation : les processus à quarks sont partoniques ; le modèle
Hadronsoffre une description effective de quelques hadrons légers. - Pas de vertex de contact à quatre fermions (chaque vertex fermionique a exactement deux fermions).
- Les largeurs entrent dans les propagateurs sans schéma de soustraction complexe ; hors résonance, l'effet est négligeable.
- Le modèle SU(5) intégré comporte les couplages fermioniques de X et Y, pas leurs auto-couplages de jauge.
- Monte Carlo : l'échantillonnage RAMBO est uniforme ; près des pôles (photons mous), le nombre de points doit augmenter — les coupures sont là pour cela.