arXiv:1101.0957v1 [quant-ph] 5 Jan 2011 ne h ffc fapotential a of effect the under h anrslsadda ocuin.I re omk hspprs Appendix. paper an this in make Wronskians to the order of In properties mathematical conclusions. known draw and results res main scattering interac the the b short–range of origin Gassian general the a discuss a consider we problem we in exactly–solvable examples specialize nontrivial illustrative we As III probability. Sec. In for asymptotic center. the of coefficients the connecting equations matrix ics h e vial xcl ovbeoedmninlmodels one–dimensional solvable exactly available few the discuss h subject the oteaayi foedmninlmdl o eoac unln r tunneling resonance for models is one–dimensional of discussion analysis Present the scattering. to potential in calc interest numerical of the for quantity suitable also are They scattering. potential h td fodnr ierdffrnilequations differential linear ordinary of study the penetratio barrier describe that quantities other and scattering potential probabilities of properties analytical the on only not h ieidpnetSh¨dne qainfrapril fmass of particle a for Schr¨odinger equation time–independent The are eerto rtne ffc soeo h otsrkn pre striking most the of one is effect tunnel or penetration Barrier nSc Iw pl rnkast h croigreuto o po for Schr¨odinger equation the to Wronskians apply we II Sec. In h ups fti ae st hwta tagtowr applica straightforward a that show to is paper this of purpose The 6–16 ASnmes 36.k 36.p 03.65.Ge 03.65.Xp, 03.65.Nk, scattering numbers: PACS the and barrier Gaussian a of penetration the for exampl illustrative As wavefunction. th the of between behavior relationships totic necessary the obtain to one enable cro igreuto ral aiiaetesuyo quan of study the facilitate greatly Schr¨odinger equation eso htWosin ewe rprycoe ieryin linearly chosen properly between Wronskians that show We swl si n–iesoa oeta cteigi general in scattering potential one–dimensional in as well as I RNKASADTESCHR THE AND WRONSKIANS II. NFA(NP C aPaaCNCT,Bv.13y6 S/N, 64 y 113 Blvd. Plata–CONICET), La CCT (UNLP, INIFTA uusl4 ail eCre 6 90L lt,Argentina Plata, La 1900 16, Correo de Casilla 4, Sucursal V ( X unu cteigb Wronskians by scattering Quantum is ) − 2 ¯ h m 2 rnic .Fern´andez M. Francisco ψ .INTRODUCTION I. ′′ ( 5,22 X + ) ral aiiae h icsino h nltclpoete of properties analytical the of discussion the facilitates greatly , V ( 17,19,20 X ) ψ ( X so h aeucinlf n ih ftescattering the of right and left wavefunction the of ms lto ftetasiso rbblt n n other any and probability transmission the of ulation sw aclt h rnmsinprobability transmission the calculate we es u loo h ueia aclto ftransmission of calculation numerical the on also but n oeta scattering potential and n nne nptnilwls nSc Vw summarize we IV Sec. In wells. potential in onances 1–5 = ) oiae ya ale plcto fWronskians of application earlier an by motivated itoso unu ehnc n ottextbooks most and mechanics quantum of dictions u cteigi n ieso.They dimension. one in scattering tum DNE EQUATION ODINGER ¨ offiinsta eemn h asymp- the determine that coefficients e eactions inaddrv qain o h transmission the for equations derive and tion eoacsfraGusa well. Gaussian a for resonances nti ora hr a engetitrs in interest great been has there journal this In . m Eψ eta cteigi n ieso n derive and dimension one in scattering tential ∗ ffiinl efcnandw olc oewell some collect we self–contained ufficiently htmvsi n ieso ( dimension one in moves that ino rnkas hc r elkonin known well are which Wronskians, of tion riradaGusa el ymaso an of means By well. Gaussian a and arrier ( X 23 7,8,17–19 (1) ) eedn ouin fthe of solutions dependent . uhitrs a enfocussed been has interest Such . 9,11,12,14–16,21 −∞ . < X < ∞ ) 2
If we define the dimensionless coordinate x = X/L, where L is an appropriate length scale, then we obtain the dimensionless eigenvalue equation
1 ϕ′′(x)+ v(x)ϕ(x)= ǫϕ(x) −2 mL2 mL2E ϕ(x)= √Lψ(Lx), v(x)= V (Lx), ǫ = (2) ¯h2 ¯h2 The length unit L that renders both ǫ and v(x) dimensionless is arbitrary and we can choose it in such a way that makes the Schr¨odinger equation simpler. We will see some examples in Sec. III. In most cases of physical interest we can write the asymptotic behavior of the dimensionless wavefunction ϕ(x) as follows
A1C1(x)+ B1S1(x) for x ϕ(x) → ∞ (3) → A3C3(x)+ B3S3(x) for x → −∞ and in the intermediate region ϕ(x)= A2C2(x)+ B2S2(x) (4) The form of the functions Cj (x) and Sj (x) depends on the problem and throughout this paper we choose all of them real. On the other hand, the coefficients Aj and Bj may be complex. We discuss some examples in Sec. III; for the time being we assume that those functions satisfy the conditions (27) at a conveniently chosen point x0. By means of the Eq. (29) given in the Appendix we can easily obtain matrix expressions connecting the coefficients Ai and Bi in the three regions: A1 A2 A2 A3 = R1 , = R2 , B1 B2 B2 B3 A1 A3 = R , R = R1 R2 (5) B1 B3 · where W (C2,S1) W (S2,S1) R1 = , W (C1, C2) W (C1,S2) W (C3,S2) W (S3,S2) R2 = (6) W (C2, C3) W (C2,S3) In these equations W (f,g) denotes the Wronskian of the functions f(x) and g(x)5,22 already defined in the Appendix. If we repeat the procedure and obtain the inverse relations we appreciate that W (C ,S ) W (S ,S ) R−1 1 2 2 1 1 = − , W (C1, C2) W (C2,S1) − W (C ,S ) W (S ,S ) R−1 2 3 3 2 2 = − (7) W (C2, C3) W (C3,S2) − 3 from which it follows that 0 1 Rt J JR−1 J j = j , = (8) 1 0 − where the superscript t stands for transpose. Since the skew–symmetric matrix J satisfies Jt = J−1 = J we conclude − that − RtJ = JR 1 (9) which resembles the symplectic condition for a canonical transformation in classical mechanics24. Any matrix that satisfies Eq. (9) is said to be symplectic24. Besides, Eq. (7) tells us that the determinant of every symplectic matrix R1, R2 and R is unity. If the potential is parity–invariant (v( x)= v(x)) then C2( x)= C2(x), S2( x)= S2(x) and − − − − W (C2,S1) = W (C2,S3), W (S2,S1) = W (S3,S2) W (C1, C2) = W (C2, C3) W (C1,S2) = W (C3,S2) (10) so that we need to calculate half the number of Wronskians. It is worth noting that Wronskians containing C1 and S1 are constant for x and those with C3 and S3 are constant for x . → ∞ → −∞ III. POTENTIAL SCATTERING We assume that lim v(x) = v x→−∞ lim v(x) = v+ (11) x→∞ where v± are finite constants. If v(x) approaches those limits sufficiently fast then we know that the asymptotic behavior of the solution is ′ ik1x ′ −ik1x A1e + B1e for x ϕ(x) → −∞ , → A′ eik1x + B′ e−ik1x for x 3 3 → ∞ k1 = 2(ǫ v−), k3 = 2(ǫ v+) (12) p − p − provided that ǫ> max v+, v− . { } If we choose sin(kj x) Cj = cos(kj x), Sj (x)= , j =1, 3 (13) kj in Eq. (3) and compare it with Eq. (12) we obtain ′ 1 Bj ′ 1 Bj Aj = Aj i , Bj = Aj + i (14) 2 − kj 2 kj 4 ′ Suppose that we want to study the scattering of a particle that comes from the left (x< 0). In such a case B3 =0 and the transmission probability is given by ′ 2 k3 A3 T = | ′ |2 (15) k1 A | 1| Noting that B3 = ik3A3 we can easily rewrite the transmission probability in terms of the coefficients Aj and Bj because ′ A3 2k1A3 ′ = (16) A k1A1 iB1 1 − In the intermediate or scattering region we write ϕ(x) as in Eq. (4) where C2(x) and S2(x) are two solutions of the dimensionless Schr¨odinger equation that satisfy Eq. (27). In order to obtain them we may resort to any available numerical integration method, like, for example, Runge–Kutta25 (see also http://en.wikipedia.org/wiki/Runge–Kutta methods). Let y(x) be either C2(x) or S2(x). Typical numerical inte- gration methods yield y(x) at a set of coordinate points x0 jh where j = 0, 1,... and h is the step size. They ± simultaneously provide the derivative of the function y′(x) at the same set of points so that the numerical calculation of the Wronskians between the intermediate solutions and the asymptotic ones is straightforward. Numerical inte- gration methods like Runge–Kutta are available in many commercial and free softwares so that it is unnecessary to write a computer program for that purpose. In the present case we resorted to the fourth–order Runge–Kutta method built in the computer algebra system Derive (http://www.chartwellyorke.com/derive.html). The starting point of the ′ integration process requires y(x0) and y (x0) that are already known for the functions C2(x) and S2(x). We propagate the solution left and right till the Wronskians appearing in the matrices (6) are constant within a given error. Then we calculate the coefficients A1 and B1 in terms of A3 by means of Eq. (5) and finally the transmission coefficient from equations (15) and (16). Note that the application of the numerical integration method is straightforward because all the functions Cj (x) and Sj (x) are real. In the simplest case of a parity–invariant potential we need half the Wronskians in order to obtain the matrix R in equations (5) and (6) as discussed at the end of Sec. II. The transmission probability reads 2 k1 T = 2 2 2 2 2 2 (17) [(W (C2, C3) + k1W (C2,S3) ) (W (C3,S2) + k1W (S3,S2) )] and we only have to integrate the differential equations for C2(x) and S2(x) from x0 = 0 to the right: xj = jh, j =0, 1,...,N. As an illustrative example we choose the Gaussian barrier −αX2 V (X)= V0e , V0,α> 0 (18) −x2 If we set L = 1/√α then v(x) = v0e and the dimensionless Schr¨odinger equation depends on just one potential 2 2 2 parameter v0 = mV0/(¯h α). Since v− = v+ = 0 then k1 =2ǫ, where ǫ = mE/(¯h α) is the dimensionless energy. Throughout this paper we choose the integration step size h = 0.01 and N = 500 integration points so that the maximum coodinate value is x500 = 5. Fig. 1 shows the behaviour of the functions C2(x) and S2(x) for ǫ = 1 and the Gaussian barrier with v0 = 2. In Fig. 2 we appreciate that the Wronskians approach constants as x and | | → ∞ that x = 5 is large enough for an accurate estimation of those limits. We thus calculated the transmission coefficients shown in Fig. 3 for three values of v0. As expected there is tunneling for all ǫ> 0 and limǫ→∞ T = 1. It is clear that 5 the oscillatory behavior of T (ǫ) found by Chalk9 is due to the truncation of the Gaussian potential with the purpose of connecting the power–series solution for the intermediate region with the asymptotic plane waves10. We appreciate that the necessary truncation of the integration interval does not produce any undesirable effect on the transmission probabilities calculated in terms of Wronskians. In fact, the great advantage of the Wronskian method is that we calculate the constant asymptotic limit of each Wronskian with a given desired accuracy as shown in Fig. 2. In passing we point out that there is no trace of the questionable tunneling condition derived by Nandi26 for the Gaussia barrier. Another simple, nontrivial, and most interesting problem is the Gaussian well −αX2 V (X)= V0e , V0,α> 0 (19) − −x2 Proceeding as in the preceding example we obtain the dimensionless potential v(x)= v0e and the same expres- − sions for v0 and ǫ. In this case it is most instructive to calculate the transmission probability T in terms of v0 for fixed values of ǫ in order to reveal the scattering resonances as shown in Fig. 4. Note that the maxima of the transmission probability T = 1 that occur at some particular values of v0 are roughly independent of the energy ǫ. This well known phenomenon is better understood by means of an exactly solvable model. The scattering resonances appearing in Fig. 4 are similar to the ones exhibited by the exactly solvable well V0 V (X)= , V0,α> 0 (20) −cosh(αX)2 2 that we easily transform into the dimensionless potential v(x) = v0/ cosh(x) . Both the dimensionless potential − parameter v0 and energy ǫ have the same expressions as in the preceding examples. In this case the transmission probability reads4 2 sinh π√2ǫ 1 T = 2 , λ = 1+ √1+8v0 (21) 2 2 sinh π√2ǫ + sin(πλ) The family of curves T vs. v0 for constant ǫ resembles the one shown in Fig. 4. Note that there is full transmission T = 1 when sin(πλ) = 0; that is to say, when λ is an integer. In order to understand the origin of these resonances we pay attention to the bound states4 1 2 ǫn = (λ 1 n) , n =0, 1,... λ 1 (22) −2 − − ≤ − We appreciate that full transmission takes place when one of the excited bound–state energies (n > 0) lies exactly 2 at the rim of the well ǫn = 0 in which case λ = n + 1 (the ground–state energy ǫ0 = (λ 1) /2 is negative for all − − v0 > 0). This exactly solvable problem also proved to be useful for testing the accuracy of our programs for the calculation of the transmission probability. IV. CONCLUSIONS In this paper we propose an alternative way of approaching quantum scattering in one dimension. We think that the method based on the Wronskians between linearly independent solutions to the Schr¨odinger equation is preferable to 6 other approaches. The relatively light effort necessary to master a few mathematical properties of the Wronskians pays generously when attacking the scattering problem either analytically or numerically. In this paper we focused mainly on the latter because we are interested in nontrivial problems that are not so widely discussed in most textbooks on quantum mechanics1–4. The derivation of all the necessary scattering equations in terms of Wronskians is straightforward as well as their practical application by means of extremely simple computer programs. The calculation of the transmission probability is quite reliable if one simply checks for constant Wronskians before truncating the propagation of the solutions towards left and right in the numerical integration routine. Thus, the error due to a finite integration interval is simply the error in the lack of constant Wronskians that is easily bounded to the desired accuracy. In this way we avoid any spurious oscillations in the transmission probability as discussed in Sec.III. In addition to all that, the Wronskians have proved to be most useful for the estimation of the complex energies that describe tunnel resonances as discussed in the paper on which we based present pedagogical presentation23. One can also apply the Wronskian method to bound states and calculate their energies by simply taking into account the appropriate asymptotic behavior of the wavefunction and requiring that it be square integrable. It is clear that the approach exhibits a wide variety of useful applications and for that reason we think that it is worth teaching in advanced undergraduate or graduate courses on quantum mechanics. V. APPENDIX In order to make this paper sufficiently self–contained in this appendix we outline some well known results about the Wronskians that are useful for the study of ordinary differential equations in general22 and also for the treatment of the Schr¨odinger equation in particular5,23. To this end, we consider the ordinary second–order differential equation L(y)= y′′(x)+ Q(x)y(x)=0 (23) If y1 and y2 are two linearly independent solutions to this equation then we have d y1L(y2) y2L(y1)= W (y1,y2)=0 (24) − dx where ′ ′ W (y1,y2)= y1y y2y (25) 2 − 1 is the Wronskian (or Wronskian determinant22). Two obvious properties are: W (f,g)= W (g,f), W (f,f)=0 (26) − By linear combination of y1(x) and y2(x) we easily obtain two new solutions C(x) and S(x) satisfying ′ ′ C(x0)= S (x0)=1, C (x0)= S(x0)=0 (27) at a given point x0 so that W (C,S) = 1 for all x. If we write the general solution to Eq. (23) as y(x)= AC(x)+ BS(x) (28) 7 then A = W (y,S), B = W (C,y) (29) This equation is quite useful for deriving relationships between the coefficients of the asymptotic expansions of the wavefunction in different regions of space as shown in sections II and III. A more detailed discussion of the Wronskians is available in Powell and Crasemann’s book on quantum mechanics5. ∗ Electronic address: [email protected] 1 D. ter Haar, Problems in quantum mechanics, (Pion Limited, London, 1975). 2 A. Messiah, Quantum Mechanics, (North–Holland, Amsterdam, 1976). 3 C. Cohen–Tannoudji, B. Diu, and F. Lalo¨e, Quantum Mechanics, (John Wiley & Sons, New York, 1977). 4 S. Fl¨ugge, Practical Quantum Mechanics, (Springer-Verlag, Berlin, 1999). 5 J. L. Powell and B. Crasemann, Quantum Mechanics, (Addison-Wesley, Reading, Massachussetts, 1961). 6 B. L. Cohen, “A simple treatment of potential barrier penetration”, Am. J. Phys. 33, 97-98 (1965). 7 J. E. Draper, “Use of |Ψ|2 and flux to simplify analysis of transmission past rectangular barriers or wells”, Am. J. Phys. 47, 525-530 (1979). 8 T. A. Weber, C. L. Hammer, and V. S. Zidell, “Resonances and Analyticity of Scattering Wave Function for Square-Well- Type Potentials”, Am. J. Phys. 50, 839-845 (1982). 9 J. D. Chalk, “A study of barrier penetration in quantum mechanics”, Am. J. Phys. 56, 29-32 (1988). 10 S. Washburn, “Comment on ’A study of barrier penetration in quantum mechanics’ by J. D. Chalk”, Am. J. Phys. 56, 679 (1988). 11 C. W. Leming and A. Van Smith, “A numerical study of quantum barrier penetration in one dimension”, Am. J. Phys. 59, 441-443 (1991). 12 A. Cedillo, “Quantum mechanical tunneling through barriers: A spreadsheet approach”, Int. J.Quantum Chem. 77, 528-531 (2000). 13 G. Rawitscher, C. Merow, M. Nguyen et al., “Resonances and quantum scattering for the Morse potential as a barrier”, Am. J. Phys. 70, 935-944 (2002). 14 H Le Rouzo, “Variational R-matrix method for quantum tunneling problems”, Am. J. Phys. 71, 273-278 (2003). 15 O. M. Probst, “Tunneling through arbitrary potential barriers and the apparent barrier height”, Am. J. Phys. 70, 1110-1116 (2002). 16 A. Dutt and S. Kar, “Smooth double barriers in quantum mechanics”, Am. J. Phys. 78, 1352-1360 (2010). 17 J. H. Eberly, “Quantum scattering theory in one dimension”, Am. J. Phys. 33, 771-773 (1965). 18 P. Senn, “Threshold anomalies in one-dimensional scattering”, Am. J. Phys. 56, 916-921 (1988). 19 D. W. L. Sprung and H. Wu, “Poles, bound states, and resonances illustrated by the square well potential”, Am. J. Phys. 64, 136-144 (1996). 20 Y. Nogami and C. K. Ross, “Scattering from a nonsymmetric potential in one dimension as a coupled-channel problem”, Am. J. Phys. 64, 923-928 (1996). 21 T. M. Kalotas and A. R. Lee, “A new approach to one-dimensional scattering”, Am. J. Phys. 59, 48-52 (1991). 22 T. M. Apostol, Calculus, Second ed. (Blaisdell, Waltham, Massachusetts, 1969). 23 W. N. Whitton and J. N Connor, L., “Wronskian analysis of resonance tunneling reactions”, Mol. Phys. 26, 1511-1520 (1973). 8 24 H. Goldstein, Classical Mechanics, Second ed. (Addison-Wesley, Reading, Massachusetts, 1980). 25 W.H. Press, B. P. Flannery, S. A. Teulosky et al., Numerical recipes. The art of scientific computing, (Cambridge University Press, Cambridge, 1986). 26 S. Nandi, “The quantum Gaussian well”, Am. J. Phys. 78, 1341-1345 (2010). 9 2 1 0 -1 -2 -5 -4 -3 -2 -1 0 1 2 3 4 5 x FIG. 1: Functions C2(x) (solid line) and S2(x) (dashed line) for ǫ = 1 and the Gaussian barrier with v0 = 2 (solid line) 3 2 W(C1,S2) 1 W(S1,C2) 0 W(S1,S2) -1 -2 W(C1,C2) -3 -5 -4 -3 -2 -1 0 1 2 3 4 5 x FIG. 2: Wronskians for ǫ = 1 and the Gaussian barrier with v0 = 2 1.0 0.8 v =1 0 v0=2 v0=3 0.6 T 0.4 0.2 0.0 0 1 2 3 4 5 ε FIG. 3: Transmission probability for three Gaussian barriers 10 1.0 T ε=0.1 0.5 ε=0.05 ε=0.01 0.0 0 1 2 3 4 5 v0 FIG. 4: Transmission probability for the Gaussian well as a function of the well depth for three values of the energy