10. Scattering from Central Force Potential

Transcription

1 University of Rhode Island Classical Dynamics Physics Course Materials Scattering from Central Force Potential Gerhard Müller University of Rhode Island, Creative Commons License This work is licensed under a Creative Commons Attribution-Noncommercial-Share Alike 4. License. Follow this and additional works at: Abstract Part ten of course materials for Classical Dynamics (Physics 52), taught by Gerhard Müller at the University of Rhode Island. Entries listed in the table of contents, but not shown in the document, exist only in handwritten form. Documents will be updated periodically as more entries become presentable. Recommended Citation Müller, Gerhard, "1. Scattering from Central Force Potential" (215). Classical Dynamics. Paper This Course Material is brought to you for free and open access by the Physics Course Materials at DigitalCommons@URI. It has been accepted for inclusion in Classical Dynamics by an authorized administrator of DigitalCommons@URI. For more information, please contact digitalcommons@etal.uri.edu.

2 Contents of this Document [mtc1] 1. Scattering from Central Force Potential Scattering from stationary central force potential [msl2] Determination of the scattering angle [mln2] Total cross section for shower of meteorites [mex 58] Rutherford scattering formula [mex56] Scattering from hard spheres [mex55] Elastic scattering from hard ellipsoids [mex6] Scattering cross section for inverse square potential [mex59] Particle experiencing soft Coulomb kick [mex1] Scattering angle in the laboratory frame [msl3] Loss of kinetic energy in elastic collision [mex57] Elastic collision: angle between scattered particles [mex24] Elastic collision: velocities of scattered particles [mex241] Mechanical refraction [mex167] Scattering from a spherical potential well [mex168] Grazing collision between flat surfaces [mex219] Absorption cross section of power-law potential [mex242] Small-angle scattering [mln15] Small-angle scattering from-power-law potential [mex246] Classical inverse scattering [mln14] Classical inverse scattering problem I [mex243] Classical inverse scattering problem II [mex244] Classical inverse scattering problem III [mex245] Decay of particle I [mln12] Decay of particle II [mln13] Decay of particle: maximum kinetic energy [mex237]

3 Decay of particle: directions in lab frame I [mex238] Decay of particle: directions in lab frame II [mex239]

4 Determination of scattering angle [mln2] Orbital integral: ϑ = r dr 2 m Periapsis: ϑ(r min ). = ψ 2ψ + θ = π. l/mr 2 [ ]. E V (r) l 2 2mr 2 ψ ψ θ ψ = Substitute u = 1/r: dr/r 2 r 2m min [E V (r)] 1 l 2 = r 2 sdr/r 2 r min 1 V (r) E. s2 r 2 ψ = umax sdu. 1 V (1/u) s E 2 u 2 Use energy conservation to determine u max : E = 1 2 mṙ2 + l2 2mr + V (r) = l 2 + V (r 2 2mrmin 2 min ) = Es 2 u 2 max + V (1/u max ). s 2 u 2 max + V (1/u max)/e = 1 u max = u max (s, E). umax Scattering angle: θ(s, E) = π 2 sdu 1 V (1/u)/E s2 u 2. Total scattering cross section: σ T = σ(θ)dω = 2π π s sin θ ds smax dθ sin θ dθ = 2π s ds. Note: In quantum mechanics σ T can be finite even if s max is infinite.

5 [mex58] Total cross section for shower of meteorites A uniform beam of small particles with mass m 1 and velocity v is directed toward a planet of mass m 2 and radius R. Calculate the total cross section σ T for the particles to be absorbed by the planet.

6 [mex56] Rutherford scattering formula Derive the scattering cross section σ(θ) = ( κ ) 2 1 4E sin 4 (θ/2), κ = ZZ e 2 4πɛ for elastic scattering of particles with electric charge Ze and energy E from stationary atomic nuclei with charge Z e. Note that σ(θ) does not depend on whether the beam is positively or negatively charged.

7 [mex55] Scattering from hard spheres (a) Calculate the scattering cross section σ(θ) for elastic scattering from hard spheres of radius a. (b) Calculate the total scattering cross section σ T = d 2 Ω σ(θ). 1

8 [mex6] Elastic scattering from a hard ellipsoid Show that the cross section for elastic scattering from a hard ellipsoid described by the equation x 2 /a 2 + (y 2 + z 2 )/b 2 = 1 with the incident beam along the x-axis is σ(θ) = 1 4 b2 a 2 b 2 [a 2 sin 2 (θ/2) + b 2 cos 2 (θ/2)] 2. z b a ψ s x

9 [mex59] Scattering cross section for inverse square potential Show that the cross section for scattering from the stationary potential V (r) = κ/r 2 with κ > is σ(θ) = κπ2 E π θ θ 2 (2π θ) 2 sin θ.

10 [mex1] Particle experiencing soft Coulomb kick A particle with charge Q 1 and mass m 1 moves at very high velocity v 1 along a (nearly) straight line that passes at a distance b from a particle with charge Q 2 and mass m 2, which is initially at rest. The assumptions are that the two particles interact via a Coulomb central force and that the second particle does not change its position significantly during the encounter. (a) Find the direction in which the second particle will move after the encounter. (b) Find the energy E transferrred from the first to the second particle during the encounter. 1 v 1 b r θ 2

11 Scattering angle in the laboratory frame [msl3] The scattering experiment is performed in the laboratory frame. observed scattering angle: θ, observed scattering cross section: σ( θ), projectile of mass m 1 and target of mass m 2. The theoretical analysis is performed in the center-of-mass frame: problem reduced to one degree of freedom, total mass M = m 1 + m 2, reduced mass m = m 1 m 2 /(m 1 + m 2 ), calculated scattering angle: θ, calculated scattering cross section: σ(θ). Task #1: establish the relation between θ and θ.. v cm v 1 θ v v 1 θ m 1 v = (m 1 + m 2 ) v cm v cm = v 1 sin θ = v 1 sin θ, v 1 cos θ = v 1 cosθ + v cm. m 1 m 1 + m 2 v = m m 2 v. 1

12 Relative velocity after collision: v = v 2 v 1 = v 2 v 1 (frame-independent). Linear momentum in center-of-mass frame: m 1 v 1 + m 2 v 2 =. v 1 = m 2 m 1 + m 2 v, v 2 = m 1 m 1 + m 2 v m 1 v 1 = mv. tan θ = v 1 sin θ v 1 cosθ + v cm = sin θ cosθ + ρ, ρ = m m 2 v v 1 = m 1 m 2 v v. cos θ = ρ(1 cos 2 θ) + cos θ 1 ρ 2 (1 cos 2 θ). Elastic scattering: T = 1 2 m v2 = 1 2 mv2 (in center-of-mass frame) v = v ρ = m 1 /m 2. Task #2: establish the relation between σ and σ. Number of particles scattered into infinitesimal solid angle: 2πIσ(θ) sin θ dθ = 2πI σ( θ) sin θ d θ. σ( θ) = σ(θ) sin θ sin θ dθ d θ = σ(θ) d cosθ d cos θ. [ ] σ( θ) = σ(θ) 2ρ cos θ ρ2 cos(2 θ). 1 ρ2 sin 2 θ Special case: elastic scattering between particles of equal mass: m 1 = m 2 cos θ = cos(2 θ) θ = θ 2, σ( θ) = 4 cos θ 2 σ(θ). 2

13 [mex57] Loss of kinetic energy in elastic collision Consider a particle of mass m 1 and incident velocity v undergoing an elastic collision via central force with a target of mass m 2 that is initially at rest. The particle emerges with velocity v 1 from the collision as viewed in the laboratory frame. The figure shows this velocity in relation to the center-of-mass velocity v cm and the final velocity v 1 of the particle in the center-of-mass frame. Also shown are the scattering angles θ (center-of mass frame) and θ (laboratory frame). Show that the ratio of the final and initial kinetic energies in the laboratory frame is T ρ cos θ + ρ2 = T (1 + ρ) 2, ρ = m 1. m 2 v cm. v 1 θ θ v 1

14 [mex24] Elastic collision: angle between scattered particles A particle of mass m 1 and incident velocity v undergoes an elastic collision via central force with a particle of mass m 2 that is initially at rest. Given the scattering angles θ 1, θ 2 = π θ 1 in the center-of-mass frame, find the sum θ 1 + θ 2 of the scattering angles in the laboratory frame as a function of θ 1 and m 1 /m 2. Show that if m 1 = m 2 then we have θ 1 + θ 2 = π/2 for < θ 1 < π.

15 [mex241] Elastic collision: velocities of scattered particles A particle of mass m 1 and incident velocity v undergoes an elastic collision via central force with a particle of mass m 2 that is initially at rest. Show that the velocities of the scattered particles depend on the scattering angles in the laboratory frmae as follows: v 2 v = 2 m m 2 cos θ 2, m. = m 1m 2 m 1 + m 2, v 1 v = m m 2 ( cos θ 1 ± ) m 2 2 m 2 sin 2 θ1, 1 where two solutions (±) exist for m 1 > m 2 and one solution (+) for m 1 < m 2.

16 [mex167] Mechanical refraction A particle of mass m moving in the xy-plane is subject to a potential energy which assumes the constant value V 1 at y and the constant value V 2 at y <. Let us assume that V 2 < V 1. Use conversation laws to show that if the particle approaches the x-axis with speed v 1 at an angle θ 1 as shown, it will proceed with a different speed v 2 at a different angle θ 2 after crossing the line where the potential energy changes abruptly. Show in particular that the relation n sin θ 1 sin θ 2 = mv1 2 (V 1 V 2 ), between the two angles holds, where n plays the role of index of refraction. y v 1 θ 1 V 1 x V 2 θ 2 v 2

17 [mex168] Scattering from a spherical potential well Consider a spherical potential well of depth U and radius a, decribed by the potential energy V (r) = UΘ(a r). According to [mex167], the path of an incident particle with energy E encountering this potential will then be that of a ray of light refracted from a sphere with refractive index n = 1 + U/E. (a) Calculate the maximum scattering angle as a function of n. (b) Show that the scattering cross section has the form σ(θ) = a2 n 2 [n cos(θ/2) 1][n cos(θ/2)] 4 cos(θ/2) [n n cos(θ/2)] 2. α β β α s a

18 [mex219] Grazing collison between flat surfaces Consider a cube of mass m in translational motion with velocity v on a frictionless airtrack. The cube is approaching a wall at a grazing angle α with one of its sides parallel to the wall. The coefficient of kinetic friction between the cube and the wall is µ. Determine the angle β describing the direction of the velocity v the cube has after the collison. Assume that the recoil motion of the wall is negligible. v y v β m α x

19 [mex242] Absorption cross section of power-law potential A uniform beam of particles of mass m and velocity v is directed toward an attractive power-law potential V (r) = κ/r α with α > 2. Depending on the energy E and the angular momentum l the orbit of the particle leads to the center of force or it passes by at a nonvanishing minimum distance. Assume that all particles that arrive at the center of force are absorbed whereas all other particles are scattered elastically. Calculate the total cross section σ T for particle absorption as a function of α, E, κ.

20 Small-Angle Scattering [mln15] Scattering angle from transverse momentum: sinθ = p y p Impact parameter: s. Impulse and transverse momentum: p y = Transverse force: F y = V y = dv dr + r y = dv dr dt F y. θ = p y mv +... y r, r = x 2 + y 2 + z 2. Amount of transverse motion during collision assumed negligible: F y = dv dr Change in speed of particle during collision assumed negligible: Eliminate dx: x = r 2 s 2 Transverse momentum: p y = 2s Scattering angle: θ(s) = s E dt = dx p y = s + dv v v dr dx dr = v s + Scattering cross section: σ(θ) = s(θ) θ s + r r2 s 2. dv dr dv dr ds dθ Application to power-law potential: [mex246] dr r2 s 2. dx r. dr r2 s 2, E = 1 2 mv2. with s(θ) from inversion of θ(s). s r. y v s m x θ

21 [mex246] Small-angle scattering from power-law potential Consider small-angle scattering from a repulsive power-law potential V (r) = κ/r α using the relations derived in [mln15]. (a) Find the scattering cross section θ(s). (b) Find the scattering cross section σ(θ). (c) Show that the small-angle results of σ(θ) for α = 1 and α = 2 are consistent with the general results from [mex56] and [mex59], respectively.

22 Classical inverse scattering [mln14] Goal: reconstruction of potential V (r) from cross section σ(θ). Assumptions: dv/dr < (repulsive force), V () > E, V ( ) =. Consequence: θ(s 1 ) > θ(s 2 ) if s 1 < s 2. Calculate s(θ) from σ(θ): 2π π dθ sin θ σ(θ) = 2π s θ Orbital integral from [mln2] with u. = 1/r: ds s = πs 2. π θ(s) 2 = um sdu, s 2 u 2 1 V (1/u) s E 2 u 2 m + V (1/u m )/E = 1. Transformation with x. = 1/s 2, θ(x) = θ(s), and w(u). = 1 V (1/u)/E: π θ(x) 2 = um du x[w(u)]2 u 2, u2 m = x[w(u m )] 2 u m (x). Transformation: 1 2 α dx π θ(x) α x = α um dx du (x[w(u)]2 u 2)( α x ). π α α dx θ (x) um(α) du α x = π w(u). [mex243] Set α = u 2 /w 2, implying u m (α) u, take d/du, and multiply by du: π ( u ) u 2 /w 2 w dw = d w θ (x) dx u2 /w 2 x. [mex244] Integrate differentials dw and d(u/w) with consistent boundary values: ( ) 1 θ(s) w(u) = exp ds. [mex245] π w/u s2 [w(u)] 2 /u 2 The solution w(u) of this integral equation for given θ(s) thus determines V (r) from σ(θ). [Landau and Lifshitz 1976]

23 [mex243] Classical inverse scattering problem I The reconstruction of the (central force) scattering potential V (r) from the observed scattering cross section σ(θ) as outlined in [mln14] involves the transformation of the orbital integral π θ(s) 2 = um sdu, s 2 u 2 1 V (1/u) E s 2 u 2 m + V (1/u m )/E = 1, where u =. 1/r and θ(s) is the scattering angle as a function of the impact parameter, into the relation π α α dx θ (x) um(α) du α x = π w(u), u2 m = α[w(u m )] 2, for the unknown function w(u). = 1 V (1/u)/E, where x. = 1/s 2 and θ(x) = θ(s). Carry out the initial steps as indicated in [mln14]. Integrate by parts on the left and interchange the order of integrations on the right.

24 [mex244] Classical inverse scattering problem II The reconstruction of the (central force) scattering potential V (r) from the observed scattering cross section σ(θ) as outlined in [mln14] involves the conversion of the relation π α α as derived in [mex243] into the differential relation dx θ (x) um(α) du α x = π w(u), u2 m = α[w(u m )] 2, π dw ( u ) u 2 /w 2 w = d θ (x) dx w u2 /w 2 x by setting α = u 2 /w 2, taking the derivative with respect to uon both sides and multiplying back by du. Show that the boundary value u m (α) becomes the unrestricted u in the process.

25 [mex245] Classical inverse scattering problem III The final step in the reconstruction of the (central force) scattering potential V (r) from the observed scattering cross section σ(θ) as outlined in [mln14] involves the integration of the differential relation π ( u ) u 2 /w 2 w dw = d θ (x) dx w u2 /w 2 x into the integral equation w(u) = exp ( 1 π w/u ) θ(s) ds s2 [w(u)] 2 /u 2 for the quantity w(u). = 1 V (1/u)/E with u = 1/r. Carry out this step by interchanging the order of integration on the right.

26 Decay of Particle I [mln12] Particle at rest decays into two particles: Decay energy: ǫ = E () int E(1) int E(2) int Masses of decay products: m 1, m 2. Momentum conservation: p 1 + p 2 =, p 1 = p 2. = p. (change in internal energy). Energy conservation: E () int = E(1) int + p2 + E (2) int 2m + p2. 1 2m 2 Reduced mass m = m 1m 2 m 1 + m 2. ǫ = p2 2m = T 1 + T 2. Kinetic energies: T 1 = p2 2m 1 = ǫm 2 m 1 + m 2, T 2 = p2 2m 2 = ǫm 1 m 1 + m 2. Decay products move in opposite directions. All directions of p 1 equally likely. Kinetic energies T 1, T 2 determined by conservation laws alone. Particle at rest decays into three particles: Decay energy: ǫ = E () int E(1) int E(2) int E(3) int. Masses of decay products: m 1, m 2, m 3. Momentum conservation: p 1 + p 2 + p 3 =. Energy conservation: E () int = E(1) int + p2 1 + E (2) int 2m + p2 2 + E (3) int 1 2m + p m 3 Relative directions between decay products constrained but not determined by conservation laws. Kinetic energies T 1, T 2, T 3 constrained but not determined by conservation laws. Maximum kinetic energy T i limited by ǫ, m i. [mex237]

27 Decay of Particle II [mln13] Particle in motion decays into two particles. View from center-of-mass frame: Momenta: m 1 v 1 = m 2 v 2. Directions of decay products: θ 1 + θ 2 = π. v 1 θ 1 View from laboratory frame: Momentum of particle before decay: (m 1 + m 2 )v. Directions of decay products: θ1, θ 2. v θ 2 2 _ θ 1 v _ v 1 θ 1 v 1 _ θ 1 v _ v 1 v 1 θ 1 v < v 1 v > v 1 Task #1: Find the relation between θ 1 and θ 1. tan θ 1 = v 1 sin θ 1 v 1 cosθ 1 + v cosθ 1 = v v 1 sin 2 θ1 ± cos θ 1 1 v2 sin 2 v1 2 θ1. Task #2: Find the relation between θ 1 and θ 2. [mex238] Task #3: Find the range of the angle θ. = θ 1 + θ 2. [mex239]

28 [mex237] Decay of particle: maximum kinetic energy. A particle of mass M at rest decays into three particles of masses m 1, m 2, m 3 by releasing a total decay energy ɛ. Assume that mass-energy conversion is negligible (M = m 1 + m 2 + m 3 ) and that the resulting momenta of the decay products are nonrelativistic. What is the maximum kinetic energy T1 max of the emerging particle with mass m 1?

29 [mex238] Decay of particle: directions in lab frame I A particle of mass M and velocity v (in the lab frame) decays into two particles of masses m 1, m 2 by releasing a total decay energy ɛ. Assume that mass-energy conversion is negligible (M = m 1 + m 2 ) and that the resulting momenta of the decay products are nonrelativistic. Show that the angles of the directions of the decay products relative to the forward direction of the original particle satisfy the following relation: 2ɛ sin 2 ( θ 1 + θ 2 ) (m 1 + m 2 )v 2 = m 1 m 2 sin 2 θ1 + m 2 m 1 sin 2 θ2 2 sin θ 1 sin θ 2 cos( θ 1 + θ 2 ). Explain what this relation implies in the two limits v and v.

30 [mex239] Decay of particle: directions in lab frame II A particle of mass M and velocity v (in the lab frame) decays into two particles of masses m 1, m 2 by releasing a total decay energy ɛ. Assume that mass-energy conversion is negligible (M = m 1 + m 2 ) and that the resulting momenta of the decay products are nonrelativistic. (a) Calculate the angle θ =. θ 1 + θ 2 between the two emerging particles in the lab frame as a function of v, v 1, v 2 and θ 1, where θ 1, θ 2 = π θ 1 are the corresponding angles in the center-of-mass frame. (b) Determine the range of θ as a function of v under the assumption that v 1 < v 2 on a map as follows: _ θ π π/2 v 1 v 2 v