arxiv:quant-ph/ v1 23 Dec 1996

Transcription

1 Displaced and Squeezed Number States quant-ph/ LA-UR Michael Martin Nieto 1 arxiv:quant-ph/961050v1 3 Dec 1996 Theoretical Division Los Alamos National Laboratory University of California Los Alamos, New Mexico 87545, U.S.A. and Abteilung für Quantenphysik Universität Ulm D Ulm, GERMANY ABSTRACT After beginning with a short historical review of the concept of displaced (coherent) and squeezed states, we discuss previous (often forgotten) work on displaced and squeezed number states. Next, we obtain the most general displaced and squeezed number states. We do this in both the functional and operator (Fock) formalisms, thereby demonstrating the necessary equivalence. We then obtain the time-dependent expectation values, uncertainties, wave-functions, and probability densities. In conclusion, there is a discussion on the possibility of experimentally observing these states. 1 mmn@pion.lanl.gov

2 1 Background The coherent states were discovered in 196 by Schrödinger 1], as an example of how his wave functions could mimic classical particles. They consisted of Gaussians in an harmonic oscillator potential. The widths of these Gaussians were those of the ground-state Gaussian. In modern notation, they were ψ = π 1/4 exp (x x 0) + ip 0 x ]. (1) It is to be noted that, this being before Born s probability interpretation of ψ ψ, Schrödinger was concerned by the complex nature of the wave function, wondering if it was the real part that was significant. Lorentz was the person who had been especially bothered by the lack of classical properties of eigenstate wave functions. Schrödinger and Lorentz exchanged many letters on this subject, and their correspondence is published ]. Soon after this, Kennard 3] published a paper on quantum motion. There he described what are, in modern parlance, squeezed states. They follow the classical motion, they are Gaussians whose widths are not that of the ground state, and the widths and uncertainty products oscillate with time. However, relatively little notice was given to this paper in later times 4]. In the 1960 s the emergence of quantum optics was a fertile background for the modern development of coherent states by Glauber, Klauder, and Sudarshan 5]-8]. Using boson operator calculus, these states could be defined as either the eigenstates of the destruction operator, a α = α α, () or equivalently as the state obtained by operating on the ground state by the displacement operator: D(α) 0 = α, (3) D(α) = expαa α a] = exp α /] expαa ] exp α a], (4)

3 the last equality coming from a Baker-Campbell-Hausdorff (BCH) calculation. These states are α = e α / n=0 α n n! n. (5) Use of a generating formula for Hermite polynomials shows that this yields the same wave function as Eq. (1), with the identifications α = α 1 + iα = x 0 + ip 0. (6) Squeezed states were rediscovered and elucidated by a number of people 9]-1]. (Ref. 11] was where the term squeeze was invented. See Ref. 13] for a discussion of this history.) In operator form, these states are created by D(α)S(z) 0 = α, z, (7) where the squeeze operator is 1 S(z) = exp za a 1 ] z aa, z = re iφ = z 1 + iz. (8) By BCH relations the squeeze operator can be written as ] ( ) 1 1 ( 1 S(z) = exp eiφ (tanh r)a a +a a) exp 1 ] cosh r e iφ (tanh r)aa (9) ] 1 = exp eiφ (tanh r)a a (cosh r) 1/ n=0 exp 1 ] e iφ (tanhr)aa The functional forms of D and S are 14] (sechr 1) n (a ) n (a) n n!. (10) where S = S 1/ exp D(α) = exp ix 0 p 0 /] expip 0 x] exp x 0 ], (11) ] ] iz sinh r iz sinh r r S (x ) exp (ln S)(x )] exp r S ( ), (1) S = cosh r + z 1 r sinh r = cosh r + cosφsinh r = er cos φ + e r sin φ, (13) 3

4 and one should recall the operator definitions on a function h: is 14] expc ]h(x) = h(x + c) (14) expτ(x )]h(x) = h(xe τ ) (15) ] expc( 1 (y x) )]h(x) = exp h(y)dy. (16) 4πc] 1/ 4c Using these functional forms for D and S the most general squeezed wave function ψ ss = D(α)S(z)ψ 0 = 1 exp ix 0 p 0 /] π 1/4 S(1 + iκ)] 1/ exp ( ) ] (x x 0 ) 1 S (1 + iκ) iκ + ip 0 x,(17) where κ z sinh r rs. (18) Setting z to be real and positive, yields the most commonly studied example: ψ ss = π 1/ s] 1/ exp (x x 0) s ip 0 x ], s = e r. (19) The Forgotten Displaced Number States In the 1950 s, before the modern work on coherent states started, there was a short flurry of activity by four authors, which has essentially been forgotten. These papers were by Senitzky 15] (perhaps the best known), Plebanski 16], Husimi 17], and Epstein 18]. These authors investigated whether there are other wave packets, besides the coherent-state packets, which keep their shapes and follow the classical motion. They found that there are. In our modern parlance, they found that any displaced number state follows the classical motion and keeps its shape. In particular, given that the number-state wave functions are ψ n = exp x ] 4 H n (x), (0) π 1/4 n n!] 1/

5 the displaced number states are found by applying the functional form of the displacement operator onto the number wave functions: ψ cs(n) = D(α)ψ n = exp ix 0p 0 /] expip 0 x] π 1/4 exp (x x 0) An equivalent result can be obtained in Fock notation by writing ] Hn (x x 0 ) n n!] 1/. (1) D(α) n = exp α /] expαa ] exp α a] n () = exp α α k n ( α ) j ] 1/ (n j + k)!n! /] n j + k.(3) k! j! (n j)!(n j)! k=0 j=0 Then by changing from the Fock state basis to the wave function basis, interchanging the order of the summations, and using the sum of Ref. 19], one has D(α) n π 1/4 exp (x α 1 ) / + i α x iα 1 α ] (n!) 1/ n ( α ) j j/ n/ j!(n j)! H n j(x α/ ). (4) j=0 Changing the sum over j to a sum over k = n j, the sum is in the form of the seventh equality on page 55 of Ref. 0]. Using this one obtains the same answer as Eq. (1) above. Time-evolution can be shown to be described by 14] x 0 x 0 (t) = (x 0 cost + p 0 sin t), p 0 p 0 (t) = (p 0 cost x 0 sin t). (5) The wave packets have the properties x(t)] = p(t)] = (n + 1/), (6) H = (n + 1/) + α. (7) In the later excitement over ordinary coherent and squeezed states, the work of the earlier authors 15]-18] was, unfortunately, mainly forgotten. In 1973 one paper appeared apllying the operator-formalism displacement operator on the Fock number state 1]. Then, in the 1980 s, papers started appearing studying these states ]- 4], one of which ] noted the early work by Senitzky and others 15]-18]. 5

6 3 Squeezed Number States Plebanski 16], in wave-function form, also looked a little into what we would call squeezed number states. By using the functional forms in Eqs. (11) and (1) of D and S 14], we find that the most general squeezed number states are where ψ ss(n) = D(α)S(z)ψ n = exp ix 0 p 0 /] π 1/4 F 1 ] 1/ exp (F 3 ) n/ 1 n n!] 1/H n (x x 0) ] F + ip 0 x ( ) ] x x0, (8) F 4 F 1 = S(1 + iκ) = cosh r + e iφ sinh r, (9) ( ) ( ) ( )] 1 F = S (1 + iκ) iκ 1 1 = iκ 1 + S (1 + 4κ ) S (1 + 4κ ) = 1 i sin φ sinh r(cosh r + e iφ sinh r) (cosh r + cosφsinh r)(cosh r + e iφ sinh r), (30) F 3 = 1 iκ 1 + iκ = ( cosh r + e iφ ) sin φ sinh r = e itan 1 (κ), (31) cosh r + e iφ sin φ sinh r F 4 = S(1 + 4κ ) 1/ = cosh r + sinh r + cosφ cosh r sinh r] 1/. (3) By using properties of the F s. such as that F 3 is a phase, F + F ] = /F 4, and F1 F 1 = F4, one can verify that 1 = dx ψ ss(n) (x) ψ ss(n)(x). (33) Note that when n = 0, the final large square bracket of Eq. (8) is unity and the wave function reduces to the normal squeezed states. Satyanarayana 3] defined the problem in Fock notation. He obtained that the squeezed number states are α, z, n = D(α)S(z) n 6

7 = exp α /] m,l G mn (z) ( ) 1/ m! L ( m) (l m) ( α )α l m l, (34) l! where the G mn (z) are given by infinite sums. Later discussions followed. One emphasized number states having the displacement operator and squeeze operator being applied separately (rather than in combination) 5]. Another was interested in these states from the point of view of the parity operator 6]. A third still addressed the topic from the viewpoint of molecular wave packets 7]. Now we will show how to sum the Fock problem in closed form, thereby reproducing the result of Eq. (8). Start with a number state and apply the squeeze operator. If then d = 1 eiφ tanh r, (35) ( )1 S(z) n = e da a 1 +a a e d aa n (36) cosh r ( )1 = e da a 1 +a a n ] cosh r j=0 n! ] 1/ ( d ) j (n j)! j! n j (37) = ( ) 1 n+1/ (n!) 1/ n ] ( d ) j (cosh r) j cosh r j=0 (n j)!j! d k (n j + k)!] 1/ n j + k, (38) k=0 k! where, in the upper limit of the sums, n/] is the greatest integer function. With the identification n j + k and using the known infinite sum of Ref. 19], k=0 e x / H n j+k (x) π 1/4 n j+k (n j + k)!] 1/, (39) v k ] ( ) 4vx k! H k+m(x) = (1 + 4v) m/ 1/ x exp H m, (40) 1 + 4v 1 + 4v] 1/ 7

8 one finds where S(z) n exp x F /](n!) 1/ π 1/4 F 1/ 1 cosh r(cosh r + e iφ sinh r)] n/ n/r, (41) and R = n ] τ j (n l)!j! H n j(xy), (4) j=0 τ = e iφ sinh r(cosh r + e iφ sinh r), y = cosh r. (43) F 1/ 1 To do the sum in Eq. (4), we start with the n is even case, n/] = n/. By changing the variable of the sum to k = n/ j and doing a little algebra, one finds R E = τn/ n/ ( n/) k (n/)! k=0 (k)! But this sum is of the form of equation of Ref. 19]: n k=0 ( ) 1 k H k(xy). (44) τ ( n) k (k)! vk H k (x) = n! (n)! ( 1 v)n H n x(1 + 1/v) 1/ ]. (45) Using this and a fair amount of further algebra one finally obtains for the n being even case (n = n E ), S(z) n = n E ] 1 exp x π 1/4 F 1 ] 1/ F (F 3) n/ ( ) x n n!] 1/H n. (46) F 4 But this is exactly the form of Eq. (8) for the special case of no displacement, x 0 = p 0 = 0. Therefore, now applying the displacement operator trivially yields the general result of Eq. (8), as it should. The n being odd case n = n O is obtained in the same manner. The only significant difference is that instead of using the sum of Eq. (45), one uses Eq of Ref. 19]: n k=0 ( n) k (k + 1)! vk H k+1 (x) = n! (n + 1)! ( v)n (1+1/v) n+1/ H n+1 x(1+1/v) 1/ ]. (47) 8

9 The identical functional form is obtained and so we have demonstrated, as it must be, that the functional and operator formalisms yield the identical squeezed number states. 4 Time-Dependent Uncertainties Calculating the time-dependent uncertainties is most easily done with the bosonoperator formalism. This is because we know that D (α)ad(α) = a α, (48) S (z)as(z) = (cosh r)a + e iφ (sinh r)a, (49) T (t)at(t) = ae it, (50) where T is the time-evolution operator T(t) = exp iht] = exp i(a a + 1/)t] = exp i( + x )t/]. (51) With this it is straight-forward to calculate the uncertainties, as a function of time, for the squeezed number states. To begin, x s (t) n = n S D T xtds n = x 0 cos t + p 0 sin t. (5) Similarly, p s (t) n = p 0 cost x 0 sin t, (53) x s (t) n = (n + 1)(cosh r) + (sinh r) + (cosh r)(sinh r) cos(t φ)] + x s (t) n, (54) p s (t) n = (n + 1)(cosh r) + (sinh r) (cosh r)(sinh r) cos(t φ)] + p s (t) n. (55) 9

10 Therefore, the uncertainty-product as a function of time is x(t)] p(t)] (n + 1/) = 1 + 4(cosh r) (sinh r) sin (t φ) (56) = ( s 1 s ) sin (t φ). (57) Eq. (56) agrees with the results of Ref. 6] for t = 0 and Eq. (57) agrees with the standard squeezed-state result when n = 0 8]. 5 Time-Evolution of the States Taking the functional definition of the time-displacement operator, T, in Eq. (51), and using BCH relations it has been shown 14] that the harmonic-oscillator timedisplacement operator can be written as T = cos t] 1/ exp i ] ] i tant(x ) exp (ln cost)(x )] exp tan t( ). (58) Using this result, along with the functional definitions in Eqs. (14-16), and combining it with the squeezed state wave function, ψ ss(n) of Eq. (8), one has Ψ ss(n) (x, t) = Tψ ss(n) (x) = e i tan t x / e ix 0p 0 / F n/ 3 I, (59) π 3/4 i sin tf 1 n+1 n!] 1/ where I = dy exp (y x/ cost) i tant 1 (y x 0) F + ip 0 y ] ( ) y x0 H n. (60) F 4 This integral can be evaluated with the aid of of Ref. 9], yielding Ψ ss(n) (x, t) = ] F n 1 3 A n ] 1/ ( ) X(t) H π 1/4 BF 1 ] 1/ n n n! F 4 BA] 1/ exp exp x p 0 F cost + i sin t B ( ) i sin t B + x (x 0F + ip 0 ) cost B x 0 ] F cost B ( F sin t x 0 p 0 B + i ) ], (61) 10

11 where the F i are given in Eqs. (9) - (3), and ( ) ( ) i sin t B i sin t/f A = 1 4, (6) F4B B B = cost + if sin t, (63) X(t) = x (x 0 cost + p 0 sin t). (64) Appropriate limits can be checked. For the case t = 0. one finds Ψ ss(n) (x, t = 0) = ψ ss(n) (x), the function given in Eq. (8). Also, for the non-squeezed case of z = 0 or F i = 1, one finds the appropriate time-evolution of displaced number states: e i(n+1/)t ] Ψ ss(n) (x, t, F i = 1) = π 1/4 n n1] 1/H n(x(t)) exp X (t) The probability density is given by expi{x (x 0 cost + p 0 sin t)/}(p 0 cost x 0 sin t)]. (65) ρ ss(n) (x, t) = Ψ ss(n) (x, t)ψ ss(n)(x, t) (66) AA ] n/ ( ) ( ) X(t) X(t) = H n (n!)πbb ] 1/ n H F 4 F 4 B A ] 1/ n F 4 BA] 1/ ] exp X (t). (67) F4BB With the aid of integral of Ref. 30], dx e px H n (bx)h n (cx) = n 1 n! π 0 p (n+1)/ (b +c 3 p) n/ bc P n, (68) p(b + c p) where the P n are the Legendre functions, it can in principle be demonstrated that the probability density is properly normalized: 1 = dx ρ ss(n) (x, t). (69) However, for small n it is easier to verify the normalization on a case by case basis. The first three probability densities are ρ ss(0) (x, t) = 1 πbb ] 1/ F 4 exp 11 X BB F4 ], (70)

12 ρ ss(1) (x, t) = ρ ss() (x, t) = ] ( ) exp X X, (71) πbb ] 1/ F 4 BB F4 BB F4 AA ] exp X πbb ] 1/ F 4 BB F4 ( X B A F 4 1 ) ( X B AF 4 1 ). (7) ρ ss(0) (x, t) is the probability density for the ordinary squeezed state. ρ ss(1) (x, t) and ρ ss() (x, t) are n = 1 and n = squeezed number states, having two and three humps, respectively, in their wave packets. In Figs. 1-4 we show three-dimensional plots of ρ ss(1) vs. x and t for various n = 1 squeezed states. Fig. 1 looks at the case x 0 = 8, p 0 = 0, and z = ln. The picket starts with no velocity at maximum displacement. It is at its broadest when t = 0, π, and π for x = 8, 8, and 8. In between the locations where the wave packet is broadest, the two humps are closer together, narrower, and higher in amplitude. In Fig, we show the case again having x 0 = 8 and p 0 = 0, but this time z = ln. This basically changes the phase by π/. The narrow, peaked wave packets are when t = 0, π, and π for x = 8, 8, and 8. The broader and shorter conditions are in between. Fig. 3 is again an example with x 0 = 8 and p 0 = 0, but now z is imaginary: z = (ln ) expiπ/] = i ln. One can see from the t = π case that now we have a phase shift such that neither extreme configuration is located at the t = an integral number times π positions. Finally, in Fig. 4, we change to a central position and maximum velocity to the right at t = 0: x 0 = 0 and p 0 = 8, with z = ln. Now the broadest wave packets are centered at x = 0 for t = 0, π, and π. The narrow peaked packets are when t = π/ and 3π/. 6 Discussion There is hope that, in the not too distant future, displaced and squeezed number states can be observed. This optimism is based on the recent work of Wineland s 1

13 group with trapped Be + ions 31, 3]. With laser cooling and a series of laser pulses to entangle the electronic and motional states they produced (displaced) even and odd coherent states 31, 33]. They also have been able to produce a squeezed (but not displaced) ground state and number states 3]. If all these techniques can be combined then in principle displaced and squeezed number states perhaps can be produced. However, this would be no easy feat since there are very complicated heating problems involved. When these states are produced, they should mimic the features displayed in our Figs. 1 to 4. Acknowledgments I gratefully acknowledge conversations with Iwo and Sophie Bialynicki-Birula, Peter Knight, Wolfgang Schleich, and especially with A. Vourdas which took place at the Humboldt Foundation Workshop on Current Problems in Quantum Optics, organized by Prof. Harry Paul. This work was supported by the U.S. Department of Energy and the Alexander von Humboldt Foundation. References 1] E. Schrödinger, Naturwiss. 14 (196) 664. ] E. Przibram, Letters on Wave Mechanics, translated into English by M. J. Klein (Philosophical Library, NY, 1967). 3] E. H. Kennard, Zeit. Phys. 44 (197) 36. See Sec. 4C. 4] D. S. Saxon, Elementary Quantum Mechanics, Preliminary Edition, (Holden- Day, San Francisco,1964), p ] R. J. Glauber, Phys. Rev. 130 (1963)

14 6] J. R. Klauder, Annals of Physics 11 (1960) 13; J. Math. Phys. 4 (1963) 1055; J. Math. Phys. 4 (1963) ] E. C. G. Sudarshan, Phys. Rev. Lett. 10 (1963) 7. 8] J. R. Klauder and E. C. G. Sudarshan, Fundamentals of Quantum Optics (Benjamin, NY, 1968). 9] H. Takahasi, Ad. Communication Systems 1 (1965) 7. 10] H. P. Yuen, Phys. Rev. A 13 (1976) 6. 11] J. N. Hollenhorst, Phys. Rev. D 19 (1979) ] C. M. Caves, J. S. Thorne, R. W. P. Drever, V. D. Sandberg, and M. Zimmerman, Rev. Mod. Phys. 5 (1980) ] M. M. Nieto, in: Frontiers of Nonequilibrium Statistical Physics, eds. G. T. Moore and M. O. Scully (Plenum, New York, 1986) p ] M. M. Nieto, Quantum Semiclass. Optics 8 (1996) ] I. R. Senitzky, Phys. Rev. 95 (1954) ] J. Plebanski, Phys. Rev. 101 (1956) 185; Bull. Acad. Polon. 11 (1954) 13; ibid. 14 (1955) ] K. Husimi, Prog. Theort. Phys. 9 (1953) ] S. T. Epstein, Am. J. Phys. 7 (1959) ] E. R. Hansen, A Table of Series and Products, (Prentice-Hall, Englewood Cliffs, N.J., 1975). 0] W. Magnus, F. Oberhettinger, and R. P. Soni, Formulas and Theorems for the Special Functions of Mathematical Physics, (Springer, New York, 1966). 14

15 1] M. Boiteux and A. Levelut, J. Phys. A 6 (1973) 589. ] S. M. Roy and V. Singh, Phys. Rev. D 5 (198) ] M. V. Satyanarayana, Phys. Rev D 3 (1985) ] F. A. M. de Oliveira, M. S. Kim, P. L. Knight, and V. Bu zek, Phys. Rev. A 41 (1990) ] M. S. Kim, F. A. M. de Oliveira, and P. L. Knight, Opt. Comm. 1 (1989) 99; Phys. Rev. A 40 (1989) 494; J. Mod. Opt. 37 (1990) ] R. F. Bishop and A. Vourdas, Phys. Rev. A 50 (1994) ] J. P. Dowling, W. P. Schleich, and J. A. Wheeler, Ann. Phys. (Leipzig) 48 (1991) 43. 8] M. M. Nieto, Quantum Optics 6 (1994) 9. The quantity in the large round brackets of Eq. (5) should be squared. 9] I. S. Gradshteyn and I. M. Ryzhik, Tables of Integrals, Series, and Products, (Academic Press, N.Y., 1965). 30] A. P. Prudnikov, Yu. A. Brychkov, and O. J. Marichev, Integrals and Series, Vol., (Gordon and Breach, New York, 1986). 31] C. Monroe, D. M. Meekhof, B. E. King, and D. J. Wineland, Science 7 (1996) ] D. M. Meekhof, C. Monroe, B. E. King, W. M. Itano, and D. J. Wineland, Phys. Rev. Lett. 76 (1996) ] Also see the analyses of these states in R. L. de Matos Filho and W. Vogel, Phys. Rev. Lett. 76 (1996) 608; M. M. Nieto, Phys. Lett A 19 (1966)

16 Figure Captions Figure 1. A three-dimensional plot of the n = 1 squeezed-state probability density ρ ss(1) as a function of position, x, and time, t, for x 0 = 8, p 0 = 0, and z = ln. The spikes near the maxima of ρ ss(1) are an artifact of the numerical routine. The true, smooth maxima follow the tops of the spikes. Figure. A three-dimensional plot of the n = 1 squeezed-state probability density ρ ss(1) as a function of position, x, and time, t, for x 0 = 8, p 0 = 0, and z = ln. Figure 3. A three-dimensional plot of the n = 1 squeezed-state probability density ρ ss(1) as a function of position, x, and time, t, for x 0 = 8, p 0 = 0, and z = (ln ) expiπ/] = i ln. The spikes near the maxima of ρ ss(1) are an artifact of the numerical routine. The true, smooth maxima follow the tops of the spikes. Figure 4. A three-dimensional plot of the n = 1 squeezed-state probability density ρ ss(1) as a function of position, x, and time, t, for x 0 = 0, p 0 = 8, and z = ln. 16

17 This figure "sns1.gif" is available in "gif" format from:

18 This figure "sns.gif" is available in "gif" format from:

19 This figure "sns3.gif" is available in "gif" format from:

20 This figure "sns4.gif" is available in "gif" format from: