Interaction of Particles with Matter (Lecture 1)
Total Page:16
File Type:pdf, Size:1020Kb
Interaction of particles with matter (lecture 1) Johann Collot [email protected] http://lpsc.in2p3.fr/collot UdG 1/66 A brief review of a few typical situations is going to greatly simplify the subject. Mean free path of a particle, i.e. average distance travelled between two consecutive collisions in matter : 1 λ= σ n where : total interaction cross-section of the particle n number of scattering centers per unit volume N A for a monoatomic element of molar mass M and example : n= M specific mass ρ . N A Avogadro number Electromagnetic interaction : 1 m (charged particles) Strong interaction : 1cm (neutrons ....) Weak interaction : 1015 m≃0,1 light year (neutrinos) A practical signal ( >100 interactions or hits ) can only come from electromagnetic interaction Particle detection proceeds in two steps : 1) primary interaction 2) charged particle interaction producing the signals Johann Collot [email protected] http://lpsc.in2p3.fr/collot UdG 2/66 typical examples : photon detection Compton scattering 2 1 Signal is induced by electrons e- e- 3 e+ Pair production Johann Collot [email protected] http://lpsc.in2p3.fr/collot UdG 3/66 neutral pion detection : A π0 decays into two photons with a mean lifetime of 8.5 10-17 s. e+ 1 - e 0 e- 2 e+ Johann Collot [email protected] http://lpsc.in2p3.fr/collot UdG 4/66 neutrino detection : A 2800 MW nuclear power station produces 130 MW of neutrinos ! A detector of 1 m3 located 20 m away from the reactor core can detect 100 neutrinos/h. e liquid scintillator ( CH + 6L ) e nuclear reactor x i e e e p 6 e L t i n e e ~0 p n e e Charged particles produce light 6 nthLi t4,8 MeV in the target scintillator Johann Collot [email protected] http://lpsc.in2p3.fr/collot UdG 5/66 Interaction of charged particles with matter For heavy particles ionization and excitation are the dominant processes producing energy loss. Particle P of Z charge state * Excitation : PZ atom atom*PZ followed by : atom atom Ionization : P(Z )+atom→atom + +e-+P(Z ) Ionization + excitation : P(Z )+atom→atom *++e-+P( Z ) Johann Collot [email protected] http://lpsc.in2p3.fr/collot UdG 6/66 max Te , Maximal kinetic energy transferred to an ionized electron : m m0 e v v = c p= m0 v atom γ=(1−β2)−1/ 2 E=( p2+m2)−1/2 Z 0 β=(1−γ−2)1/2 Hypothesis : V > <ve> = Z α c , speed of deepest atomic orbit electrons where α is the fine structure constant : α = 1/137 . 2 2 max max 2me One may show (exercise) that : T e =E e −me= 2 (In natural units , c=ħ=1) E CM /m0 2 2 1/2 where : ECM = m0me2me E total energy in center-of-mass frame Johann Collot [email protected] http://lpsc.in2p3.fr/collot UdG 7/66 Two cases : m0 >> me , i.e. the incoming particle is not an electron and if its energy is not too big m2 m2 2m E 2 0 e e E=γ m (ECM /m0) = ( 2 + 2 + 2 )≃1 with 0 m0 m0 m0 2 γ me ≪1 proton Ep < 50 GeV , muon Eμ < 500 MeV (medium energy range) m0 max max 2 2 then : T e =Ee −me=2me β γ m0 = me the incoming particle is an electron T max=(E−m ) due to undistinguishibility of electrons, max transferable e e energy = Tmax / 2 If the incoming particle is not an electron then in practice m0 >> me . Johann Collot [email protected] http://lpsc.in2p3.fr/collot UdG 8/66 Stopping power of heavy particles by excitation and ionization in matter. Average energy loss by a charged particle (other than an electron) in matter. Bethe and Bloch formula (see Nuclei and particles, Émilio Segré, W.A. Benjamin ; Principles of Radiation Interaction in Matter and Detection, C. Leroy and P.G. Rancoita, World Scientific ; Introduction to experimental particle physics, R. Fernow ) Stopping power or mean specific energy loss charge of incoming particle density effect correction Z of medium at high energy 2 2 2 max dE MeV 0.3071 z Z 1 2me T e 2 C e − [ ] = ln − − − dx g/cm2 A g mol−1 2 2 I 2 2 Z Atomic shell correction at low mean excitation energy energy Atomic mass of medium (not covered in this lecture, see Leroy & Rancoita) Surface mass density of medium dx = ρ dl (or mass thickness of medium) Johann Collot [email protected] http://lpsc.in2p3.fr/collot UdG 9/66 show that : (β γ)2 β2= 1+ (β γ)2 Johann Collot [email protected] http://lpsc.in2p3.fr/collot UdG 10/66 Particle identification in Alice TPC Johann Collot [email protected] http://lpsc.in2p3.fr/collot UdG 11/66 Johann Collot [email protected] http://lpsc.in2p3.fr/collot UdG 12/66 few remarks : dE 5 3 − ~− / non relativistic particles - for βγ < 1 : dx dE − - for βγ ~ 3-4 : dx is minimal over a large energy plateau . A particle in this state is called a minimum ionizing particle (MIP) dE MIP MeV In media composed of light elements : − ≃2 - 2 dx g cm dE - for βγ > 4 : relativistic increase of − as ln dx wich is tempered by -δ/2 correction. - I : mean excitation and ionization energy , I = 15 eV for atomic H and 19.2 eV for H2 I = 41.8 for He 0.9 I = 15 Z eV for Z > 2 2me max 2 2 At medium energy : ≪1 T e =2me m0 2 2 2 dE MeV 0.3071 z Z 2me 2 C e − [ ]= ⋅ [ln − − − ] dx g/cm2 A g 2 I 2 Z Johann Collot [email protected] http://lpsc.in2p3.fr/collot UdG 13/66 Density effect correction When energy increases, stopping power decreases to a minimum (1/β2 dependance) and then starts rising again due to logarithmic term. In fact, the max. transverse electric field increases as γ but its influence is screened by nearby atoms beyond a distance of 70 (A(g)/ρ(g)Z)1/2 Å (shown by Bohr). This density effect tempers the relativistic rise. Studies have be carried-out by Sernheimer, Peierls, Berger & Seltzer (see Leroy & Rancoita). The density correction effect term, δ is given by : md 2 S1 S β γ S0 S1 1 10 for 0 for 10 < β γ< 10 : δ=2ln(β γ)+ C+ a ln( ) β γ< 10 : δ=δ0( S ) β γ 10 0 [ ln (10) ] S I for β γ> 10 1 : δ=2ln(β γ)+ C where : C=−2ln( )−1 h νp nr c² with ν = e in which n is the density of electrons and r the classical radius of e- : r = 2.82 fm p √ π e e S show that at very high energy : β γ> 10 1 ρ(g/cm³) h ν ≃28.7 Z eV max p A(g) dE MeV z2 Z 2m T √ −( )[ ]=0.3071 [ln( e e )−1] dx 2 2 g/cm 2.A(g) (h νp) Johann Collot [email protected] http://lpsc.in2p3.fr/collot UdG 14/66 Johann Collot [email protected] http://lpsc.in2p3.fr/collot UdG 15/66 Restricted energy loss knock-on electron (delta ray) generated by a 180 GeV muon as observed by the experiment GridPix at CERN SPS. High energy transfers generate delta rays that may espace the detector if it is too thin. So average energy deposits are very often much smaler thant predicted by B&B. If T0 is the average maximal delta ray energy that can be absorbed in the detecting medium, a better estimate of the average deposited energy is given by : 2 2 dE MeV 0.3071 z2 Z 1 2m β γ T β2 T δ(γβ) C −( )[ ] = ( ln( e 0 )− (1+ 0 )− − e ) dx 2 −1 2 2 2 2 2 2 2 Z g/cm A (g mol ) β I 2me β γ S At extremely high energies, when β γ> 10 1 , stopping power reaches a constant called Fermi plateau. 2 dE MeV z Z 2me T 0 −( )[ 2 ]=0.3071 ln( 2 ) dx g/cm 2.A(g) (h νp) Johann Collot [email protected] http://lpsc.in2p3.fr/collot UdG 16/66 Fermi plateau measured in silicon Johann Collot [email protected] http://lpsc.in2p3.fr/collot UdG 17/66 Delta rays (secondary electrons) The differential probability to generate a delta ray of kinetic energy T is given by : 2 dw(T , E) z Z F (T ) − − =0.3071 MeV 1 cm2 g 1 dTdx 2.A(g)β2 T 2 F(T) is a spin-dependent factor. 2 T F (T )=F 0(T )=(1−β ) for spin-0 particles T max 1 T 2 F (T )=F (T )=F (T )+ ( ) for spin-1/2 particles 1/2 0 2 E 2 1 T me 1 T 1 T me F (T )=F 1(T )=F 0(T )(1+ 2 )+ ( ) (1+ 2 ) for spin-1 particles 3 m0 3 E 2 m0 Johann Collot [email protected] http://lpsc.in2p3.fr/collot UdG 18/66 Delta rays (secondary electrons) 2 For T << Tmax and T << m0 / me , 2 dw(T , E) z Z 1 −1 2 −1 =0.3071 MeV cm g dTdx 2.A(g)β2 T 2 This allows to compute an approximate probability to generate a delta ray of kinetic energy greater than Ts in a thin absorber of mass thickness x : z2 Z 1 w(Ts , E , x)≃0.3071 x 2 show this expression 2.A(g)β T s Johann Collot [email protected] http://lpsc.in2p3.fr/collot UdG 19/66 Energy straggling distribution So far, only the average energy loss has been considered.