Fluka advanced calculations on stopping power and multiple Coulomb scattering
|
|
- Maximilian Haynes
- 6 years ago
- Views:
Transcription
1 Fluka advanced calculations on stopping power and multiple Coulomb scattering Andrea Fontana INFN Sezione di Pavia
2 Outline Stopping power Ionization potential Range calculation: which range? Fluka options Multiple Coulomb scattering Molière theory Nuclear interactions Fluka options Tools for analysis USRBIN readout and plot with ROOT USRBIN dynamic scan
3 Case study This study is performed in a simplified geometry, consisting of a water phantom in water (xx30 cm 3 ), hit by a monochromatic beam of protons with different energies. Geometry setup Andrea Fontana INFN Sezione di Pavia 3/45
4 Case study This study is performed in a simplified geometry, consisting of a water phantom in water (xx30 cm 3 ), hit by a monochromatic beam of protons with different energies. Gaussian (σ=0.6 cm) beam profile (lin) Andrea Fontana INFN Sezione di Pavia 4/45
5 Case study This study is performed in a simplified geometry, consisting of a water phantom in water (xx30 cm 3 ), hit by a monochromatic beam of protons with different energies. Gaussian (σ=0.6 cm) beam profile (log) Andrea Fontana INFN Sezione di Pavia 5/45
6 Case study This study is performed in a simplified geometry, consisting of a water phantom in water (xx30 cm 3 ), hit by a monochromatic beam of protons with different energies. Pencil beam profile (log) Andrea Fontana INFN Sezione di Pavia 6/45
7 Stopping power de [ = Kz 2 Z dx A ϱ 1 1 β 2 2 ln 2mec2 β 2 γ 2 T max β 2 δ(βγ) ] I 2 2 Andrea Fontana INFN Sezione di Pavia 7/45
8 Ionization potential A material is completely described by a single number, the mean excitation potential I. Felix Bloch approximation (1933): I = (ev )Z whit Z atomic number of the atoms of the material. Modern tables (e.g. ICRU Report 37 and 49) yield better results than this formula. For compounds: ( ) Z Mj ZM lni = Σw j lni j / A Mj A M with ZM A M = Σw j Z Mj A Mj Z M and A M are charge number and mass number of the target. Andrea Fontana INFN Sezione di Pavia 8/45
9 "What is the ionization potential of water?" FLUKA default value is I=75 ev, but it is usually redefined. Its exact knowledge is difficult, both theoretically and experimentally. MAT-PROP 77. WATER Andrea Fontana INFN Sezione di Pavia 9/45
10 Ionization potential Bragg peak and water ionization potential for 200 MeV protons energy (a.u.) I=75 ev z (cm) Andrea Fontana INFN Sezione di Pavia /45
11 Ionization potential Bragg peak and water ionization potential for 200 MeV protons energy (a.u.) I=75 ev I=77 ev z (cm) Andrea Fontana INFN Sezione di Pavia 11/45
12 Ionization potential Bragg peak and water ionization potential for 200 MeV protons energy (a.u.) I=75 ev I=77 ev I=79 ev z (cm) Andrea Fontana INFN Sezione di Pavia 12/45
13 Range calculation: which range? Mean range Projected range CSDA range 90% distal fall-off range analytical range Andrea Fontana INFN Sezione di Pavia 13/45
14 CSDA range The CSDA (Continuous Slowing Down Approximation) ranges for protons of different energies in water is investigated. The CSDA range is calculated by integrating the inverse of the linear stopping power of the proton from zero to the initial energy: no energy straggling R CSDA = no delta ray production no nuclear interactions L no Coulomb multiple scattering 0 dx = 0 E 1 de/dx de Energy (MeV) Ranges CSDA (mm) FLUKA NIST This corresponds essentially to follow the primary particle ionization energy loss only. Andrea Fontana INFN Sezione di Pavia 14/45
15 Fluka options EMFCUT: transport and production thrs. for e +, e and γ EMFCUT -3E-05 3.E-6 1. PROD-CUT EMFCUT -3E-05 3.E PART-THR: all partice thrs. kev, neutron thr MeV PART-THR -1E PART-THR NEUTRON 0.0 THRESHOL: infinite hadron elastic scattering and inelastic interactions THRESHOL Andrea Fontana INFN Sezione di Pavia 15/45
16 DELTARAY: infinite delta ray thr. DELTARAY -0.E PRINT IONFLUCT: switch off hadrons and leptons ionization straggling IONFLUCT FLUKAFIX: fraction of the kinetic energy to be lost in a step FLUKAFIX STEPSIZE: set maximum step size STEPSIZE 0.01 Andrea Fontana INFN Sezione di Pavia 16/45
17 CSDA range Andrea Fontana INFN Sezione di Pavia 17/45
18 Analytical stopping power? Øverås has found a simple empirical relation between p(x) 2 β(x) 2 and normalized residual range (see also Gottschalk papers): ) k ( p(x) 2 β(x) 2 = p(0) 2 β(0) 2 1 x R ( p(x) 2 β(x) 2 p(0) 2 β(0) 2 1 x ) R Valid for many materials of density ρ (g/cm 3 ) and radiation length X 0 (cm): k = exp( 0.09ρX 0) Andrea Fontana INFN Sezione di Pavia 18/45
19 Analytical range calculation? For an incident proton of kinetic energy E(MeV), Ulmer has published a very accurate formula for a medium of density ρ (g/cm 3 ) is: R(cm) = 1 A M ρ Z M N α ni p i E n n=1 order α i p i e e e e where I is the mean ionization potential of the medium (in ev) and Z M and A M are the effective atomic and mass numbers. For water and protons of incident kinetic energy E < 300 MeV, a sum with N = 1..4 gives results accurate to a more than 0.5%. Andrea Fontana INFN Sezione di Pavia 19/45
20 Multiple Coulomb scattering Molière has found a probability distribution function for θ that is a solution of the transport equation and that well describes the experimental data. The expression of this solution reads as: f(θ)θdθ = f M (θ)d(cosθ)dφ/2π and, by using the approximation d(cosθ) = sinθ dθ θ dθ, can be written as the sum of three terms: f M (θ) = 1 [ f 0 (θ ) + f 1 (θ ) 2θM 2 B ) + f 2 (θ ] ) B 2 where θ M is the characteristic multiple scattering angle of the target, θ = θ/( 2θ M ) is the reduced angle and B is related to the logarithm of the effective number of collisions in the target. Andrea Fontana INFN Sezione di Pavia 20/45
21 Molière theory Molière theory depends on two parameters: χ 2 c and χ 2 α. The parameter χ 2 c is related to the RMS of the scattering angle: χ 2 c = Z 2 z 2 x 1 A p 2 β, 2 The parameter χ 2 α is related to the screening of the Coulomb potential: ( ) µ 2 = z2 Z 2 χ 2 α = µ 2 χ β 2 0 ( χ 2 ħ Z 1/3 0 = p (cm) ) 2 Andrea Fontana INFN Sezione di Pavia 21/45
22 Multi-gaussian approximation f 0 (x) = 2e x f 1 (x) = 2e x (x 1)[Ēi(x) lnx] 2(1 2e x ) f 2 (x) = e ([ψ x 2 (2) + ψ(2)](x 2 4x + 2) + 4 y 3 dy[lny/(1 y) ψ(2)] 0 ) [(1 y 2 )e xy 1 (x 2)y (x 2 /2 2x + 1)y 2 ] whit x = θ 2, ψ(n) = d{lnγ(n + 1)}/ dn the Digamma function and Ei(x) the Exponential Integral function. The first term is a gaussian function that constitutes the "core" of the distribution, while the extra terms account for the tails of the distribution that are non-gaussian: they can be evaluated by numerical integration or by using mathematical tables. 1 Andrea Fontana INFN Sezione di Pavia 22/45
23 Nuclear interactions Andrea Fontana INFN Sezione di Pavia 23/45
24 Nuclear interactions In the case of protons in water, the total nuclear cross-section for proton-nucleus interactions is often calculated only taking into account the p 16 O reactions: p O n F p O p + n O p O p + p N p O α N p O d O n O p N n O p + n N In general terms: t(x) f(x) = W pf M (x)+(1 W p) + t(x)dx f M (x) = Molière electromagnetic theory t(x) = parametrization of nuclear tails W p = fraction of events without nuclear interactions Andrea Fontana INFN Sezione di Pavia 24/45
25 Analytical description of nuclear interactions Ulmer formula for the determination of the percentage W p of events without nuclear interactions (primary protons), as a function of the traversed thickness, for protons of incident kinetic energy E and range R in water, at a certain water thickness x: [ W p = 1 ( ) E f ] [ ( )] Eth x R x erf, (1) 2 m R τ where erf is the error function, f = 1.032, m is the proton mass in MeV, E th = 7 MeV is the 16 O threshold energy of the Coulomb barrier. The parameter τ takes into account the range variation due to the straggling along the beam path and can be parametrized as: τ = R t, where t = { if R 1cm if R < 1cm (2) Andrea Fontana INFN Sezione di Pavia 25/45
26 Decrease of the fluence of primary protons in water Andrea Fontana INFN Sezione di Pavia 26/45
27 Secondaries production protons Andrea Fontana INFN Sezione di Pavia 27/45
28 Secondaries production neutrons Andrea Fontana INFN Sezione di Pavia 28/45
29 Secondaries production deuterons Andrea Fontana INFN Sezione di Pavia 29/45
30 Secondaries production alphas Andrea Fontana INFN Sezione di Pavia 30/45
31 Fluka options Fluka is the only MC code firmly based on the full Molière theory for MCS. MULSOPT: to define tracking conditions with many sub-options (Fano correction, Molière screening angle, spin relativistic corrections, finite nuclear size correction, single scattering at boundaries or very low energy). To completely suppress MCS: MULSOPT MCSTHRESh to specify accuracy for MCS of heavy charged particles, depending on DEFAULTS. Andrea Fontana INFN Sezione di Pavia 31/45
32 Fluka options Lateral profile for 200 MeV protons at z=25 cm energy (a.u.) -3 default z (cm) Andrea Fontana INFN Sezione di Pavia 32/45
33 Fluka options Lateral profile for 200 MeV protons at z=25 cm energy (a.u.) -3 default precision z (cm) Andrea Fontana INFN Sezione di Pavia 33/45
34 Fluka options Lateral profile for 200 MeV protons at z=25 cm energy (a.u.) -3 default precision optimize z (cm) Andrea Fontana INFN Sezione di Pavia 34/45
35 Fluka options Lateral profile for 200 MeV protons at z=25 cm energy (a.u.) -3 default precision optimize fano off z (cm) Andrea Fontana INFN Sezione di Pavia 35/45
36 Fluka options Lateral profile for 200 MeV protons at z=25 cm energy (a.u.) -3 default precision optimize fano off mlsh off z (cm) Andrea Fontana INFN Sezione di Pavia 36/45
37 Fluka options Lateral profile for 200 MeV protons at z=25 cm energy (a.u.) -3 default precision optimize fano off mlsh off born 1st z (cm) Andrea Fontana INFN Sezione di Pavia 37/45
38 Fluka options Lateral profile for 200 MeV protons at z=25 cm energy (a.u.) -3 default precision optimize fano off mlsh off born 1st born 2nd z (cm) Andrea Fontana INFN Sezione di Pavia 38/45
39 Fluka options Lateral profile for 200 MeV protons at z=25 cm energy (a.u.) -3 default precision optimize fano off mlsh off born 1st born 2nd size 1st z (cm) Andrea Fontana INFN Sezione di Pavia 39/45
40 Fluka options Lateral profile for 200 MeV protons at z=25 cm energy (a.u.) -3 default precision optimize fano off mlsh off born 1st born 2nd size 1st size 2nd z (cm) Andrea Fontana INFN Sezione di Pavia 40/45
41 Fluka options Lateral profile for 200 MeV protons at z=25 cm energy (a.u.) -3 default precision optimize fano off mlsh off born 1st born 2nd size 1st size 2nd off z (cm) Andrea Fontana INFN Sezione di Pavia 41/45
42 Tools for analysis USRBIN readout and plot with ROOT macro plotbrg.c : to plot 2D USRBIN ASCII output and project longitudinal profile. macro plotmcs.c : to plot 2D USRBIN ASCII output and project trasversal profile. USRBIN dynamic scan macro scan.c : to quickly inspect 2D USRBIN ASCII output and perform fits. Other files Fluka input brg.inp and brgcsda.inp : for Bragg peak and range studies. Fluka input mcs.inp : for Multiple Coulomb scattering studies. Andrea Fontana INFN Sezione di Pavia 42/45
43 References [1] J. R. Sabin et al., Advances in Quantum Chemistry, Vol. 65, 2013 [2] H. Øverås, CERN Yellow Report No , [3] B. Gottschalk, Med. Phys. 37(20)352 [4] R. Fruhwirth et al., NIM A 456(2000)369 [5] V.L. Highland, NIM 129(1975)497 [6] G.B. Arfken, Mathematical Methods for Physicists, Academic Press 2005 [7] W. Ulmer, Rad. Phys. and Chem. 76(2007)89 [8] W. Ulmer et al., Rad. Phys. and Chem. 80(2011)378 Andrea Fontana INFN Sezione di Pavia 43/45
44 Questions? Contact:
45 Acknowledgements Fluka team: A. Ferrari, P. Sala, A. Mairani et al. Pavia Fluka group: A. Rotondi, G. Magro, M. Mori, E. Bellinzona and A. Embriaco
Ionization and Transport. FLUKA Beginners Course
Ionization and Transport FLUKA Beginners Course General settings Interactions of leptons/photons Photon interactions Photoelectric Compton Rayleigh Pair production Photonuclear Photomuon production Electron/positron
More informationPhysics of particles. H. Paganetti PhD Massachusetts General Hospital & Harvard Medical School
Physics of particles H. Paganetti PhD Massachusetts General Hospital & Harvard Medical School Introduction Dose The ideal dose distribution ideal Dose: Energy deposited Energy/Mass Depth [J/kg] [Gy] Introduction
More informationPhysics of Particle Beams. Hsiao-Ming Lu, Ph.D., Jay Flanz, Ph.D., Harald Paganetti, Ph.D. Massachusetts General Hospital Harvard Medical School
Physics of Particle Beams Hsiao-Ming Lu, Ph.D., Jay Flanz, Ph.D., Harald Paganetti, Ph.D. Massachusetts General Hospital Harvard Medical School PTCOG 53 Education Session, Shanghai, 2014 Dose External
More informationIonization and Transport. FLUKA Beginners Course
Ionization and Transport FLUKA Beginners Course General settings Interactions of leptons/photons Photon interactions Photoelectric Compton Rayleigh Pair production Photonuclear Photomuon production Electron/positron
More information2. Passage of Radiation Through Matter
2. Passage of Radiation Through Matter Passage of Radiation Through Matter: Contents Energy Loss of Heavy Charged Particles by Atomic Collision (addendum) Cherenkov Radiation Energy loss of Electrons and
More informationCHARGED PARTICLE INTERACTIONS
CHARGED PARTICLE INTERACTIONS Background Charged Particles Heavy charged particles Charged particles with Mass > m e α, proton, deuteron, heavy ion (e.g., C +, Fe + ), fission fragment, muon, etc. α is
More informationParticle-Matter Interactions
Particle-Matter Interactions to best detect radiations and particles we must know how they behave inside the materials 8/30/2010 PHYS6314 Prof. Lou 1 Stable Particles Visible to a Detector Hadrons (Baryon/Meson)
More informationOutline. Chapter 6 The Basic Interactions between Photons and Charged Particles with Matter. Photon interactions. Photoelectric effect
Chapter 6 The Basic Interactions between Photons and Charged Particles with Matter Radiation Dosimetry I Text: H.E Johns and J.R. Cunningham, The physics of radiology, 4 th ed. http://www.utoledo.edu/med/depts/radther
More informationIl picco di Bragg. G. Battistoni INFN Milano. 08/06/2015 G. Battistoni
Il picco di Bragg G. Battistoni INFN Milano 08/06/015 G. Battistoni 1 Φ(z) The physics of Bragg Peak 180 MeV proton in water Longitudinal profile: Transversel profile: Φ(z,x) dominated by interaction with
More informationInteractions of Particulate Radiation with Matter. Purpose. Importance of particulate interactions
Interactions of Particulate Radiation with Matter George Starkschall, Ph.D. Department of Radiation Physics U.T. M.D. Anderson Cancer Center Purpose To describe the various mechanisms by which particulate
More informationThe interaction of radiation with matter
Basic Detection Techniques 2009-2010 http://www.astro.rug.nl/~peletier/detectiontechniques.html Detection of energetic particles and gamma rays The interaction of radiation with matter Peter Dendooven
More informationMonte Carlo radiation transport codes
Monte Carlo radiation transport codes How do they work? Michel Maire (Lapp/Annecy) 23/05/2007 introduction to Monte Carlo radiation transport codes 1 Decay in flight (1) An unstable particle have a time
More informationPHYS 352. Charged Particle Interactions with Matter. Intro: Cross Section. dn s. = F dω
PHYS 352 Charged Particle Interactions with Matter Intro: Cross Section cross section σ describes the probability for an interaction as an area flux F number of particles per unit area per unit time dσ
More informationParticle Interactions in Detectors
Particle Interactions in Detectors Dr Peter R Hobson C.Phys M.Inst.P. Department of Electronic and Computer Engineering Brunel University, Uxbridge Peter.Hobson@brunel.ac.uk http://www.brunel.ac.uk/~eestprh/
More informationIII. Energy Deposition in the Detector and Spectrum Formation
1 III. Energy Deposition in the Detector and Spectrum Formation a) charged particles Bethe-Bloch formula de 4πq 4 z2 e 2m v = NZ ( ) dx m v ln ln 1 0 2 β β I 0 2 2 2 z, v: atomic number and velocity of
More informationPHYS 571 Radiation Physics
PHYS 571 Radiation Physics Prof. Gocha Khelashvili http://blackboard.iit.edu login Interaction of Electrons with Matter The Plan Interactions of Electrons with Matter Energy-Loss Mechanism Collisional
More informationEEE4101F / EEE4103F Radiation Interactions & Detection
EEE4101F / EEE4103F Radiation Interactions & Detection 1. Interaction of Radiation with Matter Dr. Steve Peterson 5.14 RW James Department of Physics University of Cape Town steve.peterson@uct.ac.za March
More informationNuclear Physics and Astrophysics
Nuclear Physics and Astrophysics PHY-30 Dr. E. Rizvi Lecture 4 - Detectors Binding Energy Nuclear mass MN less than sum of nucleon masses Shows nucleus is a bound (lower energy) state for this configuration
More informationInteraction of Particles with Matter
Chapter 10 Interaction of Particles with Matter A scattering process at an experimental particle physics facility is called an event. Stable particles emerging from an event are identified and their momenta
More informationPHYS 5012 Radiation Physics and Dosimetry
Radiative PHYS 5012 Radiation Physics and Dosimetry Mean Tuesday 24 March 2009 Radiative Mean Radiative Mean Collisions between two particles involve a projectile and a target. Types of targets: whole
More informationBethe-Block. Stopping power of positive muons in copper vs βγ = p/mc. The slight dependence on M at highest energies through T max
Bethe-Block Stopping power of positive muons in copper vs βγ = p/mc. The slight dependence on M at highest energies through T max can be used for PID but typically de/dx depend only on β (given a particle
More informationDetectors in Nuclear Physics (40 hours)
Detectors in Nuclear Physics (40 hours) Silvia Leoni, Silvia.Leoni@mi.infn.it http://www.mi.infn.it/~sleoni Complemetary material: Lectures Notes on γ-spectroscopy LAB http://www.mi.infn.it/~bracco Application
More informationPhysics 663. Particle Physics Phenomenology. April 23, Physics 663, lecture 4 1
Physics 663 Particle Physics Phenomenology April 23, 2002 Physics 663, lecture 4 1 Detectors Interaction of Charged Particles and Radiation with Matter Ionization loss of charged particles Coulomb scattering
More informationMonte Carlo radiation transport codes
Monte Carlo radiation transport codes How do they work? Michel Maire (Lapp/Annecy) 16/09/2011 introduction to Monte Carlo radiation transport codes 1 Outline From simplest case to complete process : Decay
More informationUniversity of Oslo. Department of Physics. Interaction Between Ionizing Radiation And Matter, Part 2 Charged-Particles.
Interaction Between Ionizing Radiation And Matter, Part Charged-Particles Audun Sanderud Excitation / ionization Incoming charged particle interact with atom/molecule: Ionization Excitation Ion pair created
More informationDetectors in Nuclear Physics (48 hours)
Detectors in Nuclear Physics (48 hours) Silvia Leoni, Silvia.Leoni@mi.infn.it http://www.mi.infn.it/~sleoni Complemetary material: Lectures Notes on γ-spectroscopy LAB http://www.mi.infn.it/~bracco Application
More informationInteraction of Ionizing Radiation with Matter
Type of radiation charged particles photonen neutronen Uncharged particles Charged particles electrons (β - ) He 2+ (α), H + (p) D + (d) Recoil nuclides Fission fragments Interaction of ionizing radiation
More informationEmphasis on what happens to emitted particle (if no nuclear reaction and MEDIUM (i.e., atomic effects)
LECTURE 5: INTERACTION OF RADIATION WITH MATTER All radiation is detected through its interaction with matter! INTRODUCTION: What happens when radiation passes through matter? Emphasis on what happens
More informationInterazioni delle particelle cariche. G. Battistoni INFN Milano
Interazioni delle particelle cariche G. Battistoni INFN Milano Collisional de/dx of heavy charged particles A particle is heavy if m part >>m electron. The energy loss is due to collision with atomic electrons.
More informationFLUKA Manual and Basic input. 7th FLUKA Course NEA Paris, Sept.29-Oct.3, 2008
FLUKA Manual and Basic input 7th FLUKA Course NEA Paris, Sept.29-Oct.3, 2008 The FLUKA Manual in continuous development, just as the program more a User Guide than a Reference Manual (only a short summary
More informationHeavy charged particle passage through matter
Heavy charged particle passage through matter Peter H. Hansen University of Copenhagen Content Bohrs argument The Bethe-Bloch formula The Landau distribution Penetration range Biological effects Bohrs
More informationInteraction of Particles and Matter
MORE CHAPTER 11, #7 Interaction of Particles and Matter In this More section we will discuss briefly the main interactions of charged particles, neutrons, and photons with matter. Understanding these interactions
More informationChapter II: Interactions of ions with matter
Chapter II: Interactions of ions with matter 1 Trajectories of α particles of 5.5 MeV Source: SRIM www.srim.org 2 Incident proton on Al: Bohr model v=v 0 E p =0.025 MeV relativistic effect E p =938 MeV
More informationChapter 2 Radiation-Matter Interactions
Chapter 2 Radiation-Matter Interactions The behavior of radiation and matter as a function of energy governs the degradation of astrophysical information along the path and the characteristics of the detectors.
More informationThe next three lectures will address interactions of charged particles with matter. In today s lecture, we will talk about energy transfer through
The next three lectures will address interactions of charged particles with matter. In today s lecture, we will talk about energy transfer through the property known as stopping power. In the second lecture,
More informationInteractions of particles and radiation with matter
1 Interactions of particles and radiation with matter When the intervals, passages, connections, weights, impulses, collisions, movement, order, and position of the atoms interchange, so also must the
More informationNeutron Interactions Part I. Rebecca M. Howell, Ph.D. Radiation Physics Y2.5321
Neutron Interactions Part I Rebecca M. Howell, Ph.D. Radiation Physics rhowell@mdanderson.org Y2.5321 Why do we as Medical Physicists care about neutrons? Neutrons in Radiation Therapy Neutron Therapy
More informationChapter V: Interactions of neutrons with matter
Chapter V: Interactions of neutrons with matter 1 Content of the chapter Introduction Interaction processes Interaction cross sections Moderation and neutrons path For more details see «Physique des Réacteurs
More informationDetectors for High Energy Physics
Detectors for High Energy Physics Ingrid-Maria Gregor, DESY DESY Summer Student Program 2017 Hamburg July 26th/27th Disclaimer Particle Detectors are very complex, a lot of physics is behind the detection
More informationOutline. Physics of Charge Particle Motion. Physics of Charge Particle Motion 7/31/2014. Proton Therapy I: Basic Proton Therapy
Outline Proton Therapy I: Basic Proton Therapy Bijan Arjomandy, Ph.D. Narayan Sahoo, Ph.D. Mark Pankuch, Ph.D. Physics of charge particle motion Particle accelerators Proton interaction with matter Delivery
More informationRadiation Physics PHYS /251. Prof. Gocha Khelashvili
Radiation Physics PHYS 571-051/251 Prof. Gocha Khelashvili Interaction of Radiation with Matter: Heavy Charged Particles Directly and Indirectly Ionizing Radiation Classification of Indirectly Ionizing
More informationIonization Energy Loss of Charged Projectiles in Matter. Steve Ahlen Boston University
Ionization Energy Loss of Charged Projectiles in Matter Steve Ahlen Boston University Almost all particle detection and measurement techniques in high energy physics are based on the energy deposited by
More informationLecture 14 (11/1/06) Charged-Particle Interactions: Stopping Power, Collisions and Ionization
22.101 Applied Nuclear Physics (Fall 2006) Lecture 14 (11/1/06) Charged-Particle Interactions: Stopping Power, Collisions and Ionization References: R. D. Evans, The Atomic Nucleus (McGraw-Hill, New York,
More informationParticle Detectors. Summer Student Lectures 2010 Werner Riegler, CERN, History of Instrumentation History of Particle Physics
Particle Detectors Summer Student Lectures 2010 Werner Riegler, CERN, werner.riegler@cern.ch History of Instrumentation History of Particle Physics The Real World of Particles Interaction of Particles
More informationLET! (de / dx) 1 Gy= 1 J/kG 1Gy=100 rad. m(kg) dose rate
Basics of Radiation Dosimetry for the Physicist http://en.wikipedia.org/wiki/ionizing_radiation I. Ionizing radiation consists of subatomic particles or electromagnetic waves that ionize electrons along
More informationToday, I will present the first of two lectures on neutron interactions.
Today, I will present the first of two lectures on neutron interactions. I first need to acknowledge that these two lectures were based on lectures presented previously in Med Phys I by Dr Howell. 1 Before
More informationNeutrino detection. Kate Scholberg, Duke University International Neutrino Summer School Sao Paulo, Brazil, August 2015
Neutrino detection Kate Scholberg, Duke University International Neutrino Summer School Sao Paulo, Brazil, August 2015 Sources of wild neutrinos The Big Bang The Atmosphere (cosmic rays) Super novae AGN's,
More informationDecays and Scattering. Decay Rates Cross Sections Calculating Decays Scattering Lifetime of Particles
Decays and Scattering Decay Rates Cross Sections Calculating Decays Scattering Lifetime of Particles 1 Decay Rates There are THREE experimental probes of Elementary Particle Interactions - bound states
More informationElastic scattering. Elastic scattering
Elastic scattering Now we have worked out how much energy is lost when a neutron is scattered through an angle, θ We would like to know how much energy, on average, is lost per collision In order to do
More informationUNIVERSITY OF CAPE TOWN. EEE4101F / EEE4103F: Basic Nuclear Physics Problem Set 04. Due 12:00 (!) Wednesday 8 April 2015
UNIVERSITY OF CAPE TOWN EEE4101F / EEE4103F: Basic Nuclear Physics Problem Set 04 Due 12:00 (!) Wednesday 8 April 2015 1) Explain the concept of a cross-section, using language that a high-school teacher
More informationInteraction of Electrons with Matter and Comparison of the Obtained Results with Experimental Measurements
Australian Journal of Basic and Applied Sciences, 5(): 83-80, 0 ISSN 99-878 Interaction of Electrons with Matter and Comparison of the Obtained Results with Experimental Measurements Somayeh Almasi Bigdelo,
More informationTHE mono-energetic hadron beam such as heavy-ions or
Verification of the Dose Distributions with GEANT4 Simulation for Proton Therapy T.Aso, A.Kimura, S.Tanaka, H.Yoshida, N.Kanematsu, T.Sasaki, T.Akagi Abstract The GEANT4 based simulation of an irradiation
More information8/3/2016. Chia-Ho, Hua St. Jude Children's Research Hospital. Kevin Teo The Hospital of the University of Pennsylvania
Bijan Arjomandy, Ph.D. Mclaren Proton Therapy Center Mark Pankuch, Ph.D. Cadence Health Proton Center Chia-Ho, Hua St. Jude Children's Research Hospital Kevin Teo The Hospital of the University of Pennsylvania
More informationParticle Physics Homework Assignment 4
Particle Physics Homework Assignment 4 Prof. Costas Foudas March 01 Problem 1: Show the the momentum, p of a particle moving in a circular trajectory of radius, R, in a magnetic field, B, is given by:
More informationWeek 2: Chap. 2 Interaction of Radiation
Week 2: Chap. 2 Interaction of Radiation Introduction -- Goals, roll back the fog -- General Nomenclature -- Decay Equations -- Laboratory Sources Interaction of Radiation with Matter -- Charged Particles
More informationPhysics 736. Experimental Methods in Nuclear-, Particle-, and Astrophysics. Lecture 3
Physics 736 Experimental Methods in Nuclear-, Particle-, and Astrophysics Lecture 3 Karsten Heeger heeger@wisc.edu Review of Last Lecture a colleague shows you this data... what type of reaction is this?
More informationExercise 1 Atomic line spectra 1/9
Exercise 1 Atomic line spectra 1/9 The energy-level scheme for the hypothetical one-electron element Juliettium is shown in the figure on the left. The potential energy is taken to be zero for an electron
More informationInteraction theory Photons. Eirik Malinen
Interaction theory Photons Eirik Malinen Introduction Interaction theory Dosimetry Radiation source Ionizing radiation Atoms Ionizing radiation Matter - Photons - Charged particles - Neutrons Ionizing
More informationFLUKA simulations of selected topics regarding proton pencil beam scanning
FLUKA simulations of selected topics regarding proton pencil beam scanning C. Bäumer, J. Farr, J. Lambert and B. Mukherjee Westdeutsches Protonentherapiezentrum Essen, Germany T. Mertens, and B. Marchand
More informationGaseous Detectors. Bernhard Ketzer University of Bonn
Gaseous Detectors Bernhard Ketzer University of Bonn XIV ICFA School on Instrumentation in Elementary Particle Physics LA HABANA 27 November - 8 December, 2017 Plan of the Lecture 1. Introduction 2. Interactions
More informationPassage of particles through matter
Passage of particles through matter Alexander Khanov PHYS6260: Experimental Methods is HEP Oklahoma State University September 11, 2017 Delta rays During ionization, the energy is transferred to electrons
More informationMATR316, Nuclear Physics, 10 cr
MATR316, Nuclear Physics, 10 cr Fall 2017, Period II Pertti O. Tikkanen Lecture Notes of Tuesday, Nov. 28th and Thursday, Nov. 30th Department of Physics pertti.tikkanen@helsinki.fi 1 Interaction of radiation
More information3 Monte Carlo tools for proton dosimetry
3 Monte Carlo tools for proton dosimetry 3.1 Introduction As is clear from the previous chapter the problems in proton radiation transport are difficult to solve analytically, because of the statistical
More informationElectromagnetic and hadronic showers development. G. Gaudio, M. Livan The Art of Calorimetry Lecture II
Electromagnetic and hadronic showers development 1 G. Gaudio, M. Livan The Art of Calorimetry Lecture II Summary (Z dependence) Z Z 4 5 Z(Z + 1) Z Z(Z + 1) 2 A simple shower 3 Electromagnetic Showers Differences
More informationCompound Nucleus Reactions
Compound Nucleus Reactions E CM a Q CN Direct CN decays Time. Energy. Two-step reaction. CN forgets how it was formed. Decay of CN depends on statistical factors that are functions of E x, J. Low energy
More informationremsstrahlung 1 Bremsstrahlung
remsstrahlung 1 Bremsstrahlung remsstrahlung 2 Bremsstrahlung A fast moving charged particle is decelerated in the Coulomb field of atoms. A fraction of its kinetic energy is emitted in form of real photons.
More informationRadiation damage calculation in PHITS
Radiation Effects in Superconducting Magnet Materials (RESMM'12), 13 Feb. 15 Feb. 2012 Radiation damage calculation in PHITS Y. Iwamoto 1, K. Niita 2, T. Sawai 1, R.M. Ronningen 3, T. Baumann 3 1 JAEA,
More informationNeutron pulse height analysis (R405n)
Neutron pulse height analysis (R405n) Y. Satou April 6, 2011 Abstract A pulse height analysis was made for the neutron counter hodoscope used in R405n. By normalizing the pulse height distributions measured
More informationImprovements and developments of physics models in PHITS for radiotherapy and space applications
Improvements and developments of physics models in PHITS for radiotherapy and space applications L. Sihver 1-9, T. Sato 10, S. Hashimoto 10, T. Ogawa 10, K. Niita 11 1 Atominstitut, TU Wien, Austria, 2
More informationSimulation of low energy nuclear recoils using Geant4
Simulation of low energy nuclear recoils using Geant4 [ heavy ions, T < 1 kev/amu ] Ricardo Pinho ricardo@lipc.fis.uc.pt 12th Geant4 Collaboration Workshop LIP - Coimbra Hebden Bridge 14 Set 2007 Outline
More informationFLUKA manuals and Basic Input. FLUKA Beginner s Course
FLUKA manuals and Basic Input FLUKA Beginner s Course The FLUKA Manual FM.pdf in continuous development (as the program)! More a User Guide than a Reference Manual (only a short summary about physics)
More information2nd-Meeting. Ionization energy loss. Multiple Coulomb scattering (plural and single scattering, too) Tracking chambers
2nd-Meeting Ionization energy loss Multiple Coulomb scattering (plural and single scattering, too) Tracking chambers #2 -Particle Physics Experiments at High Energy Colliders John Hauptman, Kyungpook National
More informationLecture 2 & 3. Particles going through matter. Collider Detectors. PDG chapter 27 Kleinknecht chapters: PDG chapter 28 Kleinknecht chapters:
Lecture 2 & 3 Particles going through matter PDG chapter 27 Kleinknecht chapters: 1.2.1 for charged particles 1.2.2 for photons 1.2.3 bremsstrahlung for electrons Collider Detectors PDG chapter 28 Kleinknecht
More informationChapter 2: Interactions of Radiation with Matter
Chapter : Interactions of Radiation with Matter Slide set of 5 slides based on the chapter authored by J. R. Cunningham and D. R. Dance of the publication (ISBN 978-9-0-131010-1): Diagnostic Radiology
More informationQuantum Physics III (8.06) Spring 2008 Assignment 10
May 5, 2008 Quantum Physics III (8.06) Spring 2008 Assignment 10 You do not need to hand this pset in. The solutions will be provided after Friday May 9th. Your FINAL EXAM is MONDAY MAY 19, 1:30PM-4:30PM,
More informationRANGE OF ALPHAS. Advanced Laboratory, Physics 407 University of Wisconsin Madison, Wisconsin 53706
(revised 10/20/10) RANGE OF ALPHAS Advanced Laboratory, Physics 407 University of Wisconsin Madison, Wisconsin 53706 Abstract A silicon solid state detector is used to measure the energy of alphas which
More informationSimulation of Energy Loss Straggling Maria Physicist January 17, 1999
Simulation of Energy Loss Straggling Maria Physicist January 17, 1999 i This text should be read with a pinch of salt 1. ntroduction Due to the statistical nature of ionisation energy loss, large fluctuations
More informationRadiation Detection for the Beta- Delayed Alpha and Gamma Decay of 20 Na. Ellen Simmons
Radiation Detection for the Beta- Delayed Alpha and Gamma Decay of 20 Na Ellen Simmons 1 Contents Introduction Review of the Types of Radiation Charged Particle Radiation Detection Review of Semiconductor
More informationDetector Simulation. Mihaly Novak CERN PH/SFT
Detector Simulation Mihaly Novak CERN PH/SFT CERN Summer Student Program, 1 August 2017 Foreword This lecture is aimed to offer a simple and general introduction to detector simulation. Geant4 will be
More informationDedicated Arrays: MEDEA GDR studies (E γ = MeV) Highly excited CN E*~ MeV, 4 T 8 MeV
Dedicated Arrays: MEDEA GDR studies (E γ = 10-25 MeV) Highly excited CN E*~ 250-350 MeV, 4 T 8 MeV γ-ray spectrum intermediate energy region 10 MeV/A E beam 100 MeV/A - large variety of emitted particles
More informationEUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH GEANT4 SIMULATION OF ENERGY LOSSES OF SLOW HADRONS
EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH 2 September 1999 GEANT4 SIMULATION OF ENERGY LOSSES OF SLOW HADRONS V.N. Ivanchenko Budker Institute for Nuclear Physics, Novosibirsk, Russia S. Giani, M.G. Pia
More informationPhysics 736. Experimental Methods in Nuclear-, Particle-, and Astrophysics. Lecture 4
Physics 736 Experimental Methods in Nuclear-, Particle-, and Astrophysics Lecture 4 Karsten Heeger heeger@wisc.edu Homework Homework is posted on course website http://neutrino.physics.wisc.edu/teaching/phys736/
More informationEnergy loss of alpha particles - Prelab questions
Energy loss of alpha particles - Prelab questions 1. Write down the decay path from 226 Ra to 206 Pb. Show the intermediate nuclides and the nuclear reactions which cause each transformation (α/β ± decay).
More informationPHY492: Nuclear & Particle Physics. Lecture 24. Exam 2 Particle Detectors
PHY492: Nuclear & Particle Physics Lecture 24 Exam 2 Particle Detectors Exam 2 April 16, 2007 Carl Bromberg - Prof. of Physics 2 Exam 2 2. Short Answer [4 pts each] a) To describe the QCD color quantum
More informationChapter V: Cavity theories
Chapter V: Cavity theories 1 Introduction Goal of radiation dosimetry: measure of the dose absorbed inside a medium (often assimilated to water in calculations) A detector (dosimeter) never measures directly
More informationCHARGED PARTICLE INTERACTIONS
Physics Department, Yarmouk University, Irbid Jordan Dr. Nidal M. Ershaidat Phys. 649 Nuclear Instrumentation CHARGED PARTICLE INTERACTIONS Important note: This supplement is a re-edition of the published
More informationChapter Four (Interaction of Radiation with Matter)
Al-Mustansiriyah University College of Science Physics Department Fourth Grade Nuclear Physics Dr. Ali A. Ridha Chapter Four (Interaction of Radiation with Matter) Different types of radiation interact
More informationHadronic Showers. KIP Journal Club: Calorimetry and Jets 2009/10/28 A.Kaplan & A.Tadday
Hadronic Showers KIP Journal Club: Calorimetry and Jets 2009/10/28 A.Kaplan & A.Tadday Hadronic Showers em + strong interaction with absorber similarities to em-showers, but much more complex different
More informationIonizing radiation produces tracks defined by the geometry of the energy deposition events. An incident ion loses energy by Coulombic interactions
Track Structure Ionizing radiation produces tracks defined by the geometry of the energy deposition events. An incident ion loses energy by Coulombic interactions with electrons of the medium. These primary
More informationFUNDAMENTALS OF PHYSICS Vol. III - Interaction Of Nuclear Radiation With Matter- Arturo Menchaca Rocha INTERACTION OF NUCLEAR RADIATION WITH MATTER
INTERACTION OF NUCLEAR RADIATION WITH MATTER Arturo Menchaca Rocha Institute of Physics, Universidad Nacional Autonoma de Mexico, México Keywords: Nuclear Radiation, Nuclear Decay, Nuclear Residue, Ionization,
More information2. Determine the excess charge on the outer surface of the outer sphere (a distance c from the center of the system).
Use the following to answer question 1. Two point charges, A and B, lie along a line separated by a distance L. The point x is the midpoint of their separation. 1. Which combination of charges will yield
More informationPossible Interactions. Possible Interactions. X-ray Interaction (Part I) Possible Interactions. Possible Interactions. section
Possible Interactions X-ray Interaction (Part I) Three types of interaction 1. Scattering Interaction with an atom Deflected May or may not loss of energy 1 Possible Interactions Three types of interaction
More informationShielding Design Considerations for Proton Therapy Facilities
Shielding Design Considerations for Proton Therapy Facilities p p n π ± INC π 0 Nisy Elizabeth Ipe, Ph.D., C.H.P. Consultant, Shielding Design, Dosimetry & Radiation Protection San Carlos, CA, U.S.A. Email:
More informationHeavy Ion Interactions. Beginners FLUKA Course
Heavy Ion Interactions Beginners FLUKA Course 1 Overview The models DPMJET RQMD BME Input options Beam definition Transport thresholds 2 Heavy ion interaction models in FLUKA - 1 E > 5 GeV/n Dual Parton
More informationLewis 2.1, 2.2 and 2.3
Chapter 2(and 3) Cross-Sections TA Lewis 2.1, 2.2 and 2.3 Learning Objectives Understand different types of nuclear reactions Understand cross section behavior for different reactions Understand d resonance
More informationProperties of the nucleus. 9.1 Nuclear Physics. Isotopes. Stable Nuclei. Size of the nucleus. Size of the nucleus
Properties of the nucleus 9. Nuclear Physics Properties of nuclei Binding Energy Radioactive decay Natural radioactivity Consists of protons and neutrons Z = no. of protons (tomic number) N = no. of neutrons
More informationInteraction of charged particles and photons with matter
Interaction of charged particles and photons with matter Robert Miyaoka, Ph.D. Old Fisheries Center, Room 200 rmiyaoka@u.washington.edu Passage of radiation through matter depends on Type of radiation
More informationenergy loss Ionization + excitation of atomic energy levels Mean energy loss rate de /dx proportional to (electric charge) 2 of incident particle
Lecture 4 Particle physics processes - particles are small, light, energetic à processes described by quantum mechanics and relativity à processes are probabilistic, i.e., we cannot know the outcome of
More informationLECTURE 6: INTERACTION OF RADIATION WITH MATTER
LCTUR 6: INTRACTION OF RADIATION WITH MATTR All radiation is detected through its interaction with matter! INTRODUCTION: What happens when radiation passes through matter? Interlude The concept of cross-section
More informationInteraction of Electron and Photons with Matter
Interaction of Electron and Photons with Matter In addition to the references listed in the first lecture (of this part of the course) see also Calorimetry in High Energy Physics by Richard Wigmans. (Oxford
More information