The Scattering and Shrinking of a Gaussian Wave Packet by Delta Function Potentials Fei
Total Page:16
File Type:pdf, Size:1020Kb
The Scattering and Shrinking of a Gaussian Wave Packet by Delta Function Potentials by Fei Sun Submitted to the Department of Physics ARCHIVF in partial fulfillment of the requirements for the degree of Bachelor of Science r U at the MASSACHUSETTS INSTITUTE OF TECHNOLOGY June 2012 @ Massachusetts Institute of Technology 2012. All rights reserved. A u th or ........................ ..................... Department of Physics May 9, 2012 C ertified by .................. ...................... ................ Professor Edmund Bertschinger Thesis Supervisor, Department of Physics Thesis Supervisor Accepted by.................. / V Professor Nergis Mavalvala Senior Thesis Coordinator, Department of Physics The Scattering and Shrinking of a Gaussian Wave Packet by Delta Function Potentials by Fei Sun Submitted to the Department of Physics on May 9, 2012, in partial fulfillment of the requirements for the degree of Bachelor of Science Abstract In this thesis, we wish to test the hypothesis that scattering by a random potential causes localization of wave functions, and that this localization is governed by the Born postulate of quantum mechanics. We begin with a simple model system: a one-dimensional Gaussian wave packet incident from the left onto a delta function potential with a single scattering center. Then we proceed to study the more compli- cated models with double and triple scattering centers. Chapter 1 briefly describes the motivations behind this thesis and the phenomenon related to this research. Chapter 2 to Chapter 4 give the detailed calculations involved in the single, double and triple scattering cases; for each case, we work out the exact expressions of wave functions, write computer programs to numerically calculate the behavior of the wave packets, and use graphs to illustrate the results of the calculations. In Chapter 5, we study the parameters that determine how much the wave function shrinks, including the initial width, the initial position and the momentum of the Gaussian wave packet, and the strength of and the spacing between the delta functions; then we examine dif- ferent combinations of the parameters in order to find a pattern to achieve maximum shrinking. Chapter 6 concludes the thesis with the essential results of this research as well as its implications and potentials. Thesis Supervisor: Professor Edmund Bertschinger Title: Thesis Supervisor, Department of Physics 3 4 0.1 Acknowledgments I would like to express my gratitude to all those who helped me to complete this thesis. I am indebted most to my supervisor Professor Edmund Bertschinger from the Department of Physics for all his inspiring suggestions, valuable guidance and encouragement in the past two years. He was patient in all the time of research for and writing of this thesis, and he was of great help in difficult times. The research was more complicated than we intially expected, and Professor Bertschinger helped me go through many bottlenecks where patience and persistence really made a difference. I also want to thank MIT's UROP office and Professor Bertschinger for giving me funding to do the research. Professor David Vogan from the Department of Mathematics and Robyn Sanderson from the Department of Physics also helped me in my research work. Professor Vogan gave me assistance on contour integrals, and Robyn Sanderson helped me with numerical simulations. I would like to thank them for all their interest and support. 5 6 Contents 0.1 Acknowledgments ... ...... ...... ...... ...... 5 1 Introduction 15 2 Single Scattering 17 2.1 Fourier Transform . .... .... .... .... ..... .... ... 17 2.2 The 01 Integral .............. 19 2.3 Table of Notifications ... ..... ..... ..... ..... .... 25 3 Double Scattering 27 3.1 <A(x, t) and ')B(X, t) .... .... - . ---- . ... - ---. 27 3.2 Analysis of IA when fi(x) < 0 .... ...... ....... ..... 31 3.3 Table of Notifications ..... ..... ..... ..... ..... 36 4 Triple Scattering 39 4.1 Triple Scattering Solution ... .... ..... .... .... .... 39 4.2 Table of Notifications ... .... .... .... .... ..... ... 45 5 Results and Analysis 47 5.1 Verification ... ...... .... ..... ..... ..... .... 47 5.2 Spacetime Color Graphs . ...... ....... ...... ..... 49 6 Conclusion 57 7 8 List of Figures 2-1 This is a schematic graph of the position of the pole k = -io in the u-plane for (1) #0 > 0 (2) 0 > 0 > -2 (3)#o < - <0, att=0. 2oO The integral is along the real k-axis. When 0 > 0, the contour of integral is always above the pole k = -io; when 0 < 0, the contour of integral is always below the pole k = -io. The position of the pole is given by equation (2.24), and the line Im(k)=0 accords with equation (2.26). ... .... .... .... .... .... .... .. 21 2-2 This is a schematic graph of shifting the contour of integral from the real k-axis to the real u-axis for (1) #0 > 0 (2) 0 > po > - (3) 0 < - ) < 0, at t = 0. The integral is originally along the real k-axis. When 0 > 0, the contour of integral is always above the pole k = -io; when 0 < 0, the contour of integral is always below the pole k = -iZ30 . But we want to replace the original contour of integral, the real k-axis, by the real u-axis to make the integral easier, so we close the contour at positive and negative infinities. ... ... ... 22 9 2-3 This is a schematic graph of the pole k = -io (1) moving up and crossing the real u-axis for #o > 0 (2) moving up and crossing the real u-axis for 0 > #3 > - fx (3) moving down and crossing the real u-axis for #o < - -2-T < 0, all for t > 0. The position of the pole is given by un(x, t) = -z(x, t) as in equation (2.21). In case (1), the pole crossing the real u-axis results in a negative residue added to the solution, because the contour going around the pole is clockwise. In case (2), the pole crossing the real u-axis results in the residue removed from the solution, because the pole is no longer enclosed by the contour. In case (3), the pole crossing the real u-axis results in a positive residue added to the solution, because the contour going around the pole is counterclockwise. .... .... .... .... .... .... ... 24 3-1 This is a schematic drawing for the scattering centers and the wave- function in three different regions for the double scattering case. 28 3-2 This is a schematic graph of the position of the pole k = -iO at t = 0 in the u-plane for integral IA when fi(x) < 0. The three cases are (1) 0 > - > 0, (2) 0 < - - , and (3) #o < 0. The integral is along the real k-axis. When 0 > 0, the contour of integral is always above the pole k = -i3o; when 0o < 0, the contour of integral is always below the pole k = -io. The position of the pole is given by equation (3.30), and the line Im(k)=0 always lies above the real u-axis when fi(x) < 0, according to equation (3.32). .. .... ... ... 32 10 3-3 This is a schematic graph of shifting the contour of integral from the real k-axis to the real u-axis for (1) #0 > - > 0, (2) 0 < #0 < - (3) /3o < 0, at t = 0. The integral is originally along the real k-axis. When #0 > 0, the contour of integral is always above the pole k = -i; when #o < 0, the contour of integral is always below the pole k = -io. But we want to replace the original contour of integral, the real k-axis, by the real u-axis to make the integral easier, so we close the contour at positive and negative infinities. ...... .... 34 3-4 This is a schematic graph of the pole k = -iO (1) moving up and crossing the real u-axis for 3o > -fi(x) > 0 (2) moving down and crossing the real u-axis for 0 < 0o < - 2'(T (3) moving down and crossing the real u-axis for 0 < 0, all for t > 0. The position of the pole is given by Up,A(X, t) = -zA(x, t) as in equation (3.27). In case (1), the pole crossing the real u-axis results in a negative residue added to the solution, because the contour going around the pole is clockwise. In case (2), the pole crossing the real u-axis results in the residue removed from the solution, because the pole is no longer enclosed by the contour. In case (3), the pole crossing the real u-axis results in a positive residue added to the solution, because the contour going around the pole is counterclockwise. ... .... ..... ..... 35 4-1 This is a schematic drawing for triple scattering case. ...... ... 40 5-1 The graphs on the left are made by our computer program, wpscatt.c. The graphs on the right tell us the difference between the graphs made by our computer program and the graph made by Mathematica with a delta function potential of infinite strength. The parameters are set to be: xi = 8 nm, oo = v2 nm, ko = 3 nm1, and t = 25 fs. ...... 54 11 5-2 This is the spacetime color graph of |@(x, t) 2 for single scattering; The parameters are set to be xi = 8 nm, oo = v nm, ko = 3 nm- 1 , and Vo = 0.15 * 10-19 J*nm. White indicates 1.0, purple 0.001, and black any value below 0.001. .... ........ ...... ....... 55 5-3 This is the spacetime color graph of |@(x, t) 2 for double scattering; zO = 1 nm, Vo = 0.15 * 10-19 J*nm, and all the other parameters have the same values as in the single scattering graph . White indicates 1.0 and purple 0.001. .... ..... .... .... ..... .... 55 5-4 This is the spacetime color graph of |@(x, t)12 for triple scattering; x1 = 2 nm, V1 = 0.15 * 10-19 J*nm, and all the other parameters have the same values as in the double scattering graph.