Chapter 18 Scattering in Quantum Mechanics Application of 3D Green's function Masatsugu Sei Suzuki Department of Physics SUNY at Binghamton (Date: November 16, 2010) Max Born (11 December 1882 – 5 January 1970) was a German born physicist and mathematician who was instrumental in the development of quantum mechanics. He also made contributions to solid-state physics and optics and supervised the work of a number of notable physicists in the 1920s and 30s. Born won the 1954 Nobel Prize in Physics (shared with Walther Bothe). http://en.wikipedia.org/wiki/Max_Born ________________________________________________________________________ Julian Seymour Schwinger (February 12, 1918 – July 16, 1994) was an American theoretical physicist. He is best known for his work on the theory of quantum electrodynamics, in particular for developing a relativistically invariant perturbation theory, and for renormalizing QED to one loop order. Schwinger is recognized as one of the greatest physicists of the twentieth century, responsible for much of modern quantum field theory, including a variational approach, and the equations of motion for quantum fields. He developed the first electroweak model, and the first example of confinement in 1+1 dimensions. He is responsible for the theory of multiple neutrinos, Schwinger terms, and the theory of the spin 3/2 field. http://en.wikipedia.org/wiki/Julian_Schwinger ________________________________________________________________________ 18.1 Green's function in scattering theory We return to the original Schrödinger equation. 2 2 (r) V (r) (r) E (r) , 2m k or 2 2m 2m ( 2 Ek ) (r) 2 V (r) (r) . We assume that 2 E k 2 , k 2m 2m f (r) 2 V (r) (r) , 2 2 Lr k . 2 2 Lr (r) ( k ) (r) f (r) Suppose that there exists a Green's function G(r) such that 2 2 (r k )G(r,r') (r r') , with exp(ik r r') G(r,r') . 4 r r' Then (r) is formally given by 2m exp(ik r r') (r) (r) dr'G(r,r') f (r') (r) dr' V (r') (r') , 2 4 r r' where (r) is a solution of the homogeneous equation satisfying (2 k 2 )(r) 0 , 1 (r) r k exp(ik r) , (plane wave) (2 )3 / 2 with k k 18.2 Born approximation M r r - r uˆ r 0P L r |r - r r' rˆ r Here we consider the case of () (r) r r' r r'er k' ker ik rr' e eik(rr'er ) eikreik'r' for large r. 1 1 r r' r Then we have 1 2m 1 eikr () (r) exp(ik r) dr'eik'r'V (r') () (r') 3/ 2 2 (2 ) 4 r or 1 eikr () (r) [eikr f (k',k)] (2 )3 / 2 r The first term: original plane wave in propagation direction k. The second term: outgoing spherical wave with amplitude, f (k',k) , 1 2m f (k',k) (2 )3 / 2 dr'eik'r'V (r') () (r') . 2 4 The first Born approximation: 1 2m f (k',k) dr'eik'r'V (r')eikr' 2 4 when () (r) is approximated by 1 () (r) eikr . (2 )3/ 2 18.3 Differential cross section d Detector V(r) z 0 ikz f (,) e k eikz r incident spherical plane wave wave d We define the differential cross section as the number of particles per unit time d scattered into an element of solid angle d divided by the incident flux of particles. The particle flux associated with a wave function 1 1 (r) r k eikr eikz k (2 )3 / 2 (2 )3 / 2 is obtained as 1 k v N J [ * (r) (r) (r) * (r)] z z 2mi k z k k z k (2 )3 / 2 m (2 )3 / 2 unit area z k v rel relative velocity Fig. Note that = m in this figure. k volume = 1 m ikz 2 e 1 means that there is one particle per unit volume. Jz is the incident flux (number of particles) of the incident beam crossing a unit surface perpendicular to OZ per unit time. The flux associated with the scattered wave function 1 eikr f ( ) r (2 )3 / 2 r is 2 1 k f ( ) J ( * * ) r 2mi r r r r r r (2 )3 / 2 m r 2 dA r 2d 2 v f ( ) v 2 N J dA r 2d f ( ) d r (2 )3 / 2 r 2 (2 )3 / 2 (the number of particles which strike the opening of the detector per unit time) Detector r dS d 2 dS r d The differential cross section d N 2 d f ( ) d d N z or 2 f ( ) First order Born amplitude: 1 2m 1 2m f (k',k) (2 )3 k'V k dr'ei(k k')rV (r') 2 2 4 4 which is the Fourier transform of the potential with respect to q, where q k k' : scattering wave vector. For a spherically symmetric potential, f (k',k) is a function of q. q 2k sin( ) , 2 where is an angle between k’ and k (Ewald’s sphere). Fro simplicity we assume that ’ is an angle between q and r’. We can perform the angular integration over ’. 1 2m f (1) ( ) dr'eiqrV (r') 4 2 1 2m dr' d 'eiqr'cos ' 2r'2 sin 'V (r') 2 4 00 Note that 2 d 'eiqr'cos ' sin ' sin(qr') 0 qr' Then 1 2m 2 f (1) ( ) dr'2r'2 V (r') sin(qr') 2 4 0 qr' 1 2m dr'r'V (r')sin(qr') 2 q 0 The differential cross section is given by d 2 f (1) ( ) d 18.4 Yukawa potential Hideki Yukawa (23 January 1907 – 8 September 1981) was a Japanese theoretical physicist and the first Japanese Nobel laureate. http://en.wikipedia.org/wiki/Hideki_Yukawa ________________________________________________________________________ The Yukawa potential is given by V V (r) 0 er , r where V0 is independent of r. 1/ corresponds to the range of the potential. 2m 1 V 2m 1 V f (1) ( ) dr'r' 0 er' sin(qr') 0 dr'er' sin(qr') 2 2 q 0 r' q 0 Note that q dr'er' sin(qr') 2 2 0 q (1) 2m V0 1 f () 2 2 2 q Since q 2 4k 2 sin 2 ( ) 2k 2 (1 cos ) 2 so, in the first Born approximation, 2 2 d (1) 2mV0 1 f ( ) 2 2 2 2 d [2k (1 cos ) ] Note that as 0 , the Yukawa potential is reduced to the Coulomb potential, provided the ratio V0 / is fixed. V 0 ZZ'e2 2 2 2 d (1) 2 (2m) (ZZ'e ) 1 f ( ) 4 4 4 d 16k sin ( / 2) 2k 2 Using E , k 2m 2 2 d 2 1 ZZ'e 1 f ( ) 4 d 16 Ek sin ( / 2) which is the Rutherford scattering cross section that can be obtained classically. 18.5 Validity of the first-order Born approximation If the Born approximation is to be applicable, r () should not be too different from r k inside the range of potential. The distortion of the incident wave must be small. 2m eik rr' r () r k dr' V (r') r' () 2 4 r r' r () r k at the center of scattering potential at r = 0. 2m eikr' eikr' 1 dr' V (r') 2 3 / 2 3 / 2 4r' (2 ) (2 ) or 2m 1 eikr' dr' V (r')eikr' 1 2 4 r' ____________________________________________________________________ 18.6 Lipmann-Schwinger equation The Hamiltonian H is given by ˆ ˆ ˆ H H0 V where H0 is the Hamiltonian of free particle. Let be the eigenket of H0with the energy eigenvalue E, ˆ H0 E The basic Schrödinger equation is ˆ ˆ (H0 V ) E (1) ˆ ˆ ˆ Both H0 and H0 V exhibit continuous energy spectra. We look for a solution to Eq.(1) such that as V 0 , , where is the solution to the free particle Schrödinger equation with the same energy eigenvalue E. ˆ ˆ V (E H 0 ) ˆ Since (E H 0 ) = 0, this can be rewritten as ˆ ˆ ˆ V (E H0 ) (E H0 ) which leads to ˆ ˆ (E H0 )( ) V or (E Hˆ )1Vˆ 0 . The presence of is reasonable because must reduce to as Vˆ vanishes. Lipmann-Schwinger equation: () ˆ 1 ˆ () k (Ek H0 i ) V 2 2 by making Ek (= k / 2m) slightly complex number (>0, ≈0). This can be rewritten as r () r k dr' r (E Hˆ i )1 r' r' Vˆ () k 0 where 1 r k eikr , (2 )3/ 2 and ˆ H0 k' Ek ' k' , with 2 E k'2 . k ' 2m The Green's function is defined by 2 1 G() (r,r') r (E Hˆ i )1 r' eik rr' 2m k 0 4 r r' ((Proof)) 2 I r (E Hˆ i )1 r' 2m k 0 2 2 dk'dk" r k' (E k'2 i )1 k' k" k" r' 2m k 2m or 2 2 I dk'dk" r k' (E k'2 i )1 (k'k") k" r' 2m k 2m 2 2 dk' r k' (E k'2 i )1 k' r' 2m k 2m 2 dk' eik'(rr') 2m (2 )3 2 E k'2 i k 2m where 2 E k 2 . k 2m Then we have 2 dk' eik'(rr') I 2m (2 )3 2 (k 2 k'2 ) i 2m dk' eik'(rr') K () (r r') (2 )3 k'2 (k 2 i ) 0 In summary, we get 2m r () r k dr'G() (r,r') r' Vˆ () 2 or 2m r () r k dr'G() (r,r')V (r') r' () .