Interaction des particules

avec la matière
Table des matières
1 Introduction : ..............................................................................................................................................2
2 Interaction des particules chargées avec la matière :.................................................................................4
2.1 Énergie cinétique maximale transférée à un électron d'ionisation :...................................................4
2.2 Perte d'énergie de particules lourdes par ionisation et excitation dans la matière :...........................6
2.3 Perte d'énergie des électrons et des positons par ionisation et excitation dans un milieu :...............8
2.4 Perte d'énergie dans un milieu composé de plusieurs éléments :......................................................9
2.5 Rayonnement de freinage : Bremsstrahlung :....................................................................................9
2.6 Production directe de paires électron-positon :................................................................................11
2.7 Perte d'énergie par interactions photo-nucléaires :...........................................................................12
2.8 Perte d'énergie totale :.......................................................................................................................12
2.9 Fluctuations angulaires dues aux diffusions multiples coulombiennes :.........................................12
2.10 Parcours des particules chargées :..................................................................................................13
3 Énergie déposée :......................................................................................................................................16
3.1 Fluctuation du dépôt d'énergie :........................................................................................................16
3.2 Perte d'énergie par ionisation et excitation tronquée :......................................................................16
3.3 Électrons secondaires :......................................................................................................................17
3.4 Création de paires d'électrons et d'ions par ionisation :...................................................................17
4 Interactions des photons dans un milieu :................................................................................................18
4.1 Effet photo-électrique :.....................................................................................................................19
4.2 Effet Compton : ................................................................................................................................20
4.3 Production de paires e+e- :...............................................................................................................21
4.4 Section efficace totale d'absorption :................................................................................................22
5 Interaction forte des hadrons :..................................................................................................................23
6 Pour en savoir plus :.................................................................................................................................25

1
1 Introduction :
Eu égard à la diversité des interactions et des caractéristiques des particules, on pourrait penser à
priori que la question traitée dans ce chapitre est complexe et requiert un exposé long . En réalité, un
examen rapide des libres parcours moyens dans le matière caractérisant chaque type d'interaction, va
permettre de simplifier grandement la suite.
Le libre parcours moyen exprime la distance moyenne parcourue par une particule entre deux
« chocs » dus à un processus donné. Il correspond à :
1
=
n
où :
•  est la section efficace du ( ou des ) processus considérés ;
• n est le nombre de centres diffuseurs ( atomes, noyaux...) par unité de volume de
NA
matière . Exemple : n= représente le nombre d'atomes ou de noyaux par unité de
M
volume pour un corps monoatomique ; NA étant le nombre d'Avogadro et M la masse
molaire .
On peut donner quelques valeurs typiques de  pour les interactions connues :
• Électromagnétique : 1  m ( particules chargées )
• Forte : 1 cm ( neutrons , ....)
• Faible:  1015 m ≃ 0,1 année lumière (neutrinos)
Dun point de vue pratique , un signal exploitable pour la détection de particules ( nombre de chocs
100 } ne peut provenir que de l'interaction électromagnétique , donc de particules chargées . Cela ne
signifie pas pour autant que les particules neutres sont indétectables . Leur détection doit procéder par
étapes : 1) réaction primaire de conversion en particules chargées , 2) interaction des particules chargées
secondaires donnant lieu aux signaux dans les détecteurs . Voici quelques exemples choisis qui illustrent
bien ce principe :
détection de photons :

2
1 Diffusion Compton
e-

3 e-
Création d'une paire
e+ électron-positon

2
 détection de pions neutres :
Un 0 se désintègre en 2 photons avec une vie moyenne de 8,4 10-17 s .

e+
1
0 e-
e-
2
e+

détection de neutrinos :
Un réacteur nucléaire développant une puissance thermique de 2800 MW , émet environ 130 MW
sous la forme de neutrinos . Un détecteur d'environ 1 m 3 placé à environ 20 m du coeur peut espérer
enregistrer jusqu'à 100 neutrinos par heure .

e
e Réacteur nucléaire Liquide scintillant ( CHx + 6Li )
e

e e
e p 
Li6
t
n
e
e
~0

La réaction de détection est : 
 e p ne , qui est rapidement ( 100 µs ) suivie de la capture du
6 6
neutron sur un noyau de Li selon la réaction : nth Li  t 4,8 MeV . Les particules chargées
produisent des impulsions de scintillation en coïncidence . La signature de détection d'un neutrino
correspond à l'enregistrement de deux impulsions lumineuses induites par le positon et la paire −t .

3
2 Interaction des particules chargées avec la matière :
L'ionisation et l'excitation coulombienne d'atomes sont les processus dominants dans l'interaction de
particules chargées avec la matière . Un particule P d'état de charge Z peut provoquer les réactions
suivantes :
• PZ atome atome ∗P Z  , excitation suivie par atome ∗ atome : désexcitation
 Z −  Z
• P atome atome e P : ionisation
 Z ∗ −  Z
• P atome atome  e P : ionisation et excitation

2.1 Énergie cinétique maximale transférée à un électron d'ionisation :

Nous allons calculer l'énergie cinétique maximale qui peut être transférée par ionisation entre un
particule chargée incidente et un électron du cortège électronique :

m0 v me

p =m 0 v
atome
E Z

On fait l'hypothèse que la vitesse de la particule incidente est très supérieure à la vitesse orbitale des
électrons , c'est-à-dire : v ≫ 〈 v e 〉≃ Z  c pour l'orbite la plus profonde , où  est la constante de
structure fine ( =1/137 ) et c est la vitesse de la lumière dans le vide. L'électron sera donc considéré
au repos dans le choc .
Les règles de la cinématique relativiste permettent de trouver les quantités suivantes :
2 2 1/ 2
• ECM =  m 0me 2 m e E : énergie disponible dans le centre de masse de la collision , c'est-à-
dire là où ∑
p=0 ( en unités naturelles ℏ=c=1 ) ;
CM me
• pe = p : quantité de mouvement de l'électron dans le centre de masse ;
E CM
me
• ECM
e =Eme : énergie totale de l'électron dans le le centre de masse ;
E CM

CM Eme CM p
•  = ,  = : facteurs de la transformation de Lorentz du repère du laboratoire
E CM Eme
vers le centre de masse .
Sachant que la collision est élastique et que l'énergie maximale transférée à l'électron éjecté correspond à
4
la situation suivante dans le CM :

p CM
0 p CM
e

p CM

on peut alors écrire :

   
CM
E emax =  CM CM ECM
e
pe CM CM CM pCM
e

soit :
max  Eme  p2 m e
E e = 2
m e Em e  2
E CM E CM
me
T max max
e =Ee −m e= 2
 E2 −m02 p2 
E CM
2
2me p 2 me 2 2
= 2
=
E CM  ECM /m02
si m0 >> me et si l'énergie incidente est faible :
2 2
2 m0 me 2 me E
On a alors : E CM / m 0 =  2
 2
 ≃1 , avec E= m0
m0 m0 m 20
2  me
C'est-à-dire si la particule chargée n'est pas un électron et si la condition ≪1 est satisfaite
m0
( pour un proton si E p 50 GeV , pour un muon si E  500 MeV )
On obtient alors :
T max 2 2
e =2 me  

si m0 = me :
En d'autres termes, la particule incidente est un électron, l'expression du carré de l'énergie du centre
de masse se réduit à :
E CM  2=2 m2e 2 me E=2 me me E  ,
d'où :
max E 2−me 2
T = = E−me  .
e
meE
On voit apparaître deux situations : soit la particule incidente est beaucoup plus massive qu'un électron
2 m
et, dans l'hypothèse où son énergie incidente est faible ( m ≪1 ), on peut appliquer la première
0
e

formule, soit celle-ci est un électron et de ce fait l'énergie transférée maximale est beaucoup plus grande.
En pratique ces deux situations correspondent bien aux deux classes de particules chargées rencontrées
: électrons ( me) ou pions, muons, protons ....( mo >> me ).

5
2.2 Perte d'énergie de particules lourdes par ionisation et excitation dans la
matière :
Ici nous traitons le cas des particules dont les masses sont très supérieures à celle d'un électron . Pour
ces particules , la perte d'énergie moyenne ( par unité de longueur ) , que l'on appelle également le
pouvoir d'arrêt ( stopping power ou average energy loss en anglais ) , est donné par la formule de Bethe
et Bloch ( voir référence 1 pour une démonstration ) :
dE 2 2 Z 1 1 2 me 2 2 T max
e 2  Ce
− =4  N A r e me z [ ln  − − − ] ( avec ℏ c=1 ) ,
dx A 2 2 I2 2 Z
dans laquelle :
• me est la masse de l'électron ;
• r e est le rayon classique de l'électron ( r e= / me ) ;
•  est la constante de structure fine ( =1 / 137 ) ;
• N A est le nombre d'Avogadro ;
• z est l'état de charge de la particule incidente ;
• Z est le nombre atomique du milieu dans lequel la particule se propage ;
• A est la masse atomique de ce milieu ;
•  ,  sont les facteurs de Lorentz de la particule incidente ;
• T max
e est l'énergie cinétique maximale qui peut être transférée à un électron d'ionisation;
• I est l'énergie moyenne d'excitation qui est une caractéristique du milieu. I peut être
approximée par : I = 16 Z0,9 eV pour Z > 1, I = 15 eV pour de l'hydrogène gazeux, I =
19,2 eV pour de l'hélium gazeux et I = 21,8 eV pour de l'hydrogène liquide ;
•  est un paramètre qui décrit l'écrantage du champ électrique de la particule incidente
dans le milieu à haute énergie . C'est un paramètre qui dépend de la densité du milieu ;
• Ce / z est un terme de correction d'effet de couches atomiques qui tient compte du fait
qu'à basse énergie les particules incidentes ont peu de chance d'interagir avec les électrons
installés sur les orbites les plus profondes du cortège atomique ;
• dx est la densité surfacique du milieu , c'est-à-dire dx= dl , où  est la masse
volumique du milieu et dl est l'élément infinitésimal de distance parcourue dans le
milieu.
dE
Dans le système d'unités naturelles ℏ c=1 , −a la dimension de l'inverse d'un carré d'une
dx
énergie . Pour revenir à un système d'unités plus conventionnelles, on utilise la formule suivante :
dE MeV dE 2 2
− [ ]=−  c  ℏ c ,
dx g/cm 2
dx  ℏ c=1
−13
avec : c 2 = 5,61 1026 MeV/g et ℏ c=197 10 [Link] ,
2 2 5 MeV 3
c'est-à-dire : ℏ c  c =2,1810 .
g /cm2

6
Une forme plus pratique de la formule de Bethe et Bloch est la suivante :
dE MeV
2
0,3071 z Z 1 2 m e 2 2 T emax  C
− [ 2
] = 2
 ln  2
− 2 − − e  .
dx g / cm A g  2 I 2 Z
La figure qui suit présente le pouvoir d'arrêt pour plusieurs particules en fonction de leur quantité de
mouvement :

Remarques :
dE
• pour  1 , − ∝ −5 / 3 ;
dx
dE
• pour  ≃ 4 , − atteint un minimum assez large . Une particule qui a une quantité de
dx
mouvement (  = p / m ) correspondant à cette région est appelée une particule au minimum
d'ionisation ( Minimum Ionizing Particule en anglais ou MIP ) . Pour les éléments légers , la
perte d'énergie des particules au minimum d'ionisation est approximativement :

7
 dE MIP MeV
− ≃2 .
dx g / cm2
• pour  4 , on entre dans la région qui correspond à la remontée relativiste ( ∝2 ln  ).
Dans ce régime, le champ électrique de la particule incidente a une portée transverse de plus en
plus grande. Dans le milieu, cette portée est réduite par des effets de polarisation des électrons
atomiques. Ce mécanisme est pris en compte par le terme − / 2 ( = ln /2 où 
dépend du milieu ; pour des particules très énergétiques ).
2 me max 2 2
• à basse énergie, c'est-à-dire si ≪ 1 , on a T e =2 me   et par conséquent :
m0
2 2
dE MeV
2
0,3071 z Z 2 me    C
− [ 2
]= ⋅ 2
[ln  − 2 − − e ]
dx g / cm A g  I 2 Z
• la formule de Bethe et Bloch donnée dans ce cours est précise à quelques pour-cent . Aux
faibles vitesses ( v  Z c ) , il faut prendre en compte d'autres corrections. Pour plus de
précision , il faut consulter la référence 2.
• pour des ions très lourds dans des milieux à Z élevé, la condition v  Z c est très rapidement
insatisfaite lorsque l'on va vers les basses énergies. Dans ce régime, les ions capturent des
électrons et réduisent ainsi leur état de charge. Le terme z 2 dans la formule de Bethe et Bloch
2 2
doit alors être remplacé par une charge effective qui tient compte ce cet effet ( z Q  , I  )
(voir référence 3).

2.3 Perte d'énergie des électrons et des positons par ionisation et excitation
dans un milieu :
Pour les cas où les particules incidentes sont des électrons, le calcul doit tenir compte du fait que les
particules qui participent aux collisions sont de nature identique : transfert d'énergie plus important,
particules indiscernables dans la voie de sortie . La perte d'énergie s'exprime par :
dE MeV 0,3071 Z 1 T me  2 2 1 1 −1
2
− [ ]= ⋅ [ ln   2  1−2 −1 ln 2   ]
dx g /cm2 A g  2 2 2I2 2 16 
où T est l'énergie cinétique de l'électron incident : T = −1me =E−me .

Lorsque les particules incidentes sont des positons , on peut distinguer les produits dans la voie de
sortie. La perte d'énergie est alors donnée par :
2 2 2
dE MeV 0,3071 Z 1 T me    14 10 4
− [ 2
]= ⋅ 2 [ ln 2
−  23  2
 3
]
dx g / cm A g  2 2I 24 1  1 1
Cependant, si les positons perdent la totalité de leur énergie dans le milieu, ils s'annihilent ensuite
suivant :
e + e + 1,022 MeV

8
Les positons peuvent également s'annihiler en vol avec la section efficace suivante :
z  r e2 24 1 3
ln    −1 −
2
 Z , E= [ ]

2
1  −1  2
−1

2.4 Perte d'énergie dans un milieu composé de plusieurs éléments :

Pour une particule se propageant dans un milieu composé, la perte d'énergie est donnée par
l'expression suivante :
dE dE
≈∑ f i ∣ ,
dx i dx i
mi dE
où fi est le rapport massique de l'élément i ( , ∑ m i=m ) ;
f i= ∣ est la perte d'énergie
m i dx i
moyenne dans l'élément i . Cette formule ne tient pas compte d'effets moléculaires qui peuvent se
produire à basse énergie.

2.5 Rayonnement de freinage : Bremsstrahlung :

Toute charge décélérée (ou accélérée) émet un rayonnement électromagnétique sous la forme de
photons. La perte d'énergie occasionnée par ce phénomène pour une particule chargée à haute énergie
dans le champ électrique des noyaux (le plus intense, donc le plus apte à freiner) est donnée par la
formule suivante :
dE rad Z2 2 2 183
− =4  N A z r E ln 1/3  ,
dx A Z
dans laquelle :
NA est le nombre d'Avogadro ;
Z et A sont respectivement le numéro et la masse atomique du milieu ;
 est la constante de structure fine ( 1/137 ) ;
r= / m est le rayon classique de la particule ;
z , m, E sont respectivement l'état de charge , la masse et l'énergie cinétique de la
particule incidente .
Cette expression peut être réécrite sous une forme plus exploitable :
2
dE rad MeV 0,3071  2 2 me E 183
−  = Z z   ln  1/ 3  .
dx g /cm2
A g  m m e Z
Pour des protons traversant du fer, on a :
dE rad MeV E 1 MeV
−  =1,06⋅10−8 soit pour des protons de 50 TeV .
dx g /cm 2
me g /cm2
Pour des muons , toujours dans le fer , on a :
dE rad MeV E 1 MeV
−  =8.43⋅10−7 soit pour des muons de 600 GeV .
dx g /cm2
me g /cm 2
On voit sur ces quelques valeurs numériques que ce processus va essentiellement affecter les particules
les plus légères qui sont effectivement plus sensibles au ralentissement. En pratique, ce sont les électrons

9
qui y seront les plus sensibles. La perte d'énergie par rayonnement de freinage d'une particule d'état de
charge z et de masse m peut être calculée à partir de celle d'un électron possédant la même énergie
cinétique incidente :
2
dE rad m dE rad − .
 z , m= e  z 2 e 
dx m dx
Dans la suite de ce paragraphe, nous nous restreindrons au cas des électrons. Pour obtenir une formule
plus exacte, il faut ajouter la contribution provenant du Bremsstrahlung sur les électrons du cortège
atomique. Ici Z 2=1 , mais il y a Z centres qui contribuent par atome et par conséquent l'effet total est
proportionnel à Z . On obtient :
dE rad − Z  Z1 2 183
− e =4  N r e E ln  1 /3  , où re est le rayon classique de l'électron .
dx A Z
Cette équation peut être réduite sous la forme :
dE rad − E
− e = , où X0 est la longueur de radiation du milieu .
dx X0
L'énergie moyenne rayonnée dans un milieu par un électron ayant parcouru x est donnée par :
−x / X0
Erad e−=E 1−e  , dans laquelle x et X0 peuvent être exprimés en cm ou en g/cm2.
De l'expression de la perte d'énergie par rayonnement, on tire que :
A g 
X 0 g /cm 2=
183
4  N Z  Z1r 2e ln  
Z 1/3
Une formule plus rigoureuse pour le calcul de X0 est obtenue en prenant en compte l'effet d'écrantage du
champ électrique nucléaire introduit par les électrons atomiques :
716,4 A g
X 0 g /cm 2=
287
Z  Z1 ln 1 /2 
Z
Dans la littérature, les longueurs de radiation sont toujours calculées pour des électrons.
Pour un milieu composé de N éléments, la longueur de radiation peut être approchée par l'expression
suivante :
−1
fi
X o =[ ∑ ] , où fi et X0I sont respectivement la fraction de masse et la longueur de radiation de
i X i0
l'élément i .
À haute énergie, la perte d'énergie par rayonnement croît comme l'énergie, alors que la perte d'énergie
par ionisation et excitation est proportionnelle au logarithme de E. On appelle énergie critique, la valeur
de l'énergie pour laquelle on a égalité entre les deux pertes d'énergie :
dE rad dE ionisation
 Ec =  Ec 
dx dx

10
Une valeur approximative (à 10% près) de Ec peut être calculée à l'aide des formules suivantes :
610 MeV
Ec = pour les liquides et les solides
Z1,24
710 MeV
Ec = pour les gaz
Z0,92

milieu Z A X0 (g/cm2) X0 (cm) EC (MeV)

hydrogène 1 1.01 63 700000 350
hélium 2 4 94 530000 250
lithium 3 6.94 83 156 180
carbone 6 12.01 43 18.8 90
azote 7 14.01 38 30500 85
oxygène 8 16 34 24000 75
aluminium 13 26.98 24 8.9 40
silicium 14 28.09 22 9.4 39
fer 26 55.85 13.9 1.76 20.7
cuivre 29 63.55 12.9 1.43 18.8
argent 47 109.9 9.3 0.89 11.9
tungstène 74 183.9 6.8 0.35 8
plomb 82 207.2 6.4 0.56 7.4
air 7.3 14.4 37 30000 84
silice ( SiO2) 11.2 21.7 27 12 57
eau 7.5 14.2 36 36 83

2.6 Production directe de paires électron-positon :

À très haute énergie, dans le champ électrique des noyaux, la production directe – par des particules
chargées - de paires d'électrons et de positons peut devenir un processus quantitativement important.
La perte d'énergie moyenne par production de paires e+e- est également proportionnelle à l'énergie de la
particule incidente :
dE paires
− =b paires  Z , A , E E
dx

Le paramètre bpaires(Z,A,E) ne varie que très lentement avec E. Concernant les muons, ce processus
conduit à une perte d'énergie moyenne à haute énergie qui est plus importante que celle qui est due au
Bremsstrahlung.
Exemple : muons traversant du fer
dE paires Fe MeV
− ∣ =3⋅10−6 E  
dx g/cm 2
11
2.7 Perte d'énergie par interactions photo-nucléaires :
Par échanges de photons virtuels avec les noyaux, les particules chargées peuvent produire des
réactions inélastiques.
Exemple : électro-dissociation du deuton : e− d n pe−
La perte d'énergie induite par ces processus suit une loi similaire à celle qui traduit la production de
paires e+e- :
dE  nucl.
− =b  nucl. Z , A , E E
dx
dE  nucl. Fe MeV
Pour des muons de 100 GeV traversant du fer , on a : − ∣ 100 GeV =0,04
dx g/cm2

2.8 Perte d'énergie totale :

La perte d'énergie totale est donnée par :
dE tot dE ionisation dE rad dE paires dE nucl.
=   
dx dx dx dx dx
Contrairement aux pertes par ionisation, celles qui mettent en jeu le rayonnement de freinage, la
production de paires e+e- ou les réactions photo-nucléaires se caractérisent par de grands transferts
d'énergie et de larges fluctuations. La perte totale d'énergie peut également être exprimée sous la forme:
dE tot
− =a Z , A , E  b Z , A , E E
dx
où : a(Z,A,E) est le terme correspondant aux pertes par ionisation et excitation ;
b(Z,A,E) correspond à la somme des contributions dues au Bremsstrahlung , à la production de
paires e+e- et aux réactions photo-nucléaires.

Pour les muons, les événements à larges transferts d'énergie provenant du rayonnement de freinage, de la
production de paires e+e- ou des réactions photo-nucléaires sont parfois appelées interactions
catastrophiques.

2.9 Fluctuations angulaires dues aux diffusions multiples coulombiennes :

En traversant un milieu, une particule chargée est sujette à de nombreuses interactions dans le champ
coulombien des noyaux et des électrons. À chaque choc, cette particule est légèrement défléchie par
proj
rapport à sa direction incidente ( principalement sur les noyaux ). La distribution de l'angle  due à
la diffusion coulombienne est décrite par la théorie de Molière. C'est une distribution centrée sur 0 qui a
proj
une partie centrale (98% du total) qui peut être approchée par une Gaussienne.  est la projection
espace
de l'angle  sur un plan. Pour de petits angles, on a :
espace 2= proj 2 proj 2 proj 2
x   y  =2 
1 espace
 proj = 
2

12
La distribution de l'angle projeté est donnée par :


proj 2
1 
−  
1 2 0 13,6 MeV x x
P proj d proj= e d proj avec 0= z 10,038 ln  ;
 2  0
p X 0 X 0

p étant la quantité de mouvement de la particule incidente ;

X0 la longueur de radiation du milieu ;
x l'épaisseur du milieu ;
z l'état de charge de la particule incidente .
La distribution de l'angle d'espace est obtenue par :
espace 2
1 
−  
1 espace 2
 = proj 2 proj 2 proj proj
d  avec  x   y  et d =d  x d  y
2 0
Pespace  d = e
2 20
Les grandes valeurs de  sont plus probables que ce que prévoient ces distributions gaussiennes.
On se rapproche d'une diffusion de Rutherford due aux noyaux.

2.10 Parcours des particules chargées :

Dans un milieu suffisamment épais, une particule chargée va perdre toute son énergie cinétique et va
−5 /3
s'arrêter. Puisque la perte d'énergie par ionisation a une dépendance en  à basse énergie, le profil
de perte d'énergie en fonction de la distance parcourue aura cette forme caractéristique.

D'une manière générale, le parcours d'une particule chargée peut être défini par l'expression suivante (
approximation de ralentissement continu) :
To
dT
RT 0=∫
dE ,
0
T 
dx

13
où T0 est l'énergie cinétique de la particule incidente (en pratique on stoppe l'intégration à 10 eV ).
Cependant, R est sujet à de larges fluctuations causées par les interactions ou les chocs à grands
transferts d'énergie ( rayonnement de freinage ...) . On peut avoir recours à une autre définition de R, que
l'on appelle le parcours moyen : 〈 R〉 . Le parcours moyen est défini comme étant la distance - évaluée
en ligne droite entre les points de départ et d'arrivée - au bout de laquelle le nombre initial de particules
est réduit de moitié. Pour des énergies incidentes supérieures au MeV, on obtient par cette définition des
valeurs comparables aux parcours calculés selon l'approximation de ralentissement continu : 〈 R〉≈R
avec une très bonne précision.
Si pour le calcul de R(T0) on ne prend que la perte d'énergie par ionisation et excitation (approximation
valable pour des particules lourdes avec T0 inférieure à quelques GeV), on obtient pour un milieu donné
(c'est-à-dire Z, A, I fixés) :
dE ionisation 2
− ≈ z f =z 2 gT / M  ,
dx
où M et z sont la masse et l'état de charge de la particule incidente.
T0 T0 / M
M d T / M
RT 0= 2
z
∫ gT / M = Mz2 ∫ dx M
= hT o / M  ,
g x  z2
0 0

pour laquelle h est une fonction « universelle » du milieu (pour Z, A, I donnés). Ainsi, si on connaît le
parcours Ra d'une particule de masse Ma , de charge za en fonction de son énergie cinétique, le parcours
Rb d'une particule de masse Mb , de charge zb et d'énergie cinétique Tb sera :
2
M b za Ma
R b M b , z b , T b= 2
R a M a , z a ,T b  .
Ma z b Mb

On peut apporter à cette formule une légère correction C   b qui dépend de la vitesse de la particule b
et qui est en général proche de l'unité.
2
M b za Ma
R b M b , z b , T b= 2
R a M a , z a ,T b  C  b .
Ma z b Mb
R
Pour des particules possédant un état de charge de +/- 1, la fonction h est donnée par : h= . Cette
M
fonction est présentée pour quelques matériaux sur le graphique de la page suivante.

exemple : muons portant une quantité de mouvement de 1 GeV dans du fer ,

R/M = 6000 g cm-2 GeV-1 , m  = 0,105 GeV , R = 6 104 g cm-2 = 76 m .

14
15
3 Énergie déposée :
Il est crucial de noter la différence qui existe entre la perte d'énergie d'une particule et l'énergie qu'elle
dépose dans un milieu, par exemple dans une couche active d'un détecteur. Pour des particules rapides,
une fraction importante de l'énergie cinétique incidente est transférée à des particules secondaires
énergétiques qui peuvent ensuite sortir du milieu considéré sans avoir déposé la totalité de leurs énergies.
Du fait de la complexité de ce phénomène, qui met en jeu des particules secondaires qui sont souvent
de nature différente du projectile incident, de son interdépendance avec la géométrie et les
caractéristiques des milieux (des détecteurs), il n'existe pas de formule précise qui puisse être
simplement utilisée pour obtenir l'énergie déposée.
Pour traiter ce problème, on a maintenant recours à des calculs sur ordinateurs de type Monté-Carlo,
qui exécutent une simulation complète de l'histoire d'une particule dans un milieu : parcours, collisions,
génération de particules secondaires, énergie déposée. Le programme de ce type le plus connu est le code
de simulation GEANT ( version 4 en C++ ou 3 en FORTRAN ).

3.1 Fluctuation du dépôt d'énergie :

Celle-ci résulte de la nature stochastique des collisions des particules primaires ainsi que des produits
secondaires dans un milieu matériel. Jusqu'ici, les formules qui ont été indiquées pour les pertes
d'énergie ne donnent que les valeurs moyennes de ces grandeurs. La perte d'énergie par des processus
concurrentiels (ionisation, rayonnement de freinage, ...) est un phénomène très fluctuant, avec des
amplitudes de fluctuations d'autant plus grandes que les transferts d'énergies considérés sont importants.
Pour des épaisseurs traversées fines ou intermédiaires (nombre de collisions réduit), on peut considérer
que la distribution de perte d'énergie est asymétrique et présente des queues à haute énergie assez
prononcées. Une telle distribution est souvent désignée sous le nom de distribution de Landau, du fait
que Landau fut l'un des premiers physiciens à proposer un modèle analytique pour ce phénomène. Même
si la théorie de Landau a été par la suite améliorée, aucune des formules proposées ne reproduit les
distributions expérimentales avec suffisamment de précision pour être utilisée en pratique. Il est
préférable d'avoir recours à une simulation de type de Monté-Carlo.
Pour des épaisseurs traversées fortes, ou lorsque le nombre de collisions devient très grand, la
distribution de perte d'énergie tend vers une Gaussienne.
L'origine des queues à grande énergie des distributions de pertes d'énergie, réside dans les chocs
intermittents à grands transferts d'énergie.

3.2 Perte d'énergie par ionisation et excitation tronquée :

Dans les détecteurs qui utilisent des couches actives très fines, on observe fréquemment que le dépôt
d'énergie fluctue beaucoup moins que prévu. Les queues à grande énergie apparaissent moins peuplées.
Cela provient du fait que les électrons secondaires les plus énergétiques peuvent quitter les couches
actives en ne déposant qu'une fraction assez faible de leur énergie.
Pour tenir compte de cet effet dans le processus d'ionisation et d'excitation, on emploie la forme
tronquée de la formule de Bethe et Bloch en remplaçant T emax par T supe avec T supe  T emax . On
obtient :
dE ionisation MeV 0,3071 z 2 Z 1 2 me 2 2 T supe 2 T  C
− [ 2
]= 2
 ln 2
−  1 supe
e
− − e 
dx supe g/cm A(g)  2 I 2 T max 2 Z

16
 
Du fait de la présence du terme qui prend en compte les effets de polarisation à haute énergie, la
2
perte d'énergie tronquée tend vers un plateau (plateau de Fermi) à   élevé.

3.3 Électrons secondaires :

Les électrons créés par ionisation avec des énergies suffisamment grandes vont à leur tour pouvoir
ioniser le milieu. Ils se comportent alors comme des particules secondaires. On les appelle électrons
(parfois rayons) delta ou électrons secondaires.
Leur distribution par unité d'énergie cinétique et par unité de longueur parcourue par la particule
primaire, est :
d 2 N 0,3071 z2 Z F T 
= ,
dT dx 2 A g  T  MeV
2 2

où :
• T est l'énergie cinétique de l'électron δ produit ;
• z, β sont l'état de charge et la vitesse du projectile ;
• A est la masse atomique du milieu ;
• Z est le nombre atomique du milieu ;
• F(T) est un facteur qui dépend du spin de la particule incidente :
• F T ≃1 si T ≪T emax ;
T
• F T =1−  pour des particules incidentes de spin = 0 .
T emax
Cette formule est dérivée dans l'hypothèse où T ≫ I (énergie moyenne d'excitation). Lorsque les
e
T
particules incidentes sont des électrons, T ne peut atteindre que 2
du fait que les particules dans la
max

voie de sortie sont de natures identiques.

3.4 Création de paires d'électrons et d'ions par ionisation :

Lorsque le signal mesuré par un détecteur est un courant ou une charge électrique créés par
ionisation, il est utile de calculer la moyenne du nombre total de paires d'électrons et d'ions libérés par le
dépôt d'énergie :
−
−ion  Edéposée
ne = où W est l'énergie moyenne nécessaire pour créer une paire électron-ion
W
dans le milieu considéré.
En général, W I (énergie moyenne d'excitation) car une fraction de l'énergie déposée est absorbée
par des excitations qui ne conduisent pas à la libération d'électrons.
Pour les gaz courants, W est voisin de 30 eV. Dans les détecteurs utilisant des semi-conducteurs en
tant que milieu actif, W est beaucoup plus faible d'où leur intérêt : W=3,6 eV pour Si et 2,85 eV pour Ge.

17
4 Interactions des photons dans un milieu :
Les photons sont détectés indirectement par création de particules chargées (majoritairement e -, e+)
dans le milieu où ils se propagent. La nature de leurs interactions diffère fondamentalement du processus
d'ionisation décrit jusqu'ici pour les particules chargées, car dans leurs interactions, les photons sont soit
totalement absorbés (effet photo-électrique, création de paires e+e-), soit diffusés (effet Compton) sous
des angles qui peuvent être assez larges.
Du fait de leur trajectoires chaotiques, il est impossible de définir un parcours moyen pour les
photons. On a recours à une loi d'atténuation qui s'exprime par :

I = I 0 e− x

I0 I(x)

où :
• I0 est le flux de photons incidents ;
• x l'épaisseur (en g/cm2) de la lame absorbante ;
• I(x) est le flux de photons qui émergent de la lame ;
NA
• µ est le coefficient d'atténuation massique (en cm2/g),  = ∑
A i i
;

• A est la masse atomique du milieu ;

•  i est la section efficace du processus i par atome ;
• NA est le nombre d'Avogadro.

1
On peut remarquer que : = où  est le libre parcours moyen des photons.


18
 Johann Collot Cours de physique expérimentale des hautes énergies de l'Université de Grenoble
Interaction des particules avec la matière
[Link] Année : 2011

4.1 Effet photo-électrique :

γ
e-

atome

Ce processus est permis lorsqu'il se produit sur des électrons atomiques. Il ne peut se produire sur des
électrons libres. Pour des photons suffisamment énergétiques, l'électron qui participe au processus est
pris sur la couche K à cause de sa proximité avec le noyau. Le site laissé vacant peut être repeuplé par
des électrons des couches externes donnant ainsi lieu à l'émission de rayons X ou d'électrons Auger
lorsque ces rayons sont « auto-absorbés ».

Ee = EK - 2 EL e-
-EL

-EK Ex = EK - EL si Ex > EL

Réorganisation interne conduisant à l'émission d'un électron Auger

Il n'existe pas de forme analytique simple de la section efficace qui pourrait être employée sur tout le
domaine en Z et en énergie du photon incident.
Dans le domaine non-relativiste, mais suffisamment loin des seuils de couches (EK , EL ... ) , celle-ci
s'exprime par :
K 32 21 4 5 e
 = 7   Z Th
photo (par atome)


19
 Johann Collot Cours de physique expérimentale des hautes énergies de l'Université de Grenoble
Interaction des particules avec la matière
[Link] Année : 2011

où :
• α est la constante de structure fine (1/137)
• Z est le numéro atomique du milieu
• =E  / me

•  eTh est la section efficace de diffusion élastique de Thomson d'un photon sur un
e 8 2 
électron , Th =  r e avec re = , rayon classique de l'électron .
3 me
Au voisinage du seuil de la couche K ( ~ EL ) , la section efficace présente un accident (dent de scie) de
passage à la couche L.
K 24 5
Pour des énergies plus grandes, c'est-à-dire pour ε >>1, on a plutôt :  photo = 4 r Ze .

Dans les formules précédentes, la dépendance en Z a été approchée par Z5 . En fait l'exposant de Z varie
de 4 à 5 pour 0,1 MeV < Eγ < 5 MeV .
À basse énergie (Eγ < 100 keV), c'est l'effet photo-électrique qui domine dans la section efficace totale
d'interaction des photons .

4.2 Effet Compton :

Ee

Eγ θε

θγ
Eγ'

Ce processus est la diffusion élastique d'un photon sur un électron du cortège atomique considéré
comme quasi-libre. La conservation de la norme des quadri-vecteurs permet d'obtenir :
E  m e =E e E '  − pe E '  cos e   .
D'autre part la conservation de la quantité de mouvement analysée sur la direction de propagation du
photon diffusé donne :
E  cos =E '   pe E '  cos  e  .
En combinant ces deux relations on obtient :

20
 Johann Collot Cours de physique expérimentale des hautes énergies de l'Université de Grenoble
Interaction des particules avec la matière
[Link] Année : 2011
E' me
= avec E e = E me−E '
E E e − pe cose
E' me
=
E E me−E  '−E  cos E '
1 1 E
= = avec  =
E 11−cos  me
1 1−cos  
me
E  ' est minimale pour = , c'est-à-dire dans le cas de la rétro-diffusion :
E ' min. 1
=
E 12 
L'angle de diffusion de l'électron peut être obtenu par l'équation suivante :

cotg e = 1 tg  
2

La section efficace totale est fournie par la formule de Klein-Nishina :

1 2 1 1 1 13 
 ec = 2  r2e   { − ln12 }  ln12  −  ( par électron )
 2
12   2 12 2
À basse énergie pour ε << 1 , l'expression se réduit à :
8 2
ec = eTh 1−2  avec eTh = r
3 e
qui est la section efficace de Thomson ( formule classique de diffusion sur un électron ).
À haute énergie (ε > > 1), on obtient :

e  r 2e 1
 =c  ln2  
 2
La section efficace différentielle ( pour tous les régimes ) est donnée par :
d 2 1 2 1cos2  2 1−cos  2
= re      1 
d 11−cos   2 1cos 2  1 1−cos  
Pour obtenir la section efficace par atome, il faut multiplier les expressions précédentes par Z :
atom.
c = Z  ec

4.3 Production de paires e+e- :

21
 Johann Collot Cours de physique expérimentale des hautes énergies de l'Université de Grenoble
Interaction des particules avec la matière
[Link] Année : 2011

Εγ Εγ e-
e-
e+
e+
atome
atome

Dans le champ du noyau Dans le champ d'un électron

2me2
E  2me  E  4 me
mN
mN ≫me ⇒ E 2 me
La production de paires dans le champ électrique de l'électron est fortement supprimée (champ plus
faible que le champ nucléaire) .
1 atom. 2 2 7 105
Si 1≪  ≪ ⇒  paire = 4 r e Z  ln 2− 
 Z 1/ 3 9 54
1 atom. 2 2 7 183 1
Si ≫ ⇒  paire = 4  r e Z  ln  1/ 3 − 
 Z1 /3 9 Z 54

atm. 7 2 2 183
À haute énergie, si l'on néglige le petit terme 1/54, on obtient :  paire ≃ 4  r e Z ln  1 /3  , ou encore
9 Z
:
atom. 7 A 1
 paire ≃ où X0 est la longueur de radiation du milieu .
9 N X0
La production de paires sur les électrons atomiques et proportionnelle à Z, si bien que la section efficace
totale varie en Z(Z+1) . Ceci est inclus dans la définition de X0 .

4.4 Section efficace totale d'absorption :

Dans la diffusion Compton, le photon incident n'est pas totalement absorbé. On définit une section
efficace de diffusion en énergie par :
E '
 atom.
cs =  atom.
c , où E'γ et Eγ sont respectivement les énergies des photons diffusés et incidents
E
atom. atom. atom.
De même, on peut définir une section efficace Compton d'absorption par : ca = c − cs .
Les différents coefficients massiques (en cm2/g) sont alors :
N atom. N atom.
cs =  ; ca =
 ; c = cs ca
A cs A ca
N
 paire = atom.
paire ;  photo =  photo
A
a =  photo paireca coefficient total d'absorption massique
 =  photo pairec coefficient total d'atténuation massique

22
 Johann Collot Cours de physique expérimentale des hautes énergies de l'Université de Grenoble
Interaction des particules avec la matière
[Link] Année : 2011

5 Interaction forte des hadrons :

L'interaction forte des hadrons joue un rôle important pour comprendre leur détection à haute énergie
dans les calorimètres hadroniques. Pour des énergies supérieures à quelques GeV, la section efficace
totale de diffusion proton-proton dépasse 50 mb. La section efficace totale est donnée par :
 tot =  élastique  inélastique .

On définit la longueur d'absorption nucléaire a :

A
a = g/cm2
N inélastique

23
 Johann Collot Cours de physique expérimentale des hautes énergies de l'Université de Grenoble
Interaction des particules avec la matière
[Link] Année : 2011

a est parfois désignée sous le nom de longueur d'interaction nucléaire dans la littérature. De même
on trouve dans la littérature la longueur de collision nucléaire T qui correspond à :
A
T = g/cm2 .
N tot

Pour le calcul de a et T , on fait l'hypothèse que tot et inélastique ne dépendent ni de l'énergie

incidente ni de l'espèce des hadrons considérés. En toute rigueur, cela n'est pas tout à fait correct.
Pour Z  6 , a et T sont très nettement supérieures à la longueur de radiation. Pour une même
énergie incidente, il faudra beaucoup plus de matière pour stopper un hadron qu'un électron (vrai à haute
énergie).

matériau Z A σtot (barn) σinel (barn) λT (g/cm2) λa (g/cm2)

hydrogène 1 1,01 0,0387 0,033 43,3 50,8
hélium 2 4,0 0,133 0,102 49,9 65,1
béryllium 4 9,01 0,268 0,199 55,8 75,2
carbone 6 12,01 0,331 0,231 60,2 86,3
azote 7 14,01 0,379 0,265 61,4 87,8
oxygène 8 16,0 0,420 0,292 63,2 91,0
aluminium 13 26,98 0,634 0,421 70,6 106,4
silicium 14 28,09 0,660 0,440 70,6 106,0
fer 26 55,85 1,120 0,703 82,8 131,9
cuivre 29 63,55 1,232 0,782 85,6 134,9
tungstène 74 183,85 2,767 1,65 110,3 185
plomb 82 207,19 2,960 1,77 116,2 194
uranium 92 238,03 3,378 1,98 117,0 199

24
Johann Collot Cours de physique expérimentale des hautes énergies de l'Université de Grenoble
Interaction des particules avec la matière
[Link] Année : 2011

6 Pour en savoir plus :

• [Référence 1] : Nuclei and particles , Émilio Segrè , W.A. Benjamin, inc.
• [Référence 2] : Stopping powers and ranges for protons and alpha particles (ICRU Report 49,
1993) Library of Congress US-Cataloging-in-Publication Data
• [Référence 3] : GEANT - Detector description and simulation tool - CERN program library - phys
431
• Particle Detectors , Claus Grupen , Cambridge monographs on particle physics
• Detectors for Particle radiation , Konrad Kleinknecht , Cambridge University Press
• Radiation detection and measurement, G.F. Knoll , J. Wiley & Sons
• Single Particle Detection and Measurement, R. Gilmore, Taylor & Francis
• Review of Particle Physics : [Link]

25