Arxiv:1609.07380V1 [Quant-Ph] 23 Sep 2016 Where Nntl Ag Oeta Nry Uha H Ia Et Ucinpo Function Delta Conﬁ [1]

arXiv:1609.07380v1 [quant-ph] 23 Sep 2016 where nntl ag oeta nry uha h ia et ucinpo function delta confi [1]. Dirac for models the model as these standard such on a energy, is potential and large mechanics, infinitely quantum in model basic etlmdl ncnrs oteDrcdlafnto oeta,th potential, as function notation delta the Dirac using Simply the to described. contrast clearly In model. mental hoe nti oe ss a ak h esni htw a o ca not can we that is reason The lack. term a the far to so corresponds is which model operator, this in theorem rcs aclto o xetto aus nti ae,w invest we the confirmed. paper, of be this infiniteness In directly of can values. degree expectation the for describe calculation to precise formula precise more a h ieidpnetSh¨dne qainis Schrödinger equation independent time The ∗ For ehave We lcrncades [email protected] address: Electronic uiul,frifiiesur el esild o aeadrc manifes direct a have not do still we well, square infinite for Curiously, huhw aeeatsltosadteesltosaesml;howe simple; are solutions these and solutions exact have we Though h oeta energy potential The h nnt qaewl samdluigifiieyhg oeta barrier potential high infinitely using model a is well square infinite The V ( ietmnfsaino hefs’ hoe nteinfinite the in theorem Ehrenfest’s of manifestation Direct k x n endi q 1,teeeg iefntosΨ eigenfunctions energy the (1), Eq. in defined ) = nπ/L qaewl stelmto h nt el o ietmanife direct a For diverge the well. of infiniteness finite term, of the term squar degree energy of finite potential the the limit describe in to the done formula as first well is theorem square this of manifestation ewrs nnt qaewl;Ernetterm ieind condition; Time Boundary theorem; function; Ehrenfest delta well; Dirac solu square 03.65.Ge, Infinite mechanics; Keywords: Quantum 03.65.-w, numbers: PACS but vague that not precisel is calculated be model directly. expres can this values be in expectation can that infinity results which the function that means defined This well a is eigenfunction, hefs’ hoe nteifiiesur eli pt o on now to up is well square infinite the in theorem Ehrenfest’s , n 1 = V , 2 ( x , 3 fteifiiesur eli ecie by described is well square infinite the of ) ... . V E Ψ ( x − n n − )Ψ ( = 2 dV x ~ n m V = ) 2 ( .INTRODUCTION I. ~ x (ˆ ( Ψ eateto Physics, of Department 2 ∞ apiTia,R.O.C. Taipei,Taiwan, x x ,wihi h rdc fteptnileeg n h energy the and energy potential the of product the is which ), k oco University, Soochow ′′ ) = ) ( n 2 /d ( hiLn Lin Chyi-Lung svge ehv ogieo o ohnl hsifiiy eneed We infinity. this handle to how on guide no have we vague, is / x 0 q 2 + ) x ˆ , m, ent httefrei nntl ag ttetosdsof sides two the at large infinitely is force the that note We . L 2 0 ∞ V n , sin( ( ( , x x 0 otherwise )Ψ( n h eigenvalues the and ) k L, < x < n x x n hefs’ hro a econfirmed be can thereom Ehrenfest’s and y ) pnetSh¨dne qain tpfunction; Step Schrödinger equation; ependent = ) gt hspolmadso htErnettheorem Ehrenfest that show and problem this igate , ∗ ino aeeutos on states bound equations: wave of tion iegn oeta nry ota ecnperform can we that so energy, potential divergent 0 otherwise eta.BloiadRbnt aegvnareview a given have Robinett and Belloni tential. nnt nteifiiesur elhsntbeen not has well square infinite the in infinity e tptnileeg.W hwta the that show We energy. potential nt E L, < x < igprils eas aeohrmdl using models other have also We particles. ning e ntrso ia et functions. delta Dirac of terms in sed . el n hncnie h infinite the consider then and well, e Ψ( a bandaseicto.This specification. a obtained has e,tedrc aietto fEhrenfest’s of manifestation direct the ver, aino hefs’ hoe nti funda- this in theorem Ehrenfest’s of tation x tto,w edamr precise more a need we station, ) ocnn atce nieawl.Ti sa is This well. a inside particles confine to clt h xetto au fteforce the of value expectation the lculate . ymnfse niety The indirectly. manifested ly . E n fE.()aewl-nw [2–8]. well-known are (2) Eq. of qaewl model well square (2) (4) (3) (1) 2 the well, while the probability density Ψ∗(x)Ψ(x) is zero there. Hence, calculating the expectation value of the force operator, we encounter the ambiguity as the product of infinity and zero [9]. To solve this ambiguity, one may consider the infinite square well as a limiting case of a finite well [9, 10]. Rokhsar calculated the expectation value of the force in a finite well, and then took the limit of the height of the potential energy V0 to confirm Ehrenfest theorem in the infinite square well [9]. We note→ that ∞ the difficulty for a direct manifestation lies in lacking informations about the specification of the infinity used in the potential energy. For the case of Dirac delta function δ(x), whcih is usually defined as

0, x =0, δ(x)= (5) , x 6=0. ∞ The infinity in the delta function is specified by the formula ∞ δ(x)dx =1. (6) Z−∞ Eq. (6) is the specification formula for the delta function. This formula shows how to handle the infinity in the delta function. In the same spirit, we seek a specification formula for the V (x) in the infinite square well model. We discuss this in Sec. II. We show the direct manifestation of Ehrenfest theorem in Sec. III, and we discuss more about the function form of V (x)Ψ(x) in Section IV.

II. THE DERIVATION OF THE FUNCTIONAL FORM OF V (x)Ψn(x)

To explore such a specification formula of V (x), we work in an opposite way. We start from wave functions Ψn(x) and substitute these wave functions into the Schrödinger equation, we can then determine the corresponding functional form of the potential energy term. For the infinite square well, the wave function in Eq. (3) can be rewritten in a more compact form as

Ψn(x)= un(x)θ(x)θ(L x). (7) − 2 un(x)= sin(knx). (8) rL As we are to derive the form of V (x) from the wave function, we then rewrite the Schr¨odinger equation as: ~2 V (x)Ψ(x)= Ψ′′(x)+ EΨ(x). (9) 2m For Ψ(x) of the form as that in Eq. (7), the Ψ′′(x) can be calculated by the following formula Ψ(x)= f(x)θ(x)θ(L x) − Ψ′′(x)= f ′′(x)θ(x)θ(L x)+ f ′(x)δ(x) f ′(x)δ(L x). (10) − − − Substituting Eq. (7) into Eq. (9) then yields

~2 ′ V (x)Ψn(x)= δ(x) δ(L x) U (x). (11) 2m − − n h i Substituting Eq. (8) into Eq. (11), we have

2 ~2 V (x)Ψn(x)= kn δ(x) cos(knL)δ(L x) . (12) L 2m − − r h i We have then obtained the functional form of V (x)Ψn(x). It had already been argued that the potential energy term, V (x)Ψn(x), should contain a delta function at each side of the well [10, 11]. Eq. (12) shows this result explicitly. Eq. (12) is an amazing formula. If we look back into the definition of V (x) in Eq. (1), then we see that the infinity in the V (x) is not that unmanageable. That is, although V (x) is a divergent quantity outside the well, however V (x)Ψn(x) is indeed a well-defined function which is described in terms of Dirac delta functions, as shown in the right side of Eq. (12). This result should be a stepping stone for further investigations on the infinite square well model when there are ambiguities in calculating expectation values. One of the application of this result is that we can directly confirm Ehrenfest’s theorem. 3

III. THE VERIFICATION OF EHRENFEST’S THEOREM IN THE INFINITE SQUARE WELL

We now show that Ehrenfest’s theorem for time-evolved wave packets in the inﬁnite square well can be manifested. The time evolution of a general wave packet Ψ(x, t) is as follows

∞ −iωnt Ψ(x, t)= anΨn(x)e . (13) n=1 X ∞ 2 where ωn = En/~, and an = 1. To verify Ehrenfest’s theorem, we need to verify the following formula: n=1 | | P d dV (ˆx) Ψ(t) pˆ Ψ(t) = Ψ(t) Ψ(t) . (14) dth | | i −h | dxˆ | i

The calculation of the right side of Eq. (14) is related to the calculation of Ψn(x)(dV (x)/dx)Ψj (x). We note that

dV (x) Ψ (x) Ψ (x) n dx j d dΨn(x) dΨj (x) = [Ψn(x)V (x)Ψj (x)] [V (x)Ψj (x)] [Ψn(x)V (x)] . (15) dx − dx − dx For the right side of Eq. (15), the ﬁrst term makes no contribution when taking the integration. This is because after the integration of the ﬁrst term over the range of x , [0,L], the result is the value of the boundary term Ψn(x)V (x)Ψj (x) calculated at the two sides of the well. And this yields zero by using Eq. (12) and Eqs. (7-8). Using again Eq. (12), the other two terms can be calculated, and we obtain

~2 dΨn(x) n+j V (x)Ψj (x)= knkj θ(x)θ(L x)[δ(x) ( 1) δ(L x)]. (16) dx mL − − − − And then we have ∞ dΨ (x) n V (x)Ψ (x)dx dx j Z−∞ ~2 L n+j = knkj [δ(x) ( 1) δ(L x)]dx mL − − − Z0 ~2 = k k β . (17) 2mL n j nj

n+j where βnj =1 ( 1) . Above, we have used the following results − −

L 1 δ(x) dx = , (18) 2 Z0 L 1 δ(L x) dx = . (19) − 2 Z0 These results are due to the even function property of the delta function. We then have

dV (ˆx) ~2 Ψn(x) Ψj (x) = knkj βnj . (20) h | dxˆ | i −mL And then we obtain the ﬁnal result dV (ˆx) Ψ(t) Ψ(t) h | dxˆ | i ∞ dV (x) = Ψ∗(x, t) Ψ(x, t)dx dx Z−∞ ~2 ∞ ∞ ∗ i(ωn−ωj )t = a aj knkj βnj e . (21) −mL n n=1 j=1 X X 4

We also easily have the result

Ψ(t) pˆ Ψ(t) h | | ∞ i ∞ 2 ∗ knkj i(ωn−ωj )t = ( i~) a aj βnje . − L n k2 k2 n=1 j ,j6 n n j X =1X= − (22)

It then follows that d Ψ(t) pˆ Ψ(t) dth | | i ~2 ∞ ∞ = a∗ a k k β ei(ωn−ωj )t. (23) mL n j n j nj n=1 j=1 X X Comparing Eqs. (21) and (23), we see that Ehrenfest’s theorem is confirmed. We have thus directly verified Ehrenfest’s theorem in the infinite square well. Our results, Eqs. (20), (21), are the same as those of Rokhsar [9]. For an example of Eq. (21) , we choose a1 = a2 = 1/√2, and an = 0, otherwise. Then from Eq. (21), we have Ψ(t) dV (ˆx)/dxˆ Ψ(t) = (8E /L) cos(ω t), where ω = (ω ω ). h | | i − 1 12 12 2 − 1

IV. MORE ABOUT THE FUNCTIONAL FORM OF V (x)

In what follows, we try to obtain a functional form of V (x) from Eq. (12). We note that Eq. (12) is not in a form that is symmetrical with respect to the two edges. We search for a more symmetrical one. We rewrite Eq. (12) as

~2 V (x)Ψ (x)= (A1+ A2). n 2m 2 A1= knδ(x). rL 2 A2= kn cos(knL)δ(L x). (24) −rL −

We can rearrange the two terms A1 and A2 expressed in terms of un(x). For A1, using knδ(x) = [sin(knx)/x]δ(x), we then have A1 = [un(x)/x]δ(x). We also need to rewrite A2 in a similar way. This can be done, as un(x) = 2/L sin(knx) can also be written as p 2 un(x)= sin[kn(x L)] cos(knL). (25) rL −

Above, we have used sin(knL) = 0 and cos(knL) = 1. Using knδ(L x) = [sin(kn(L x))/(L x)]δ(L x), then ± − − − − we have: A2 = [un(x)/(L x)]δ(L x). Finally, Eq. (12) can be rewritten as − − ~2 δ(x) δ(L x) V (x)Ψ (x)= + − u (x). (26) n 2m x (L x) n −

Substituting Ψn(x) in Eq. (7) into Eq. (26) then yields

~2 δ(x) δ(L x) V (x)θ(x)θ(L x)un(x)= + − un(x). (27) − 2m x (L x) −

As the functions un(x) form a complete set, Eq. (27) can be extended to a general function ψ(x), which is composed from un(x), that is

∞ ψ(x)= cnun(x). (28) n=0 X 5

Then we have ~2 δ(x) δ(L x) V (x)θ(x)θ(L x)ψ(x)= + − ψ(x). (29) − 2m x (L x) − This formula describes more precisely the property of the V (x) in Eq. (1). Eq. (29) together with Eq. (28) are important results, as they provide a more precise starting point for the infinite square well model. From Eq. (29), we see that, naively, V (x)θ(x)θ(L x) behaves like having a form as (~2/2m)[δ(x)/x + δ(L x)/(L x)]. The factor θ(x)θ(L x)− accompanied with V (x) seems unavoidable. This factor is in fact needed,− due− to the divergence nature of the− V (x) outside the well. We may be curious on how to describe the degree of the infiniteness of the V (x) in Eq. (1). Because V (x) is divergent outside the well, Eq. (29) shows that this divergence is eliminated when multiplied by the function θ(x)θ(L x), as we see that the right side of Eq. (29) is well-defined. Thus, the θ(x)θ(L x) term is needed to accompany− with the V (x) in order to have a regular result. In other words, the θ(x)θ(L −x) factor shows how to handle the infinity contained in the V (x). Eq. (29) may then be viewed as a specification− formula to specify the infinity in the V (x). Eq. (29) describes the property of the infinite barrier. This formula is obtained from the known eigenfunctions Ψn(x), which are continuous at the two sides of the well. There are other types of solutions of the infinite square well, in which the continuity of Ψ(x) at boundaries is not required [12]. We do not consider that case at present.

Acknowledgments

The author is indebted to Prof. Young-Sea Huang for bringing in the attention and interests of this subject, and also for much help in discussions. The author would also like to thank Prof. Tsin-Fu Jiang and Prof. Wen-Chung Huang for much help.

[1] M. Belloni and R. W. Robinett, ”The infinite well and Dirac delta function potentials as pedagogical, mathematical and physical models in quantum mechanics”, Physics Report, 540, 25-122 (2014) [2] Leonard I. Schiff, ”Quantum Mechanics”, 3rd edition.McGraw-Hill, New York (1968). [3] A. Messiah, ”Quantum Mechanics”,North-Holland, Amsterdam,(1964). [4] E. Merzbacher ”Quantum Mechanics”, 2nd edition, J. Wiley, New York (1970). [5] R. Shankar, ” Principle of Quantum Mechanics”, Plenum Press, New York (1980). [6] G. Baym, ”Lectures on Quantum Mechanics”, W. A. Benjamin Inc., New York (1969). [7] David J. Griffiths, Introduction to quantum mechanics, 2nd edition. Pearson Education, Inc., New Jersey (2005). [8] David H. Mclntyre, ”Quantum Mechanics”, Pearson Education, Inc., San Francisco (2012). [9] D. S. Rokhsar, ”Ehrenfest’s theorem and the particle-in-a-box”, Am. J. Phys., 64, 1416-1418 (1996). [10] Ryoichi Seki, ”On boundary conditions for an infinite square well potential in quantum mechanics”, Am. J. Phys., 39, 929-931 (1971). [11] H. Bethe and J. Goldstone, ”Effect of a Repulsive Core in the Theory of Complex Nuclei”, Proc. Roy. Soc. (London) A, 238, 551-567 (1957). [12] Young-Sea Huang, ”Reexamination on the problem of the infinite square well in quantum mechanics”, DOI: 10.13140/RG.2.1.1277.5843.