A Note on the Normalization of the Momentum Eigenfunctions and Dirac Delta Function Mehdi Hage-Hassan

A note on the normalization of the momentum eigenfunctions and Dirac delta function Mehdi Hage-Hassan

To cite this version:

Mehdi Hage-Hassan. A note on the normalization of the momentum eigenfunctions and Dirac delta function. 2014. hal-00978136v2

HAL Id: hal-00978136 https://hal.archives-ouvertes.fr/hal-00978136v2 Preprint submitted on 22 Apr 2014

HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés.

A note on the normalization of the momentum eigenfunctions and Dirac delta function

M. Hage-Hassan Université Libanaise, Faculté des Sciences Section (1) Hadath-Beyrouth

Abstract We determine the generating function of the harmonic oscillator by a new method. Using this generating function we derive the eigenfunctions of the moment p. We find that the normalization of these eigenfunctions is a real and not complex number with phase factor chosen equal one (standard books of quantum mechanics). We prove that the integral of the delta function is equal to one and we derive the oscillator propagator.

1. Introduction

Despite the success of quantum mechanics many problems remain unsolved or solved by difficult methods. Among these problems the normalization Np of the eigenfunctions of the moment which is assumed to be a complex number but « the phase factor is chosen to be equal to unity » [1.p:330, 2.p:54, 3.p:101]. We need also simple methods to determine: the generating function of the basis of the oscillator [3], the integration of delta function and the Feynman propagator of the oscillator [1-3]. We use the generating function method that we develop and that has enough variety of applications [4] to solve these problems. So we show that normalization is real and not arbitrary and that Np =√1/2πħ and the calculations of other problems occur without particular difficulty using the Gaussian integral. In the second part we give a quick review of the wave functions in position and momentum. The third and the fourth part are devoted to a revision of the linear harmonic oscillator and the derivation of the generating function of the oscillator. We devote the sixth part to the derivation of the wave function of the momentum. We treat the properties of Dirac delta function in part sixth. In part seven, we give the derivation of the Feynman propagator of the harmonic oscillator without the use of Mehler’s formula [1,2].

2. The wave functions in position and momentum

In quantum mechanics the wave function is expressed in terms of coordinate {x} or momentum{p}. A- In the coordinate representation we write: h d x, p = , and x,[ px ] = ih (2.1) i dx

Following Dirac the eigenfunctions of the observable xo can be written: x x xx , and ' xx x xx ' x x x'' x o = o = =

Therefore: (x − x )' x' x = 0 (2.2)

The normalization of this space is the Dirac delta function which is defined by: (x x )' x' x (2.3) δ − = B- In moment representation (p, x = hdi / dp ) the wave function is xp )( = x p , h d h d and px p = pp(x), x p = xp )( (2.4) i dx i dx h d The solution of the equation xp )( = pp x)( is: i dx xp )( N eixp / h (2.5) = p

We will determine the normalization Np in part five. C- We write the unitary operators in the two representations by:

I x dx x , I p dp p (2.6) = ∫ = ∫

3. The linear Harmonic oscillator

We will give only a summary of known results on the harmonic oscillator. 3.1 Schrödinger equation waves functions A- The Schrödinger equation of the harmonic oscillator in one dimension is: 2 2 2 h d ψ x)( mω Hψ x)( −= 2 + ψ x)( = Eψ x)( (3.1) 2m dx 2 The energies En= ħω (n+1,2) and the waves functions are [1,2,3]:

2 −q 2/ mω ψ n x)( = N n H n q)( e , q = x (3.2) h The Hn (q) are the Hermite polynomials of degree n and the normalization is: 2/1 ⎛ 1 mω ⎞ N n = ⎜ ⎟ (3.3) ⎜ 2n n! h ⎟ ⎝ π ⎠ The parity of the wave function is given by: n ψ n (−x) = (− )1 ψ n x)( . (3.4)

B- In momentum space x = ih∂ / ∂p the Schrodinger equation becomes: 2 2 p 2 2 d φ p)( φ p)( − h mω = Eφ p)( (3.5) 2m dp 2

The solution is:

2 n −t 2/ p φn p)( = α Cn H n )( et where t = (3.6) mωh With αn is the phase factor and

2/1 1 N C ⎛ ⎞ n (3.7) n = ⎜ n ⎟ = ⎝ 2 n! mωhπ ⎠ mω

3.2 The Harmonic oscillator in Dirac notation The Schrodinger equation for the harmonic oscillator is:

2 2 ⎛ p mω ⎞ Hψ x)( = ⎜ + ⎟ψ x)( = Eψ x)( (3.8) ⎜ 2m 2 ⎟ ⎝ ⎠

This equation is analogous to the product ( ⋅ zz ), z ∈ C , therefore we put:

mω 1 d + mω 1 d a = q + , a = q − 2h 2mωh dq 2h 2mωh dq where a ,a + are adjoint operators with:

+ + + ,[ aa ] = aa ],[,1 = [,0 a ,a ] = 0 (3.9) Comparing [a, a + ]= 1and []d /( dz), z = 1we simply deduct by analogy and using Dirac notation the orthonormal basis of the harmonic oscillator. It is well known that the basis of polynomials P(z) is: n d 1, z, ... …, z ,… and 1 = 0 dz a + n therefore 1,0 ,..., n = 0 ,.. and a 0 = 0 . n! with a n = n n − ,1 a + n = n n + 1 + So a and a are the ladders operators, m n = δ ,nm , and 0 is the vacuum state. The unitary operator of the harmonic oscillator basis is:

I n n = ∑n=0 . (3.10)

Using Dirac transformation [5], the wave function may be written as: a + n u n (q) = q n = q 0 (3.11) n! 4/1 h 4/1 And u (q) ⎛ ⎞ ( x) , n u t)( m h ( p) (3.12) n = ⎜ ⎟ ψ n α n = ( ω ) φ n ⎝ mω ⎠

4. The generating functions of the harmonic oscillator

We will find the expression of the generating function by a new simple method different from the other methods [1-3] and closely related to Dirac method. This function is the trace of a product of the base ψ n x)( and the Fock space of analytic functions.

4.1 The generating functions Using Dirac transformation [1-5] and (3.9) we find:

za+ 1 d q ae 0 = (q + )G qz ),( (4.1) 2 dq za+ Using also (3.9) we find: q ae 0 = zG qz ),( By comparison of the above expressions we derive: d G qz ),( = ( 2z − q)G qz ),( (4.2) dq The solution of this equation is: q 2 G qz ),( = c exp{( 2qz − ) ϕ+ (z)} (4.3) 2 To determineϕ z)( we use the creation operator, we find:

+ za+ 1 d ∂ aq e 0 = (q − )G qz ),( = G qz ),( (4.4) 2 dq ∂z WE the help of (4.4) and (4.3) we obtain: ' ϕ z)( −= z Therefore the generating function of the oscillator is: q 2 z 2 G qz ),( = c exp{ 2qz − − } 2 2

For z = 0, we have G q),0( = q 0 = u0 q)( and it follows that:

n 1 2 2 − ∝ z 4 q z Gx qz ),( = ∑n=0 un q)( π= exp{ 2qz − − } (4.5) n! 2 2 By analogy we find the generating function of the basis { φn p)( }

n 1 2 2 − ∝ (αz) 4 t (αz) G p z,( αt) = ∑n=0 un t)( π= exp{ t(2 αz) − − } (4.6) n! 2 2 n α)( is the phase factor of the momentum basis.

4.2. The Fock space n In the expression (4.5) the function f n z)( = z / n! constitutes a basis of analytic Hilbert space known by Fock space and defined by:

z m z n f m f n = f m )( fz n )( dz μ z)( = dμ z)( δ= ,nm (4.7) ∫∫ ∫∫ m! n!

dμ z)( is the cylindrical measure or Gaussian measure:

( vu 22 ) dμ z)( = π)/1( e− + dudv, z = u + iv (4.8)

5. The wave function of momentum

We want to determine the expression xp )( = x p using the generating function of the harmonic oscillator and the orthogonality of Fock space (4.7). We write: x p Ix p x n n p x)( ( p) = = ∑n=0 = ∑n ψ n φn =

1 1 z n ( z) m u q)( ( ) mu t)( u q)( α u )( dt z)( ∑ ,mn n × δ ,mn × α m = ∫ ∑ ,mn n m μ h h n! m!

1 G q ,( z )G t,( z )d z )( (5.1) = ∫ x p α μ h

Using the Gaussian integral

1 n 2/ 1 1 ⎛ ⎞ n t t −1 (5.2) ⎜ ⎟ ∏i=1 dxi exp(−x Xx + Jx) = exp( J X J ) 2 ∫ (det(X ) 2 ⎝ π ⎠ We find 2 2 ⎛3 + α i − iα ⎞ X ⎜ ⎟, = ⎜ 2 2 ⎟ ⎝i − iα 1 α− ⎠ We obtain using (5.2), (3.2) and (3,6):

2 2 2 2 2 2 1 ⎡q α + t α − 2tqα q t ⎤ xp )( = exp − − (5.3) ⎢ 1 2 2 2 ⎥ 2πh ⎣ − + α ⎦

2 2 if 2α −= 1+ α , α ±= i we find α = i if:

1 i xp )( = exp(+ xp) (5.4) 2πh h A simple calculation gives: 1 i ∫ exp(− pq)Gx xz ),( dq =G p z p),( 2πh h

We derive after the development of the above expression:

1 i φn p)( = ∫ exp(− xp)ψ n x)( dx (5.5) 2πh h

We proved that Np is real and not complex [1,2,3] and we find the same the phase factor of the momentum-space wave function (-i)n of ref.[1.p:144].

6. The Dirac delta function

A- With the help of (2,6) we find the expression of the delta function.

q' q q'( q) ' qIq δ= − = 1 = ∫ q ' p dp p q = ∫ exp( + ik q '( − q )) dk (6.1) 2 π B- Using the basis of the oscillator (3.10) and its generating function (4.5) we find a very useful expression of delta function:

q' q = ∑ q' n n q = ∑un q )'( un q)( (6.2) n n

We deduce also from (3.4) that the delta function is an even function. Repeating the same method above (5.1), we find

q' q = ∫ G x (z, )q' Gx qz ),( dμ z)(

Replacing Gx qz ),( and G x (z, )q' by the expression (4.5) we obtain 2 2 2 2 1 q' +q z + z q' q = ∫ exp[− − + '(2 zq + zq )]dμ z)( (6.3) π 2 2 The arrangement of this expression gives

1 1 2 2 2 q' q = exp(− (q − q )' ∫ exp([− (2 u − (q + q ))' ] + 2iv(q − q ))' dudv π π 4 4 Using the change of variables and performing the integration we find that

1 1 2 q' q = exp(− (q − q ))' ∫ exp(+ik(q − q ))' dk (6.4) 2π 4

Comparing (6.4) and (6.1) and with the help of (2.3) we find:

1 2 exp(+ (q − q ))' q' q = q' q δ= q'( −q) 4 C- With the help of the Gauss integral we find that

∝+ ∝+ ∝+ 1 1 2 ∫ δ q'( −q)dq'= ∫∫exp[− (q − q )' + ik(q − q )]' dq'dk = 1 (6.5) ∝− 2π ∝− ∝− 4

Therefore we find that the integral of the Dirac delta function is equal to one.

∝+ 1 ∝+ ∫ δ q'( −q)dq'= ∫ exp[+ik(q − q )]' dq'dk = 1 (6.6) ∝− 2π ∝−

Finally, we write:

1 q' q (q q )' ∝ e − )'( kqqi dk and ∝+ (q q )' dq' 1 (6.7) δ= − = ∫ ∝− ∫ ∝− δ − = 2π

So, there is no need to impose (6.6) as constraint [6].

7. The Harmonic oscillator propagator

Finally, we note that we obtain by the same simple calculation as above (5.1) the Feynman propagator, [1.p:164, 2.p:119], of the harmonic oscillator:

K((x,t), t,(x' )) (x,t) exp iH (t t /) h )t,(x' (7.1) 0 = [− − 0 ] 0

We write the propagator in term of the generating function:

m 2/1 ⎛ ω ⎞ ω 0 2/)t-(ti- ω 0 2/)t-(ti- ω 0 2/)t-(ti- K((x, t), t,(x' 0 )) = ⎜ ⎟ e ∫ Gx e( qz ), Gx e( qz )', dμ z)( (7.2) ⎝ h ⎠

with q = (mω /) h x , q'= (mω /) h x' and α = ω(t − t0 ) . The calculations of (7.2) are carried out quickly by using:

iα− iα− 2/ iα− iα− 2/ 1( + e ) = 2e cos(α ,)2/ 1( − e ) = 2e sin(α )2/ 2 ∝+ e −az dz a,/ with Re(a) .0 and e iα− 2/ /1 e iα− ∫ ∝− = π > = We obtain the propagator:

mω i 2 2 K((q,t), t,(q' 0 )) = exp[ [(q + q )' cos −α 2qq ]]' (7.3) 2πhisin α 2sin α

Consequently we do not encounter the difficulties of the method proposed by many others [1,2,7]. Therefore our method [4] can be used to calculate the Feynman propagator of charged harmonic oscillator in uniform magnetic field and many other problems.

References [1] E. Merzbacher, (1970). Quantum Mechanics, New York:Wiley [2] J. J. Sakurai and J. Napolitano “Modern Quantum Mechanics” (Addison Wesley 2011) [3] W. Greiner, “Quantum Mechanics (An introduction) “Ed. Springer 1994 [4] M. Hage-Hassan, “Generating function method and its applications to Quantum, Nuclear and the Classical Groups”, arXiv:1203.2892v1 [5] A. Messiah, Mécanique Quantique Tomes I et II Dunod 1965 [6] I.M. Gel'fand and G.E. Shilov, Generalized functions I, Academic Press, New. York, 1968. [7] F.A.Barone, H. C. Boschi-Filho, and Farina,’ Three methods for calculating the Feynman propagator’ Am. J. Phys., Vol. 71, No. 5, May 2003, 490.