arxiv:quant-ph/ v1 21 Feb 2001

Transcription

1 Explicit summation of the constituent WKB series and new approximate wavefunctions arxiv:quant-ph/ v1 21 Feb Introduction Vladimir V Kudryashov and Yulian V Vanne Institute of Physics, National Academy of Sciences of Belarus, F Skaryna avenue 70, Minsk, Belarus Department of Theoretical Physics, Belarusian State University, F Skaryna avenue 4, Minsk, Belarus kudryash@dragon.bas-net.by and vanne@dragon.bas-net.by Abstract. The independent solutions of the one-dimensional Schrödinger equation are approximated by means of the explicit summation of the leading constituent WKB series. The continuous matching of the particular solutions gives the uniformly valid analytical approximation to the wavefunctions. The detailed numerical verification of the proposed approximation is performed for the exactly solvable problems with different kinds of potentials. PACS numbers: Lt, Sq The perturbation theory, the variational method and the WKB approximation are very extensively used in quantum mechanics. If we deal with the perturbation theory or with the variational method then the similar questions arise. How to find the unperturbed Hamiltonian or how to find the trial function for an arbitrary given potential? The universal answers are absent. In this sense both mentioned methods are not complete. At the same time the WKB approximation is directly determined by a given potential. But the conventional WKB approximation is not valid near turning points. An old problem in quasiclassical analysis is the development of global uniform approximations to the wavefunctions. In previous works [1, 2] the essential improvement of the WKB approach was done with the help of the new quasi-uniform approximation. In the present paper we shall construct the second order quasi-uniform approximation to the wavefunctions. The quality of the approximate wavefunctions will be verified by means of the comparison with the exact solutions for different kinds of potentials. We shall consider the linear one-dimensional Schrödinger equation d 2 Ψ(q, h) = Q(q) dq 2 2 Ψ(q, h) (1) h where Q(q) = 2m (V (q) E) for an arbitrary potential V (q). A logarithmic derivative d lnψ(q, h) Y (q, h) = dq

2 of a wavefunction Ψ(q, h) satisfies the nonlinear Riccati equation dy (q, h) dq + (Y (q, h)) 2 = Q(q) h 2. (2) The WKB approach deals just with functions Y (q, h). In this approach two independent solutions Y ± (q, h) of the Riccati equation are represented by their asymptotic expansions ( ) Y as ± (q, h) = h 1 ±Q 1/2 + h n Y n ± (q) (3) n=1 in powers of Plank s constant h. The usual WKB approximation contains a finite number of leading terms Y n ± ± (q) from the complete expansion Y as (q, h). As it is well known, the WKB series is divergent. The direct summation of a divergent series does not exist. By summing one means finding a function to which this series is the asymptotic expansion [3]. In recent years many studies have been devoted to extracting some useful information about the exact eigenfunctions from the divergent WKB series (see for instance [4] and references therein). There are several examples of the summation for eigenvalues [5, 6, 7]. Unlike entirely exact but very complicated methods for some classes of potentials (see for example [8]) our new way of using the WKB series gives the approximate but very simple and universal method of solving the Schrödinger equation The explicit summation of the constituent WKB series In [2] the asymptotic WKB series was reconstructed as an infinite sum Y as(q, ± h) = ± h 1 Q 1/2 + Z as,j(q, ± h) (4) j=1 of new constituent (partial) asymptotic series Z as,j ± (q, h) = ( h2/3 ) j 2 ( Q n=j h 2/3 ) j 3n+1 2 A ± n,j (Q (q), Q (q),...q (j) (q)) (5) in powers of the ratio Q/ h 2/3 where Q (q) = dq(q)/dq, Q (q) = d 2 Q(q)/dq 2 and Q (j) (q) = d j Q(q)/dq j. We have simple expressions A ± n,1 = (Q ) n B ± n,1, A ± n,2 = Q (Q ) n 2 B ± n,2 for two leading sequences of coefficients A ± n,j where B n,1 ± and B n,2 ± are numbers. The complete series Y as ± (q, h) was approximated by a finite number of leading constituent series Z as,j(q, ± h) in contrast with a finite number of leading terms Y n ± (q) in the conventional WKB approach. If we can find functions Z j ± (q, h) which are represented by asymptotic expansions Z as,j(q, ± h) then we obtain new approximations to solutions of the Riccati equation. We shall consider the second order approximation. Introducing a dimensionless variable a(q, h) = 1 Q(q) (6) h 2/3 Q (q) 2/3

3 we are able to rewrite leading constituent expansions in the form ± h 1 Q 1/2 + Z as,1 ± (q, h) + 1 Q Q Z± as,2 (q, h) = h 2/3 Q 2/3y± as,1 (a) + Q y± as,2 (a). (7) Here we separate the asymptotic series in a y as,1 ± (a) = ±a1/2 + B n,1 ± a1 3n+1 2, (8) y ± as,2(a) = n=2 n=1 B ± n,2a 2 3n+1 2. (9) The major terms of these series may be restored from [3, 5] y ± m,1 (a) = ±a1/2 1 4 a 1, (10) y ± m,2 (a) = ±1 8 a 3/ a 3. (11) Our aim is to sum constituent series (8) (9). In other words we must find functions y ± j (a) which are represented by these expansions. In order to perform the identification we substitute the approximate function Y ap ± 1 Q Q (q, h) = h 2/3 Q 2/3y± 1 (a) + Q y± 2 (a) (12) into the Riccati equation (2). As a result we get following equations dy 1 da + y2 1 = a, (13) dy 2 da + 2y 1y 2 = 1 ( 2a dy ) 1 3 da y 1 (14) for functions y j (a). The equation (13) is the Riccati equation for logarithmic derivatives of the well studied Airy functions [9]. We select particular solutions by means of the known asymptotics (10). In classically allowed region where Q(q) < 0 (a < 0) we derive explicit expressions y 1 ± (a) = d ln (Bi(a) iai(a)) (15) da and in classically forbidden region where Q(q) > 0 (a > 0) we get other solutions ỹ1 (a) = d da ln Ai(a), ỹ+ 1 (a) = d ln Bi(a). da (16) Finally, we can obtain solutions of the linear equation (14) with asymptotics (11) in the closed form y 2 ± (a) = 1 [ 8a 2 (y 1 ± (a)) 2 4ay 1 ± (a) + 8a 3 3 ], 30 (17) ỹ 2 ± (a) = 1 [ 8a 2 (ỹ 1 ± (a)) 2 4aỹ 1 ± (a) + 8a 3 3 ]. 30 (18) 3

4 Although the above functions (15) (18) have asymptotic expansions (8) (9) when a is large it should be stressed that the obtained functions have other expansions when a is small. Replacing y ± j by ỹ ± j in the expression (12) we get the second pair Ỹ ap ± 1 Q (q, h) = h 2/3 Q 2/3ỹ± 1 of approximate solutions. Q (a) + Q ỹ± 2 (a) (19) It is not surprising that asymptotics of our approximation coincides with WKB asymptotics far away from turning points. At the same time our approximation reproduces the known [3] satisfactory approximation near turning points. Naturally, our approximation gives the exact result for the linear potential of the uniform field. Note that this potential presents an example of the explicit summation of the WKB series for a logarithmic derivative of a wavefunction Approximate wavefunctions for two-turning-point problem So we have the uniformly valid approximation to solutions of the Riccati equation. Now we can construct approximate wavefunctions. We shall consider the problem with two real turning points q and q +. The potential has the minimum at point q m. The first derivative and the second derivative of the smooth potential are continuous at point q m. Two pairs of independent solutions of the Schrödinger equation are approximated by functions ( q ) Ψ ± ap (q) = exp Y ± ap (q )dq and Ψ ± ap(q) = exp ( q Ỹ ± ap(q )dq ). In classically forbidden regions we must retain only decreasing solutions Ψ (q). By matching particular solutions at points q, q + and q m we obtain the continuous approximate wavefunction Ψ ap (q) which is represented by following formulas Ψ 1 (q) = C cos π ( 3 exp q ) Ỹap (q )dq (20) when q < q, ( q Y ap(q + ) + Yap(q ) ( ) q Ψ 2 (q) = C exp dq cos q 2 q q Q Q Y + ap(q ) Y 2i ap(q ) dq π ) 3 when q < q < q + and Ψ 3 (q) = C( 1) n cos π ( ) ( q 3 exp q+ Ỹap (q )dq Y ap + exp (q ) + Y ) ap (q ) dq q + q 2 when q > q +. Here we have the new quantization condition q+ q Q Q Y + ap (q, E) Yap (q, E) 2i (21) (22) dq = π(n + 2 ), n = 0, 1, 2... (23) 3

5 which determines the spectral value E sp (n) of energy implicitly. We denote the wavefunctions with E = E sp (n) as Ψ ap (q, n). Then we may choose value of an arbitrary constant C in order to ensure the usual normalization Ψ ap (n) Ψ ap (n) = 1 where Ψ ap (n) is the vector in the Hilbert space which corresponds to the function Ψ ap (q, n). The proposed approximation is an alternative to the well known Langer approximation [10] which employs a h-expansion different from the WKB series. So the approximate eigenfunctions are defined completely. But the question arises about the optimal approximate eigenvalues because the value E sp (n) is not unique choice. Since explicit expressions for wavefunctions are already obtained we are able to calculate expectation values of the Hamiltonian Ē(n) = Ψ ap (n) Ĥ Ψ ap(n) (24) Ĥ = h2 2m dq + V (q). 2 In accordance with the eigenvalue problem d 2 Ĥ Ψ E Ψ = 0 we construct the discrepancy vector D(e, n) = Ĥ Ψ ap(n) e Ψ ap (n) with an arbitrary parameter e while Ĥ and Ψ ap(n) are given. It is natural to require that the discrepancy vector should not contain a component proportional to the approximate eigenvector. In other words we consider the orthogonality condition Ψ ap (n) D(e, n) = Ē(n) e = 0 as a criterion for selection of the optimal approximate eigenvalue. As a result we get just Ē(n) while E sp(n) does not fulfil the above requirement. It should be noted also that the scalar product D(e, n) D(e, n) is minimized at e = Ē(n) Verification of the quasi-uniform approximation Now we must verify our approximation numerically for the exactly solvable problems. We shall compare the normalized approximate wavefunctions Ψ ap (q, n) with the normalized exact wavefunctions Ψ ex (q, n). In order to estimate the closeness of two functions f 1 (q) and f 2 (q) we consider two corresponding vectors f 1 and f 2 in the Hilbert space. Then we construct a deviation vector f = f 1 f 2 and a scalar product f f = f 1 f 1 + f 2 f 2 f 1 f 2 f 2 f 1. Now we can define the relative deviation δf = 1 f 1 f 2 + f 2 f 1 f 1 f 1 + f 2 f 2 (25)

6 as the numerical estimation of the closeness of two functions. Note that δf = 0 when f 1 (q) = f 2 (q). Thus, we get the following estimation δψ(n) = 1 Ψ ex (n) Ψ ap (n) (26) in the case of normalized real functions Ψ ex (q, n) and Ψ ap (q, n). The same numerical comparison may be performed δψ (n) = 1 2 Ψ ex (n) Ψ ap (n) Ψ ex (n) Ψ ex (n) + Ψ ap (n) Ψ ap (n) (27) for first derivatives Ψ (q) = dψ(q)/dq. Naturally, we can define analogous estimations for higher derivatives. Moreover, we estimate the closeness of functions ĤΨ ap(q, n) and ĤΨ ex(q, n) with the help of the following quantity δĥψ(n) = Ψ ap(n) Ĥ2 Ψ ap (n) + E ex (n) 2 (1 2 Ψ ex (n) Ψ ap (n) ) Ψ ap (n) Ĥ2 Ψ ap (n) + E ex (n) 2 (28) where E ex (n) is the exact energy value. Besides, we compare two functions ĤΨ ap(q, n) and Ē(n)Ψ ap(q, n). This comparison may be performed when we do not know exact solutions. As a result we get the relative discrepancy d(n) = Ψ ap(n) Ĥ2 Ψ ap (n) Ē(n)2 Ψ ap (n) Ĥ2 Ψ ap (n) + Ē(n)2 (29) which is directly connected with the initial Schrödinger equation under consideration. Finally, we characterize our approximation by the usual relative energy error E(n) = Ē(n) E ex (n) 1 (30) and by the nondiagonal matrix elements I ij Ψ ap (i) Ψ ap (j) i < j. (31) The verification is performed for three potentials with different asymptotics. They are the harmonic oscillator potential the Morse potential V (q) = kq 2, V (q) = γ2 h 2 α 2 ( e 2αq 2e αq) 2m and the modified Pöschl-Teller potential V (q) = λ(λ 1) h2 α 2 2m cosh 2 (αq). The exact solutions of the Schrödinger equation for these potentials may be found in [11]. 6

7 7 n δψ(n) δψ (n) Table 1. Numerical estimation of the quasi-uniform approximation δĥψ(n) d(n) E(n) Matrix elements The harmonic oscillator potential I 02 = I 13 = The Morse potential (γ = 4.5) I 01 = I 02 = I 03 = I 12 = I 13 = The modified Pöschl-Teller potential (γ = 5) I 23 = I 02 = I 13 = Only four levels exist for this value of parameter. Only nonzero matrix elements are presented. Table 1 show that the quasi-uniform approximation gives fairly accurate results for all considered potentials and for all considered quantities. Hence, the performed reconstruction of the WKB series and subsequent explicit summation of the leading constituent (partial) series yield the satisfactory (qualitative and quantitative) description of wavefunctions. References [1] Kudryashov V V 1997 Quantum Systems: New thends and methods (Proceeding, Minsk, 1996) ed Y S Kim et al (Singapore: World Scientific) pp [2] Kudryashov V V 1998 Doklady of the National Academy of Sciences of Belarus 42 (No 6) 45 (in Russian) [3] Bender C M and Orszag S A 1978 Advanced Mathematical Methods for Scientists and Engineers (New York: McGraw-Hill) [4] Delabaere E, Dilinger H and Pham F 1997 J. Math. Phys [5] Bender C M, Olaussen S A and Wang P S 1977 Phys. Rev. D [6] Barclay D T 1994 Phys. Lett. A [7] Romanovski V G and Robnik M 2000 J. Phys. A: Math. Gen [8] Voros A 1999 J. Phys. A: Math. Gen [9] Abramovitz M and Stegun I A 1970 Handbook of Mathematical Functions (New York: Dover) [10] Langer R E 1937 Phys. Rev [11] Landau L D and Lifshitz E M 1977 Quantum Mechanics (Oxford: Pergamon)