arXiv:quant-ph/0609187v3 24 Apr 2007 ipa cen hnw a bantedffato nest o h d the for intensity to diffraction the proportional obtain directly can we then screen, display ihteKrho’ a.Scin4i ueia nlss eto sasu a is 5 Section analysis, numerical is 4 ele Section the calculate law. we Kirchhoff’s 3 the section funct In with wave these approach. mechanical calculate quantum will with elec we the following, where func electronic the area, wave In the diffraction electronic law. the area, the is Kirchhoff’s incoming area where the third the area, The is double-slit equation. area the wave first is area The second areas. three into ipewv-pia prahgvsadsrpino atrwv diff wave matter ag of good description kno yielded a well has gives and As approach optics, wave-optical wave. classical simple to plane part applied by in successfully approximations, described been and be methods can wave-theoretical usually standard molecules or atoms oiin o obesi ircin fw a aclt h electronic the calculate can we if diffraction, double-slit For position. i.e., meaning, statistical fitrs hogottels er n ti xedt tm mor atom, to extend is it and matter-wave year, The last the experiments. molecules throughout diffraction many interest explain of can idea this nteve fqatmmcais h lcrnhstentr fw of nature the has electron elect the of mechanics, diffraction quantum d double-slit function of fully the view its study the we and In work, phenomena, this quantum In are approach. diffraction and interference ∗ -al [email protected] E-mail: h ae sognzda olw.I eto ecluaeteelectr the calculate we 2 section In follows. as organized is paper The h aentr fsbtmcpril lcin n etoswsp was neutrons and elections particle subatomic of nature wave The Introduction 1. ircinpten.()we h ai of ratio the when (2) patterns. diffraction h eut:()We h ltwidth slit me the quantum When with (1) diffraction so results: double-slit by the explained electronic are study experiments we these and experiments, many ircinpten 4 eas n e unu ehnc e mechanics quantum new a find also We (4) pattern. diffraction ndffato atr.()hntertoof ratio the (3)When pattern. diffraction in etdb h lcrncdul-ltdffato experimen think We diffraction double-slit patterns. electronic diffraction the electronic by tested the to affect large ewrs atrwv;Eeto obesi diffraction double-slit Electron Matter-wave; Keywords: ψ h hnmn feeto,nurn tmcadmlclrdi molecular and atomic neutron, electron, of phenomena The ASnmes 37.b 61.14.Hg 03.75.-b; numbers: PACS ( I ~r, t 2 , ,wihcnb acltdwt croigrswv qain h wa The equation. Schr¨odinger’s wave with calculated be can which ), C 60 .Isiueo ihEeg hsc,P .Bx983,Beijin 918(3), Box O. P. Physics, Energy High of Institute c. unu hoyo lcrncdul-ltdiffraction double-slit electronic of theory Quantum .Dprmn fPyis uuNra nvriy uu2731 Qufu University, Normal Qufu Physics, of Department b. and .nttt fPyis ii omlUiest,Spn 13 Siping University, Normal Jilin Physics, of a.Institute C | ψ | 70 iWang Li in-a Wu Xiang-Yao ( ψ hc eefudi experiments in found were which , ~r, t ( ~r, t ) | ) 2 ntedul-ltdffato,teeeto aefntoscnb d be can functions wave electron the diffraction, double-slit the In . | a 2 igLiu Bing , a eepanda h atcespoaiiydniya h definite the at density probability particle’s the as explained be can a ∗ a a-u Zhang Bai-Jun , si h ag f3 of range the in is a d iHiFan Xi-Hui , + a a d + a = a n 6= ( n n 1 = ( n t. 1 = b , 2 a n iQn Guo Yi-Qing and , ioJn Liu Xiao-Jing , 12345 3 l h rdcin norwr a be can work our in predictions the all , clrtoeo ugn n ichff has Kirchhoff, and Huygens of those icular , λ 2 etertclwrs nti paper, this In works. theoretical me · · · hnclapoc.W a obtain can We approach. chanical , ul-lt ic h ircinitniyis intensity diffraction the since ouble-slit, 3 o ihqatmmcaia approach. mechanical quantum with ron ∼ incnb acltdb Schr¨odinger’s by calculated be can tion rncwv ucincnb bandby obtained be can function wave tronic etta h ltthickness slit the that ffect toi aefnto ndffato area diffraction in function wave ctronic , v,adtewv sdsrbdb wave by described is wave the and ave, ,odr2 order ), nsc xeiet,teincoming the experiments, such In . asv,cmlxojcs ielarge like objects, complex massive, e emn ihmn xeiet.This experiments. many with reement · · · ato also raction srpinnesqatmmechanical quantum needs escription 50 rcinhv ensuidby studied been have ffraction aefnto sapaewv.The wave. plane a is function wave aefunction wave suae yd rgi n12 and 1923 in Broglie de by ostulated ncwv ucini h double-slit the in function wave onic ions. ircinhsbcm ag field large a become has diffraction ,teeintmsigodrin order missing isn’t there ), λ mr frslsadconclusion. and results of mmary ecnoti h obvious the obtain can we n h lsia piswt its with optics classical the wn, 00 China 6000, 009 China 100049, g n, a 5 China 65, 3 , n, c 67 4 oee,matter-wave However, . n, efunction ve ψ · · · ( ~r, t r missing are itiuigon distributing ) c a a has ψ ( ~r, t ivided has ) 2 ✻x

✟✯ a ✟✟ z d ✟ b ✟✟ ✟✟ ✲ o y

FIG. 1: The double-slit geometry, the width a, the length b and the two slit distance d.

2. Quantum theory of electron diffraction in double-slit

In an infinite plane, we consider a double-slit, its width a, length b and the two slit distance d are shown in FIG. 1. The x axis is along the slit length b and the y axis is along the slit width a. In the following, we calculate the electron wave function in the first single-slit (left slit) with Schrodinger equation, and the electron wave function of the second single-slit(right slit) can be obtained easily. At a time t, we suppose that the incoming plane wave travels along the z axis. It is

i ~ (pz Et) Ψ1(z,t)= Ae − , (1) where A is a constant. The potential in the single-slit is

V (x,y,z)=0 0 ≤ x ≤ b, 0 ≤ y ≤ a, 0 ≤ z ≤ c = ∞ otherwise, (2)

where c is the thickness of the single-slit. The time-dependent and time-independent Schrodinger equations are ∂ ~2 ∂2 ∂2 ∂2 i~ ψ(~r, t)= − ( + + )ψ(~r, t), (3) ∂t 2M ∂x2 ∂y2 ∂z2

∂2ψ(~r) ∂2ψ(~r) ∂2ψ(~r) 2ME + + + ψ(~r)=0, (4) ∂x2 ∂y2 ∂z2 ~2 where M(E) is the mass(energy) of the electron. In Eq. (4), the wave function ψ(x,y,z) satisfies boundary conditions

ψ(0,y,z)= ψ(b,y,z)=0, (5)

ψ(x, 0,z)= ψ(x,a,z)=0, (6)

The partial differential equation (4) can be solved by the method of separation of variable. By writing

ψ(x,y,z)= X(x)Y (y)Z(z), (7)

Eq. (4) becomes

1 d2X 1 d2Y 1 d2Z 2ME + + + =0, (8) X dx2 Y dy2 Z dz2 ~2 and Eq. (8) can be written into the following three equations

1 d2X + λ =0, (9) X dx2 1

1 d2Y + λ =0, (10) Y dy2 2 3 1 d2Z + λ =0, (11) Z dx2 3

where λ1, λ2 and λ3 are constants,which satisfy 2ME = λ + λ + λ . (12) ~2 1 2 3 From Eq. (5) and (6), with X(x) and Y (y) satisfying the boundary conditions X(0) = X(b)=0 Y (0) = Y (a)=0, (13) we can obtain the equations of X(x) and Y (y) d2X + λ X =0 dx2 1 X(0) = 0 X(b)=0, (14) and d2Y + λ Y =0 dy2 2 Y (0) = 0 Y (a)=0, (15) their eigenfunctions and eigenvalues are nπ X (x)= A sin x (n =1, 2, ··· ) n n b nπ λ = ( )2, (16) 1 b and mπ Y (y)= B sin y (m =1, 2, ··· ) m m a mπ λ = ( )2. (17) 2 a The solution of Eq. (11) is

2ME n2π2 m2π2 i ~2 b2 a2 z Zmn(z)= Cmne q − − , (18) and the particular solution of the wave equation (4) is

ψmn = Xn(x)Ym(y)Zmn(z) nπx mπy i 2ME n2π2 m2π2 z = A B C sin sin e ~2 − b2 − a2 n m mn b a q π2 nπx mπy i 2ME n2π2 m2 z = D sin sin e r ~2 − b2 − a2 . (19) mn b a The time-dependent particular solution of Eq. (3) is

i ~ Et ψmn(x,y,z,t)= ψmn(x,y,z)e− . (20) The general solution of Eq. (3) is

ψ2(x,y,z,t) = ψmn(x,y,z,t) Xmn ME n2π2 m2π2 nπx mπy i 2 z i Et = D sin sin e ~2 − b2 − a2 e ~ . (21) mn b a q − Xmn Eq. (21) is the electronic wave function in the first single-slit. Since the wave functions are continuous at z = 0, we have

ψ1(x,y,z,t) |z=0= ψ2(x,y,z,t) |z=0, (22) 4 from Eq. (22), we obtain nπx mπy D sin sin = A, (23) mn b a Xmn

where Dmn is a coefficient, which is 4 a b nπξ mπη Dmn = A sin sin dξdη ab Z0 Z0 b a 16A = m, n, odd mnπ2 = 0 otherwise, (24) substituting Eq. (24) into (21), we can obtain the electronic wave function in the first single-slit.

∞ 16A (2n + 1)πx (2m + 1)πy ψ (x,y,z,t) = sin sin 2 (2m + 1)(2n + 1)π2 b a m,nX=0

2ME (2n+1)2π2 (2m+1)2π2 i i z ~ Et e q ~2 − b2 − a2 e− . (25) The electron wave function in the second single-slit can be obtained by making the coordinate translation:

′ x = x ′ y = y − (a + d) ′ z = z, (26)

on substituting Eq. (26) into (25), we can obtain the electron wave function ψ3(x,y,z,t) in the second single-slit

∞ 16A (2n + 1)πx (2m + 1)π(y − (a + d)) ψ (x,y,z,t) = sin sin 3 (2m + 1)(2n + 1)π2 b a m,nX=0

2ME (2n+1)2π2 (2m+1)2 π2 i i z ~ Et e q ~2 − b2 − a2 e− . (27)

3. The wave function of electron diffraction

With the Kirchhoff’s law, we can calculate the electron wave function in the diffraction area. It can be calculated by the formula8 ikr −→ 1 e −→ ′ 1 r ψout(~r, t)= − n · [▽ ψin + (ik − ) ψin]ds, (28) 4π Zs r r r

where ψout(~r, t) is diffraction wave function on display screen, ψin(~r, t) is the wave function of slit surface (z=c) and s is the area of the aperture or slit. The Kirchhoff formula (28) is approximate, since it neglects the effect of diffraction aperture or slit on the incoming wave ψin(~r, t). However, when the diffraction aperture or slit is larger than the electron wave length the effect can be neglected. For the double-slit diffraction, Eq. (28) becomes

ikr −→ 1 e −→ ′ 1 r ψout(~r, t)= − n · [▽ ψ2 + (ik − ) ψ2]ds 4π Zs1 r r r ikr −→ 1 e −→ ′ 1 r − n · [▽ ψ3 + (ik − ) ψ3]ds. (29) 4π Zs2 r r r In Eq. (29), the first and second terms are corresponding to the diffraction wave functions of the first slit and the second slit. In the following, we firstly calculate the diffraction wave function of the first slit, it is ikr −→ 1 e −→ ′ 1 r ψout1 (~r, t)= − n · [▽ ψ2 + (ik − ) ψ2]ds, (30) 4π Zs1 r r r

2ME −→′ The diffraction area is shown in FIG. 2, where k = ~2 , s1 is the area of the first single-slit, r is the q position of a point on the surface (z=c), P is an arbitrary point in the diffraction area, and −→n is a unit vector, which is normal to the surface of the slit. ✏✏✟ P −→ ✏✏✟ 5 r ✏✟✟ ✏✏✶✏✟ −→r ′ ✏✏ ✟✯✟ ✻ ✟ −→ o ✟✟✲ R −→n

c

FIG. 2: The diffraction area of single-slit

From FIG. 2, we have −→ R ′ r = R − · −→r R −→ r ′ ≈ R − · −→r r −→ k ′ = R − 2 · −→r , (31) k then,

−→ ′ ikr ik(R r r ) e e − r ·−→ = ′ r −→r −→ R − r · r ′ ikR i−→k r e e− 2·−→ = ′ −→r −→ R − r · r ′ ikR ik r e e −→2 −→ ′ ≈ − · (|−→r | ≪ R), (32) R ~ ~r with K2 = K r . Substituting Eq. (31) and (32) into (30), one can obtain

ikR ′ −→ e ′ ′ ′ ′ −→ R ′ ′ ′ ′ ψ ~r, t e i−→k2 −→r −→n ψ x ,y ,z i k ψ x ,y ,z ds . out1 ( )= − − · · [▽ 2( ) + ( 2 − 2 ) 2( )] (33) 4πR Zs0 R In Eq. (33), the term −→′ −→′ −→′ −→ ′ ′ ′ ′ ∂ψ2( r ) ∂ψ2( r ) ∂ψ2( r ) n · ∇ ψ (x ,y ,z )| = n ′ + n ′ + n ′ 2 z=c x ∂x y ∂y z ∂z −→′ ∂ψ2( r ) = n ′ z ∂z ∞ ∞ 16A = (2m + 1)(2n + 1)π2 mX=0 nX=0 2ME (2n + 1)π (2m + 1)π i − ( )2 − ( )2 r ~2 b a n π m π 2ME (2 +1) 2 (2 +1) 2 i ~2 ( b ) ( a ) c e q − − · (2n + 1)π ′ (2m + 1)π ′ sin x sin y , (34) b a then Eq. (33) is

ikR ′ e i ∞ ∞ A ~ Et i−→k2 r 16 ψout (~r, t) = − e− e− ·−→ 1 4πR Z (2m + 1)(2n + 1)π2 s0 mX=0 nX=0 n π m π i 2ME ( (2 +1) )2 ( (2 +1) )2 c (2n + 1)π ′ (2m + 1)π ′ e ~2 − b − a · sin x sin y q b a −→ −→ 2ME (2n + 1)π (2m + 1)π −→ n · R ′ ′ [i − ( )2 − ( )2 + i−→n · k − ]dx dy . (35) r ~2 b a 2 R2 6 −→ π π −→ Assume that the angle between k2 and x axis (y axis) is 2 − α ( 2 − β), and α(β) is the angle between k2 and the surface of yz (xz), then we have

k2x = k sin α, k2y = k sin β, (36)

−→ −→ n · k2 = k cos θ, (37) −→ where θ is the angle between k2 and z axis, and the angles θ, α, β satisfy the equation π π cos2 θ + cos2( − α)+cos2( − β)=1. (38) 2 2 Substituting Eq. (36) - (38) into (35) gives

ikR n π m π i ∞ ∞ 2ME (2 +1) 2 (2 +1) 2 e Et 16A i ( b ) ( a ) c ψ (x,y,z,t) = − e− ~ e ~2 − − · out1 4πR (2m + 1)(2n + 1)π2 q mX=0 nX=0 2ME (2n + 1)π (2m + 1)π 1 [i − ( )2 − ( )2 + (ik − ) cos2 α − sin2 β] r ~2 b a R q b a ′ ′ ′ ′ ′ ′ ik sin α x (2n + 1)π ik sin β y (2m + 1)π e− · sin x dx e− · sin y dy . (39) Z0 b Z0 a Eq. (39) is the diffraction wave function of the first slit. Obviously, the diffraction wave function of the second slit is

ikR n π m π i ∞ ∞ 2ME (2 +1) 2 (2 +1) 2 e Et 16A i ( b ) ( a ) c ψ (x,y,z,t)= − e− ~ e ~2 − − · out2 4πR (2m + 1)(2n + 1)π2 q mX=0 nX=0 2ME (2n + 1)π (2m + 1)π 1 [i − ( )2 − ( )2 + (ik − ) cos2 α − sin2 β] r ~2 b a R q b a d ′ ′ ′ 2 + ′ ′ ′ ik sin α x (2n + 1)π ik sin β y (2m + 1)π e− · sin x dx e− · sin (y − (a + d))dy , (40) Z0 b Za+d a where d is the two slit distance. The total diffraction wave function for the double-slit is

ψout(x,y,z,t)= ψout1 (x,y,z,t)+ ψout2 (x,y,z,t) (41)

From the diffraction wave function ψout(x,y,z,t), we can obtain the relative diffraction intensity I on the display screen, it is

2 I ∝ |ψout(x,y,z,t)| . (42)

4. Numerical result

In this section we present our numerical calculation of relative diffraction intensity. The main input 31 8 parameters are: M =9.11×10− kg, R =1m, A = 10 , α =0.01rad, E =0.001eV and the Planck constant 34 ~ = 1.055 × 10− Js. We can obtain the relation between the diffraction angle β and relative diffraction d+a intensity I. In double-slit diffraction, we can obtain the results: (1) When the ratio of a = n(n =1, 2, 3, ···), d+a order 2n, 3n, 4n, ··· are missing in diffraction pattern. (2)When the ratio of a 6= n(n =1, 2, 3, ···), there isn’t missing order in diffraction pattern. In FIG. 3 and FIG. 4, we take a = λ, b = 1000λ and c = λ, the d+a diffraction patterns are not obvious. In FIG. 3, the ratio of a = 6, the order 6 is missing. In FIG. 4, the d+a ratio of a =6.5, there isn’t missing order. In FIG. 5, FIG. 6 and FIG. 7, we take a = λ, b = 1000λ, c = λ, the diffraction patterns are obvious, where λ = 2π~ is electronic wavelength. In FIG.5, the ratio of √2ME d+a d+a a = 3, the orders 3, 6, ··· are missing. In FIG. 6, the ratio of a = 3.4, there isn’t missing order. In d+a FIG. 7, the ratio of a = 6, the orders 6, 12, ··· are missing. In FIG. 8, FIG. 9 and FIG.10, the slit width a are corresponding to 20λ, 30λ and 50λ, their diffraction patterns are obvious. In FIG. 8, FIG. 9 and FIG. d+a 10, the ratio of a = 3, the orders 3, 6, ··· are missing. From FIG. 11 to FIG. 14, the slit thickness c is corresponding to 0, 10λ, 100λ and 1000λ. We can find that the thickness c can make a large impact on the double-slit diffraction pattern. When the slit thickness c increases the peak values of diffraction pattern increase also. 7 5. Conclusion

In conclusion, we studied the double-slit diffraction phenomenon of electron with quantum mechanical approach. We give the relation between diffraction angle and the relative diffraction intensity. We find the following results: (1) When the slit width a is in the range of 3λ ∼ 50λ we can obtain the obvious diffraction d+a patterns. (2) when the ratio of a = n(n = 1, 2, 3, ···), order 2n, 3n, 4n, ··· are missing in diffraction d+a pattern. (3)When the ratio of a 6= n(n =1, 2, 3, ···), there isn’t missing order in diffraction pattern. (4) We also find a new quantum mechanics effect that the slit thickness c has a large affect to the electronic diffraction patterns. We think all the predictions in our work can be tested by the electronic double-slit diffraction experiment.

Acknowledgement

We are very grateful for the valuable discussions with Yang Mao-Zhi.

1 O. Carnal and J. Mlynek, Phys. Rev. Lett. 66, 2689 (1991). 2 W. Sch¨ollkopf, J. P. Toennies, Science 266, 1345 (1994). 3 M. Arudt, O. Nairz, J. Voss-Andreae, C. Kwller, G. Vander Zouw, and A. Zeilinger, Nature 401, 680 (1999). 4 O. Nairz, M. Arudt and A. Zeilinger, J. Mod. Opt. 47, 2811 (2000). 5 S. Kunze, K. Dieckmann and G. Rempe, Phys. Rev. Lett. 78, 2038 (1997). 6 B. Brezger, L. Hackermuller, S. Uttenthaler, J. Petschinka, M. Arndt, A. Zeilinger, Phys. Rev. Lett. 88, 100404 (2002). 7 A.S. Sanz, F. Borondo and M.J. Bastiaans, Phys. Rev. A 71, 042103 (2005). 8 M. Schwartz, Principles of Electrodynamics, Oxford University Press, 1972. 8 I I 0.6 0.6 0.5 0.5 0.4 0.4 0.3 0.3 0.2 0.2 0.1 0.1 β(rad) β(rad)

-1.5 -1 -0.5 0.5 1 1.5 -1.5 -1 -0.5 0.5 1 1.5

FIG. 3: The relation between β and I FIG. 4: The relation between β and I with a = λ, b = 1000λ, c = λ and d = 5λ. with a = λ, b = 1000λ, c = λ and d = 5.5λ.

I 25 I 25

20 20

15 15

10 10

5 5 β(rad) β(rad)

-1 -0.5 0.5 1 -1 -0.5 0.5 1

FIG. 5: The relation between β and I FIG. 6: The relation between β and I with a = 5λ, b = 1000λ, c = λ and d = 10λ. with a = 5λ, b = 1000λ, c = λ and d = 12λ.

I I 400 20

300 15

10 200

5 100 β(rad) β(rad)

-0.75 -0.5 -0.25 0.25 0.5 0.75 -0.2 -0.1 0.1 0.2

FIG. 7: The relation between β and I with FIG. 8: The relation between β and I a = 5λ, b = 1000λ, c = λ and d = 25λ. with a = 20λ , b = 1000λ, c = λ and d = 40λ. 9 I 800 I 2500

2000 600 1500 400 1000 200 500 β(rad) β(rad)

-0.15 -0.1 -0.05 0.05 0.1 0.15 -0.1 -0.05 0.05 0.1

FIG. 9: The relation between β and I FIG. 10: The relation between β and I with a = 30λ, b = 1000λ, c = λ and d = 60λ. with a = 50λ, b = 1000λ, c = λ and d = 100λ.

I I 100 120

80 100 80 60 60 40 40 20 20 β(rad) β(rad)

-0.6 -0.4 -0.2 0.2 0.4 0.6 -0.4 -0.2 0.2 0.4

FIG. 11: The relation between β and I FIG. 12: The relation between β and I with a = 10λ, b = 1000λ, c = 0 and d = 20λ. with a = 10λ, b = 1000λ, c = 10λ and d = 20λ.

I 2000 I 17500 1750 15000 1500 12500 1250 10000 1000 7500 750 500 5000

250 β(rad) 2500 β(rad)

-0.4 -0.2 0.2 0.4 -0.4 -0.2 0.2 0.4

FIG. 13: The relation between β and I FIG. 14: The relation between β and I with a = 10λ, b = 1000λ, c = 100λ and d = 20λ. with a = 10λ, b = 1000λ, c = 1000λ and d = 20λ.