arxiv:quant-ph/ v1 22 Jan 2003

Transcription

1 Explicit solution or a Gaussian wave packet impinging on a square barrier arxiv:quant-ph/ v1 22 Jan 2003 A. L. Pérez Prieto*, S. Brouard*, and J. G. Muga *Departamento de Física Fundamental II, Universidad de La Laguna, Spain Departamento de Química-Física, Universidad del País Vasco, Apdo. 644, Bilbao, Spain Abstract. The collision o a quantum Gaussian wave packet with a square barrier is solved explicitly in terms o known unctions. The obtained ormula is suitable or perorming ast calculations or asymptotic analysis. It also provides physical insight since the description o dierent regimes and collision phenomena typically requires only some o the terms. PACS numbers: w

2 Gaussian wave packet impinging on a square barrier 2 1. Introduction A paradigmatic textbook example [1] o a quantum scattering process is the collision o a particle o mass m, initially represented by a minimum-uncertainty-product Gaussian wave packet, with a square barrier potential, V 0 i d/2 x d/2 Vx) = 1) 0 otherwise. In spite o its simplicity, this model allows to observe and study a number o interesting quantum phenomena, such as tunneling [2, 3, 4, 5, 6, 7, 8, 9], resonances [10], incidence-relection and incidence-transmission intererences [11], the Hartman eect [12, 4, 13, 14, 15, 16], and time delays or relection due to a well when the height o the potencial, V 0, is negative [17]. It is also used as a standard model to check the validity o approximate numerical methods [18, 19], or to exempliy and test dierent theories or temporal quantities such as arrival, dwell, or decay times, the asymptotic behaviour at long times [20, 21, 22], and quantum transition state theories [23]. Except or approximate analytical treatments [12], the solution o this timedependent scattering model has been always obtained by numerical methods, typically using Fast Fourier Transorms, tridiagonal systems, or by linear combination o the solutions o the stationary Schrödinger equation. Goldberg, Schey, and Schwartz were the irst to solve the model by means o a dierence equation [17]. While or many purposes these methods may be suicient, or other applications the numerical approach may be cumbersome, time consuming or even useless. Thereore, an analytical solution in terms o known unctions is o much practical interest. In general an analytical solution is not only useul as a ast computational tool; it is also valuable because o the physical insight than can be gained rom it, very oten by means o approximations that extract the dominant contributions in dierent limits, parameter ranges and regimes. In this paper we present an explicit expression in the momentum representation o the time dependence o a Gaussian packet incident on a square barrier in terms o known unctions. There are a ew other time dependent solvable scattering models: a Lorentzian state impinging on a delta unction potential [24, 25, 26] or on a separable exponential potential [26, 27]; and a cuto plane wave impinging on a square barrier [28, 29], or on a step barrier [30]. While all have been useul and illustrative o several scattering properties, the present model is the only one that combines simultaneously a local potential, resonances, and a physical Hilbert-space initial state with inite moments. 2. Explicit solution The initial wave unction is assumed to be a minimum-uncertainty-product Gaussian wave packet located in the let hal-line and with negligible overlap with the potential.

3 Gaussian wave packet impinging on a square barrier 3 In momentum representation it is given by ) 1/4 p 2δx ψt = 0) = π h 2 exp [ δ x p p c ) 2 / h 2 ip x c / h ], 2) where x c and p c are the initial mean position and momentum respectively, and δ x is the variance o the state in coordinate representation. Its time evolution can then be written in terms o eigenstates o H, φ p [31, 32], as ψp,t) = + p φ p e ie t/ h p ψt = 0) dp, 3) with E p 2 /2m. These states are given explicitly in coordinate representation by x φ p = 1 Ie ip x/ h +Re ip x/ h x < d/2 Ce ip x/ h +De ip x/ h d/2 x d/2 4) h 1/2 Te ip x/ h x > d/2, where p p 2 2mV 0. The value o the coeicient I will be taken as 1. Note that these delta-normalized) eigenstates are associated with an incident plane wave o momentum p or p > 0 and with an outgoing plane wave with momentum p or p < 0. The coeicients R, C, D and T are determined by continuity o the wave unction and its irst derivative, where Rp ) = i ) p p p p sinp d/ h)e ip d/ h 2Ωp ) ) 1+ p Cp p e ip +p )d/2 h ) = 2Ωp ) ) 1 p Dp p e ip p )d/2 h ) = 2Ωp ) T p ) = e ip d/ h Ωp ), 5) Ωp ) cosp d/ h) i ) p + p sinp d/ h). 6) 2 p p Fourier transormation o x φ p in 4) gives ive terms proportional to the coeicients I, R, T, C, D respectively generically denoted by A hereater). Substituting these terms into 3) with the initial state in 2), we obtain ψp,t) = iτ hh 1/2 [g I p )+g R p )+g C p )+g D p )+g T p )]e φp ) dp,7) where τ 2π h) 1/2 2δ x π h 2 ) 1/4. The term in the exponential φp ) is given by φp ) = ip 2 t 2m h δ xp p c ) 2 h 2 ip x c h ip d 2 h ; 8)

4 Gaussian wave packet impinging on a square barrier 4 and g I p ) e ipd/2 h p p +i h0 + g R p ) Rp )e ip d/ h e ipd/2 h p+p +i h0 + g C p ) 2iCp )e ip d/2 h p p g D p ) 2iDp )e ip d/2 h p+p sin[p p )d/2 h] sin[p+p )d/2 h] g T p ) Tp )e ip d/ h e ipd/2 h p p i h0 +. 9) Note the three explicit structural [33]) poles at p I p + i h0 +, p T p i h0 +, and p R p i h0+, in addition to the poles o the unctions C, D, R and T, p j j = 1,..., ), which are zeros o Ωp ). All the poles lie in the lower hal complex plane except or the incidence pole, p I. The integral in Eq. 7) can be solved by completing the square in Eq.8) and introducing the variable u as with u p z, 10) ) 1/2 δx h 2 +i t. 11) 2m h u is zero at the saddle point z, deined by z m[ 4mp c δ 2 x x c +d/2) h 2 t ] i2m h[mδ x x c +d/2)+p c δ x t] 4m 2 δ 2 x +t 2 h 2, 12) andbecomes real along thesteepest descent path, astraight linewith slope t h/2mδ x ). We now deorm the contour and integrate along the steepest descent path, namely along the real u axis, where ψp,t) = iτ hh 1/2 e δxp2 c/ h 2 )+η 2 [g I u)+g R u)+g C u)+g D u)+g T u)]e u2 du, 13) Γ u η 2pc δ x h 2 i x )[ )] 1/2 c +d/2) δx 4 h h 2 +i t, 14) 2m h g A u) g A p = u+z), and Γ u goes rom to + including a circle around the poles that have been crossed by the contour deormation. Since the g I, g R, g T and g CD g C + g D are meromorphic unctions with simple poles, it is useul to extract explicitly the singularities and leave the remainder as an entire unction, h, g I u) = R I u u I

5 Gaussian wave packet impinging on a square barrier 5 g R u) = R R u u R + g CD u) = j=1 g T u) = R T u u T + j=1 R Rj u u j +h R u) R Cj +R Dj +h CD u) u u j j=1 R Tj u u j +h T u). 15) Again the poles in the u complex plane may be separated into structural [33], u I = 1 p+i0 + z) u R = 1 p i0 + z) u T = 1 p i0 + z), 16) and resonance poles see the Appendix), u j = 1 p j z) j = 1,...,. 17) The residues R are given by where R I = 1 e ipd/2 h R R = 1 Rp R)e ip R d/ h e ipd/2 h R Rj = Rp j )eip j d/ h e ipd/2 h F j p+p j +i h0 + ) R Cj = 2iCp j )eip j d/2 h sin[p p j )d/2 h]f j p p j ) R Dj = 2iDp j)e ip j d/2 h sin[p+p j)d/2 h]f j p+p j ) R T = 1 Tp T )eip T d/ h e ipd/2 h R Tj = Tp j )eip j d/ h e ipd/2 h F j p p j i h0 + ), 18) F j Ωp ) dωp )/dp 19) p =p j. Taking into account the integral expression o the w-unction, wz) = e z2 erc iz) [34], wz) = 1 e iπ Γ u2 du, 20) u z where Γ goes rom to passing below the pole, and the relation w z) = 2e z2 wz), the wave unction may inally be written as ψp,t) = iτ hh 1/2 e δxp2 c / h2 )+η 2{ iπr I wu I ) iπr R w u R ) iπr T w u T ) iπ [R Rj +R Cj +R Dj +R Tj ]w u j ) j=1

6 Gaussian wave packet impinging on a square barrier ψp,t) 2 a.u) Momentum a.u.) Figure 1. ψp,t) 2 as a unction o momentum at t = 5 a.u., when the collision has been completed. The solid line corresponds to the exact solution, whereas the dots correspond to the relection and transmission structural) terms and three resonance poletermsrom21). Atomicunitsareusedinallthecalculations,andm = 1. p c = 20, d = 2.5, V 0 = 188, δ x = 100 and x c is located 50 atomic units to the let o the center o the barrier potential. + [h R u)+h C u)+h D u)+h T u)]e u2 du }, 21) which is the main result o this paper. The remaining Gaussian integrals are quite generally only a minor correction and may be evaluated with a ew terms o the series hu) = h n) 0) un n!, 22) which gives n=0 hu)e u2 du = π [ hu = 0)+ n= n 1) 2 n 2n)! h 2n) u = 0) Note the orm o the solution in 21). There are structural terms associated with incidence, transmission and relection I, T, and R terms or short), resonance terms, and the h-corrections o 23). ].23) 3. Examples In many cases just a ew terms in 21) provide an accurate approximation to the exact result. I only the asymptotic behaviour, ater the collision has been completed, is o interest, the structural terms R and T, associated with relection and transmission, plus a ew resonance pole contributions again or relection and transmission) are quite suicient. The number o resonance poles that have to be added depends on the energy o the collision. In igure 1, or example, we show a case where only three resonance poles have been used in the R and T terms, and no correction term rom Eq. 23) has been included. There are also cases where only two terms are required during the collision: the collision with a very opaque barrier R and I terms), see Fig. 2, and collisions or

7 Gaussian wave packet impinging on a square barrier 7 2 ψp,t) 2 a.u.) Momentum a.u.) Figure 2. ψp,t) 2 as a unction o momentum, at t = 2.5 a.u., during the collision process. The solid line corresponds to the exact solution whereas the dots correspond to incidence and relection structural terms o 21) only, and no resonance poles taken into account. p c = 20, d = 3, V 0 = 400, δ x = 100, and the initial wave packet is centered 50 atomic units to the let o the barrier center ψp,t) 2 a.u.) Momentum a.u.) Figure 3. ψp,t) 2 as a unction o momentum, at t = 2.7 a.u., during the collision process. The solid line corresponds to the exact solution, whereas the dots correspond to taking rom 21) incidence, relection, transmission, and barrier C and D) terms, plus three resonance pole terms. When one correction term 23) is added, the approximate result is indistinguishable rom the exact one. Same barrier and wave packet parameters as in Fig. 1 very broad an energetic wave packets I and T terms). This later case has been used recently to explain a transient, classically orbidden enhancement o the high momentum components [35, 36] as an incidence-transmission intererence eect [11]. In Figure 3, the case o a very general typical) situation is depicted. Incidence, relection, transmission, and barrier C and D) terms rom 21), plus three resonance pole terms provide a airly good approximation or the wave unction at an intermediate time during the collision, although the approximate solution clearly overestimates the exact value at the relection and transmission peaks, and or values o momentum

8 Gaussian wave packet impinging on a square barrier 8 between 20.3 and 21 a.u.. The inclusion o the irst correction term in 23) leads to a curve which is indistinguishable in the scale o the igure rom the exact result. Let us also note that 21) is suitable or asymptotic analysis e.g. short and large times) by means o the asymptotic series o the w-unction, see e.g. [28, 27]. The w- unctions may be considered the elementary propagators o the Schrödinger transient modes [37] and play a prominent role in the explicit solution given in 21). Acknowledgments AP acknowledges support by CajaCanarias; SB acknowledges support by the Ministerio de Ciencia y Tecnología o Spain Project BFM ); JGM acknowledges support by the Ministerio de Ciencia y Tecnología BFM C03-03), UPV-EHU /2001), and the Basque Government PI ). Appendix A. Calculation o the resonance poles The resonance poles are obtained by solving the transcendental equation Ωp ) = 0. A systematic way to do it is based on the dierential equation [38] dp = idv 0 m p dv 0 V idp 0 m 2 ). A.1) As boundary condition or the integration o A.1) we take p = 0 at V 0 = n2 π 2 ; n = 0,±1,±2,... A.2) 2d2 The integration o A.1) has to be carried out until the value o V 0 corresponds to the physical barrier we want to examine. The integration however cannot be perormed along the real axis because o the actor V 0 in the denominator o A.1). To avoid the singularity at V 0 = 0, an imaginary term is added to the potential, analytically continuing the dierential equation. A parabolic path has been used across the complex V 0 plane rom the initial values o V 0 in A.2) to the inal real value o V 0. Reerences [1] Schi L I 1981 Quantum Mechanics, 3rd edition New York: McGraw-Hill) [2] Huang Z H et al IEEE Transactions on Elec. Dev [3] Leavens C R 1990 Solid State Comm [4] Brouard S, Sala R and Muga J G 1994 Phys. Rev. A [5] Krenzlin H M, Budczies, and Kehr K W 1996 Phys. Rev. A [6] Xavier A L and Aguiar M A M 1997 Phys. Rev. Lett [7] Aquino V M, Aguilera-Navarro V C, Goto M and Iwamoto H 1998 Phys. Rev. A [8] Abolhasani M and Golshani M 2000 Phys. Rev. A [9] Petrillo V and Realdi L 2000 Optics Commun [10] Bramhall M H and Casper B M 1970 Am. J. Phys [11] Pérez A, Brouard S and Muga J G 2001 Phys. Rev. A [12] Hartman T E 1962 J. Appl. Phys

9 Gaussian wave packet impinging on a square barrier 9 [13] Delgado V and Muga J G 1996 Annals o Physics NY) [14] Palao J P, Muga J G, Brouard S and Jadczyk A 1997 Phys. Lett. A [15] Ruschhaupt A 1998 Phys. Lett. A [16] Ruschhaupt A 2000 in Decoherence: Theoretical, Experimental and Conceptual Problems ed Ph Blanchard et al Berlin: Springer) [17] Goldberg A, Schey H M and Schwartz J L 1967 Am. J. Phys [18] Turner R E and Snider R F 1987 J. Chem. Phys [19] Muga J G 1991 J. Phys. A [20] Leon J, Julve J, Pitanga P and Urríes F J 2000 Phys. Rev. A [21] Baute A D, Egusquiza I L and Muga J G 2001 Phys. Rev. A [22] Muga J G, Sala Mayato R and Egusquiza I L ed) 2002 Time in Quantum Mechanics Berlin: Springer) [23] Muga J G, Delgado V, Sala R and Snider R F 1996 J. Chem. Phys [24] Ebereldt W and Kleber M 1988 Am. J. Phys [25] Muga J G, Brouard S and Snider R F 1992 Phys. Rev. A [26] Muga J G, Delgado V and Snider R F 1995 Phys. Rev. B [27] Muga J G and Palao J P 1998 J. Phys. A [28] Brouard S and Muga J G 1996 Phys. Rev. A [29] García-Calderón G and Villavicencio J 2001 Phys. Rev. A [30] Delgado F, Cruz H and Muga J G 2002 J. Phys. A 35 1? Preprint quant-ph/ [31] Hammer C L, Weber T A and Zidell V S 1997 Am. J. Phys [32] Baute A, Egusquiza I L and Muga J G 2002 Int. J. Theor. Phys. Group Theor. and Nonl. Optics. 8 1 Preprint quant-ph/ [33] Wei G W and Snider R F 1995 Phys. Rev. A [34] Abramowitz M and Stegun I Aed) 1972 Handbook o Mathematical Functions New York: Dover) [35] Brouard S and Muga J G 1998 Phys. Rev. Lett. 81, 2621 [36] Brouard S and Muga J G 1998 Ann. Phys. Leipzig) [37] Nussenzveig H M 1992 in Symmetries in Physics ed A Franck and K B Wol Berlin: Springer- Verlag) p 293 [38] Sokolovski D, Brouard S and Connor J N L 1994 Phys. Rev. A