Diffusion et interaction entre

particules
 Rappels sur la relativité restreinte

 Cinématique d’une réaction

 Matrice de diffusion S – section efficace

 Largeur de désintégration

1
 Relativité & Formalisme
quadridimensionnel
• Vecteur de position espace-temps : x µ = ( x 0 , x1 x 2 x 3 ) = (t , x ) ((ct , x ))
composantes contravariantes
3
• Espace à 4 dimensions : A = ∑ A µ eµ =A µ eµ
µ =0

A ⋅ B = Aµ eµ ⋅ Bυ eυ = A µ Bυ g µυ
gμν : tenseur métrique
g µυ = eµ ⋅ eυ
g µυ = 0 si µ ≠ υ

• Dans l’espace de Minkowski : x µ ⋅ xµ = (c 2 )t 2 − x 2 − y 2 − z 2 (c = 1)

1 0 0 0  xυ = g µυ x µ = ( x 0 ,− x1 ,− x 2 ,− x 3 )
 
 0 − 1 0 0 
g µυ = composantes cov ariantes
0 0 −1 0 
 
 0 0 0 − 1 ds 2 = dt 2 − dx 2 − dy 2 − dz 2 2
 
• Quadrivitesse :

 
( )
1
d ds dt −
v= (t , x ) = ( yc, γv ) , dτ = = , γ = 1− β 2 2
dτ c γ

dt 1
dτ = = dxµ dx µ : temps propre
γ c
µ dx µ 
v = = (cγ , γv ) : forme tensorielle
dτ


• Quadri-impulsion : p µ = mv µ = (γmc, γmv )

 
p = ( E , p ) ⇒ p 2 = g µυ p µ pυ = E 2 / c 2 − p 2 = m 2 c 2

mc 2 E 2 = p 2c 2 + m2c 4
E = γmc 2 = 
1− β 2 T = p 2 c 2 + m 2 c 4 − mc 2
T = E − mc 2 = (γ − 1)mc 2

3
• Transformations de Lorentz :

Référentiels R (fixe) et R’ (vitesse ve // ox) : (c = 1)

 x′ = γ e ( x − β e t )  γe − β eγ e 0 0
 
 y′ = y −γ β γe 0 0
 x′υ = Λυµ x µ , Λυµ =  e e
 0 0 1 0
 z′ = z  
ct ′ = γ e (t − β e x)
 0
 0 0 1 

Composition des vitesses :

v′υ = Λυµ v µ
v ′1 γ v′ v ′x v x − ve
 v′ 
0
 1 − βe 0 0  v  0
= v′x =
 1   1 v′0 γ v′ 1 − ve v x
 v′   − βe 1 0 0 v 
= γ v′2 γ v′v′y
 v′ 2  0   v2 
e
0 0 1/ γ e = v y 1 − ve2
     ′ 0
γ v′ v′y =
 v′3   0 0 1 / γ e   v 3 
v
   0 1 − ve v x

4
 Composition des 4-p :  p′0 = γ e ( p 0 − β e p1 )
 1
 p′ = γ e ( p − β e p )
1 0

 2
 p′ = p
2

 p′3 = p 3

On vérifie que : p′µ p′ µ = m 2 (la masse est un invariant relativiste)

Cas d’une transformation avec une direction quelconque :

 
 v −β 
v′// = //  e : composante parallèle à β e
1 − βe ⋅ v
 1 − β e2  
v⊥′ =   v⊥ : composante orthogonale à β e
1 − βe ⋅ v

 E′   γ e − γ e β e  E   
   =     , p⊥′ = p⊥
 p′//   − γ e β e γ e  p// 

Si R’ est le réf. propre de la particule : p = (m,0) (βe = β)
E = mγ
  
p = γmβ = Eβ 5
 Rapidité :
→ utilisée pour représenter la direction de la trajectoire d’une particule

 E ′ = γ ( E − βp3 ) E ′ − p3′ = λ ( E − p3 )
 p′ = p  1
 1 1 v //(oz ) E ′ + p3′ = ( E + p3 )
 λ
 p2′ = p2 1+ β
 p3′ = γ ( p3 − βE ) λ=
1− β

p3
• Particule au repos dans son réf. : (p’1=0,..,E’=m) β=
E
1 1 + β 1 E + p3
rapidité : ζ = ln λ = ln = ln
2 1 − β 2 E − p3
ζ = tanh −1  p3 E  ζ ∈ ]− ∞,+∞[
 


p = ( E , p ) = (m cosh ζ , pT = p12 + p22 = 0, p3 = m sinh ζ )

• Cas général : impulsion p n’est pas // à l’axe (oz) :


p = ( E , p ) = (mT cosh ζ , pT = p12 + p22 , p3 = mT sinh ζ )
mT2 = m 2 + pT2 = m 2 + p12 + p22 : masse transverse 6
Transformation de Lorentz en fonction de ζ :

−1
 cosh ζ e − sinh ζ e 
γ e = (1 − β )
2 2 2
−1
2 2
=   = cosh ζ e
cosh ζ e
e 2
 

 E ′   cosh ζ e − sinh ζ e  E   
  =    , p⊥′ = p⊥
 p3′   − sinh ζ e cosh ζ e  p3 

β3 − βe
β 3′ = , ta nhζ ′ = tanh(ζ − ζ e ) ⇒ ζ ′ = ζ − ζ e
1 − β3βe

→ même loi de composition qu’en mécanique newtonienne

• Pseudorapidité :
Cas ultra-relativiste : m << p , E ≈ p , mT = ET = pT
 
1  E + p3  1  + p cos θ 
 = − ln tan θ  = η
p
ζ − = ln  ≈ ln  
2  E − p3  2  p − p cos θ   2
7
 Référentiel du centre de masse :

Soit la collision dans R de deux particules m1 et m2 :

 
Centre de masse défini par : ∑ pi′ = 0 i

Energie disponible dans le CM :

(
ECM = ( p1 + p2 ) 2 = E1CM + E2CM = m12 + m22 + 2 E1 E2 (1 − β1β 2 cos θ ) ) 1
2  
, θ = ( p1 , p2 )

Si m2 est au repos (cible fixe) :

(
ECM = ( p1 + p2 ) = m + m + 2 E1m2
2 2
1
2
2 )
1
2

Cas d’un collisionneur : cosθ = -1

= (m + 2 E E (1 + β β ) )
1
ECM = E1CM + E2CM 2
1 +m 2
2 1 2 1 2
2
8
• Cible m2 au repos :
   
ptot p1L E + m2
β CM = βe = = ⇒ γ CM = 1L
Etot E1L + m2 ECM
   m
⇒ p1CM = p2CM = p1L 2
ECM

→ Eq. au cas classique : m1 + m2 ↔ ECM

* Cas m1 << m2 (e + p) :

γ 1m1v1 v1
β CM = =
γ 1m1 + m2 1+
m2
γ 1m1
m2 >> γ 1m1 ⇒ β CM négligeable

1
1+ γ1  2
* m1 = m2 : γ CM = 
 2 

9
Réactions avec 2 particules dans la voie de
sortie : variables de Mandelstam
Soit la réaction : 1+2→3+4
p3
p1 p3
p1 p2
p4
p2 p4

Syst. Labo. Syst. C.M.

 
s = ( p1 + p 2 ) 2 = ( p3 + p 4 ) 2 = m12 + 2 E1 E 2 − 2 p1 ⋅ p 2 + m22
 
t = ( p3 − p1 ) 2 = ( p 4 − p 2 ) 2 = m12 − 2 E1 E3 + 2 p1 ⋅ p3 + m32
 
u = ( p3 − p 2 ) 2 = ( p 4 − p1 ) 2 = m12 − 2 E1 E 4 + 2 p1 ⋅ p 4 + m42
10
• s : carré de la somme des énergies (i ou f) dans le CM
  
s=E 2
CM : p1CM + p2CM = 0

• u : carré du transfert d’énergie impulsion

• s, t et u : invariants de Lorentz
4

On montre que (TD) : s + t + u = ∑ mi

2
•
i =1
• Energies et impulsions dans le C.M.:
s + m12 − m22 s + m22 − m12
E1CM = , E2CM =
2 s 2 s
  λ ( s, m12 , m22 )
pCM = p1CM = p2CM =
2 s
s + m32 − m42 s + m42 − m32
E3CM = , E4CM =
2 s 2 s
  λ ( s , m 2
, m 2
)
'
pCM = p3CM = p4CM = 3 4

2 s
λ ( x, y, z ) = x 2 + y 2 + z 2 − 2 xy − 2 xz − 2 yz 11
• Energies et impulsions dans le laboratoire (cible fixe & p1 // (oz)) :

p1 = (E1L ,0,0, p1L ) s = m12 + m22 + 2m2 E1L

p2 = (m2 ,0,0,0 ) 
t = m12 + m32 − 2 E1L E3 L + 2 p1L p3 L cos θ L

p3 = (E3 L , p3 L ) 4

p4 = (E4 L , p4 L ) u = ∑ mi2 − s − t
i =1

12
 Section efficace

Matrice de diffusion S :
1 + 2 → 3 + 4 + 5 + --- + n (i → f)
Système dans les états |a> à ta et |b> tb :
C fi (t a , tb ) = b U (t a , tb ) a , UU + = I
U ↔ V : H = H0 + V
t a → − ∞ et tb → +∞ : a et b sont des états propres de H 0
Opérateur S : S = lim U (t a , tb ) , S fi = f U (−∞,+∞) i = f S i
t a ,t b → ∞

On sépare la matrice S en deux termes : S fi = δ fi + iT fi

S fi = δ fi + iT fi = δ fi + i (2π ) δ 4 (Pf − Pi )M fi
4

= (2π ) δ 4 (Pf − Pi )δ 4 (0)∑ M fi

2 2 2
i ≠ f : P = C fi (−∞,+∞) = f S i
8

f
13
 Taux de transition par unité de volume :

= (2π ) δ 4 (Pf − Pi )∑ M fi
1 dP 2
W fi = 1 V
4
δ 4 (0) = ∫d x = ∫
4
dt
V dt f (2π ) 4
(2π ) 4

• Section efficace :
(2π )
∑ δ 4 ( p f − pi ) M fi ⋅ dN 3dN 4 ⋅ ⋅ ⋅ ⋅dN n
4
W fi 2
dσ = =
φ φ f

d4 pf
( )
.δ p 2f − m 2f θ (E f )
n
dN 3 dN 4 ⋅ ⋅ ⋅ ⋅dN n = ∏
f =3 (2π )
3

 
φ : flux incident % mod(v1 - v2)

φ =2 ( p1. p2 )2 − m12 m22 : in var iant de Lorentz

 n d4p 
σ=
1
2 2 ∫
∏ ⋅ δ p 2
− m 2
 ⋅ (2(
π )4 4
)
δ ( p f − pi ) M fi
2

2 ( p1 − p2 ) − m1 m2  f =3 (2π )
3 f f

2

14
 Remarque : Espace des phases

• Cas non relativiste :

dN : nombre d’états quantiques dans dV = d3x.d3p

d 3x ⋅ d 3 p 1 1
∫ dN = ∫d p = (2π ) ∫
dpdΩ
3 2
dN = p
(2π ) 3
(2π ) 3 3

• Cas relativiste : d3p → d4p

2
⋅ − − m2 )
( ) ( )
3 2
 d p d ( E p 2
∫ = ∫∫ ⋅ δ − − = ∫∫ δ
⋅ − − m2
4 3 2 2 2 2
d p d pdE E p m E p
2E

( )
3
d p 
=∫ δ E 2 − p 2 − m2
2E

d3p n’est pas un invariant de Lorentz

d3p/E : invariant de Lorentz p′x = γ ( p x − βE ) , E ′ = γ (E − β p x )
dp′x dp x d 3 p d 3 p′
= ⇒ =
E′ E E E′
15
• Cas d’une diffusion 1 + 2 → 3 + 4 :

d 4 p3 ⋅ d 4 p4
σ (1 + 2 → 3 + 4) =
1
∫ (2π )6 [δ ( p3 − m3 )θ ( E3 )] ⋅
2 2

(
2 λ s, m12 , m22 )
[δ ( p42 − m42 )θ ( E4 )](2π ) δ 4 ( p3 + p4 − p1 − p2 ) M fi
4 2

 
d 3 p3 d 3 p4 4
σ (1 + 2 → 3 + 4) = ∫ 2 E3 2 E4 δ ( p3 + p4 − p1 − p2 ) M fi
1 2

2(2π )
2
λ (s, m12 , m22 )
On obtient finalement :

p3CM
σ (1 + 2 → 3 + 4) = ∫ d (cosθCM )
1 2
 M fi
32π s p1CM s

dσ 1 p3CM 2
=  M
d cos θ CM 32πs p1CM
fi

En utilisant l’expression de t :

dσ 1 1 2
= M fi
dt 64πs p1CM 2

16
• Domaine de variation de t :

 
t = ( p1 − p3 ) 2 = (E1CM − E3CM ) − ( p1CM − p3CM )
2 2

θ
= (E1CM − E3CM ) − ( p1CM − p3CM )
   
− 4 p1CM p3CM sin 2 CM
2 2

) = (m ) − ( p
2
− m22 − m32 + m42
− ( p1CM − p3CM )
2
  
t (θ CM = 0 ) = (E1CM − E3CM ) − p3CM
2 2 1 2
1CM
4s
 
t (θ CM = π ) = t (θ CM = 0 ) − 4 p1CM p3CM

• Canaux croisés :
A partir du canal a + b → c + d : prédictions pour les canaux croisés
a +b→c +d (canal s )
a +d →c +b (canal u )
d +b→c +a (canal t )
Eléments de matrice :
       
M ab→cd ( pa , pb , pc , pd ) = M ad →cb ( pa ,− pd , pc ,− pb )
   
= M d b→ca (− pd , pb , pc ,− pa )
M ab→cd ( s, t , u ) = M ab →cb (u , t , s ) = M d b→ca (t , s, u ) 17
 Désintégrations de particules

• Situation analogue au cas précédent avec 1 seule particule dans la voie

d’entrée : 1 → 2 + 3 + 4 + …..+ n

• Taux de désintégration donné par la largeur Γ : Γ.τ ≥ 1 (Γ = ΔΕ)

τ : vie moyenne de la particule

• Taux de désintégration :
n d4p
Γ(1 → 2 + 3 + ⋅ ⋅ ⋅ + n ) = 2m1ϖ fi =
1
∫ ∏ f
⋅ δ 4
( p − p ) M
2
= τ −1
2m1 f = 2 (2π )3 f i fi

Γ = ∑ Γi , Γi : largeur partielle de la voie de désintégration i

i

18
 Particule au repos :
t
− 0

probabilité que la particule vive un temps supérieur à t0 : P(t 0 ) = e τ

 Particule en mouvement (E , p) :
m1 t 0 m1 x0 m x
− − − 1 0
E1 τ
P(t0 ) = e = P( x0 ) ; x0 = β1t0 P( x0 ) = e E1 β1τ
=e
p1 τ

τ : mesuré dans le référentiel lié à la particule

Désintégrations en deux corps :

Dans le réf. lié à la particule (C.M.) :


m12 + m22 − m32 1 p 2CM
=
2
E2CM dΓ = M fi dΩ 2
2m1 32π 2
m1
m12 − m22 + m32 dΩ 2 = dϕ 2 d (cos θ 2 )
E3CM =
2m1

 
p2CM = p3CM =
{(
m − (m2 + m3 ) m − (m2 − m3 )
2
1
2
)( 2
1
2
)}
1
2

2m1
19
 Désintégration en trois corps :

Espace des phases pour trois particules :


 4 1 d 3 pi 
dN = dN 2 dN 3 dN 4 = ∏  ⋅ (2π ) δ (Pi − Pf )
4 4

 i = 2 (2π ) 2 Ei 
3

 2  2
1 p2 p3 dΩ 2 dΩ 3
∫ dN = ∫ dp 2 dp3 δ (E1 − E2 − E3 − E4 )
8 E2 E3 E4
 2  2
p2 p3 d (cos θ 23 )
=π ∫2
dp2 dp3 δ (E1 − E2 − E3 − E4 )
E2 E3 E4

  dN
si p2 et p3 sont fixes : = cste
dE2 dE3

→ en absence des effets de Mfi : distribution plate de E2 et E3

20
• Cas général :
on définit les invariants de Lorentz :
s 24 = ( p 2 + p 4 ) 2 , s34 = ( p3 + p 4 ) 2

∫ ds δ (s
24 24 ) (
− ( p 2 + p 4 ) 2 = ∫ ds34δ s34 − ( p3 + p 4 ) 2 = 1 )
 
d 3 p2 d 3 p3
( ) (
∫ dN = ∫ (2π )5 8E2 E3 E4 ds24 ds34 .δ (E1 − E2 − E3 − E4 ).δ s24 − ( p2 + p4 ) .δ s34 − ( p3 + p4 )
2 2
)

Dans le C.M . :
    
E1 = s = m1 , p1 = 0 , p4 = −( p2 + p3 )
 2   2 
d 3 p2 = 2π p2 d p2 d (cos θ 23 ) , d 3 p3 = 4π p3 d p3

2 2
1 M fi 1 M fi
dΓ = dE 2 dE3 = ds 24 ds34
(2π ) 3
8m1 (2π ) 3
32m1

21
• Si Mfi ne dépend pas explicitement de s24 et s34 : distribution uniforme
2
2 M0
• Si Mfi dépend explicitement de s24 et s34 : M fi ∝
s 24 − m 2 + iΓm

Diagramme de Dalitz

( )
2
(s34 )min,max = E + E −  E3* 2 − m32  E 4* 2 − m42 
* * 2

 
3 4

m12 − s 24 − m32 * s 24 − m22 + m42

E3 =
*
, E4 =
2 s 24 2 s 24
22
 Exercices :
4
1- Montrer que : s + t + u = ∑ mi2
i =1

2- Calculer l’énergie seuil des réactions : ν e + p → n + e + , p + p → p + p + p + p

le proton cible est au repos. Que devient l’énergie seuil de la deuxième
réaction lorsque le proton cible est en mouvement.

3- Un méson π0, de vitesse v, se désintègre en vol, en deux photons.

Déterminer la valeur des angles de diffusion des deux photons par
rapport à la direction de la particule incidente.

23