arxiv: v2 [math-ph] 18 Mar 2011

Transcription

1 SU(, ) Coherent States For Position-Dependent Mass Singular Oscillators arxiv: v [math-ph] 8 Mar 0 Sara Cruz y Cruz and Oscar Rosas-Ortiz Sección de Estudios de Posgrado e Investigación, UPIITA-IPN, Av. IPN 580, CP México DF, Mexico Departamento de Física, Cinvestav, AP 4-740, México DF, Mexico Abstract The Schrödinger equation for position-dependent mass singular oscillators is solved by means of the factorization method and point transformations. These systems share their spectrum with the conventional singular oscillator. Ladder operators are constructed to close the su(, ) Lie algebra and the involved point transformations are shown to preserve the structure of the Barut-Girardello and Perelomov coherent states. Introduction There has been much interest in the construction of coherent states since the origin of quantum mechanics []. Following Glauber, the concept is useful in the quantum description of the correlation and coherence properties of light []. The Glauber states deal with the dynamical properties of the harmonic oscillator and posses three basic properties: they () are eigenfunctions of the annihilation operator belonging to complex eigenvalues, () are displaced versions of the ground state wave-function and (3) minimize the uncertainty relation between position and momentum. For systems other than the harmonic oscillator, such properties are used as different definitions of the coherent states and, in general, they are not equivalent to each other. Thus, the term coherent states (CS) has been used for a wide class of mathematical objects over the years. For instance, the generalized CS studied by Barut and Girardello [3] and by Perelomov [4] are based on the first and second properties of the Glauber states. In turn, the construction of CS as minimizing the uncertainty relation of a pair of observables is incidentally found as a secondary result for some special systems. Nevertheless, the most valuable property of the CS is that they can be studied for many systems in terms of the definition leading to the desirable result. The purpose of this work is to draw the CS applications to the study of positiondependent mass (PDM) systems in quantum mechanics. Some previous results include the CS for the harmonic [5 8] and nonlinear [8,9] oscillators. We are interested in the construction of PDM Hamiltonians such that their spectrum is exactly the same as that of the conventional singular oscillator. We shall show that the description of these new systems obeys the su(,) algebra. The structure of the paper is as follows. In Section, expressions for the mappings and the operators leading to the relevant algebras are given.

2 The algebraic structure of the PDM linear oscillators is revisited in Section 3. Departing from the ladder operators of the linear case, in Section 4 we construct a new class of operators and show that they are connected in a natural way to the singular oscillator. It is shown that these operators are the generators of the su(,) Lie algebra (the subject studied also in [0]). In Section 5 some applications are given and the corresponding CS are constructed a la Barut-Girardello and a la Perelomov. Some final remarks close the paper. Factorization and point transformations Consider a one-dimensional quantum system of mass m(x) acted by the potential V(x). A convenient expression for the Hamiltonian reads H a = ma Pm b Pm a +V, a+b = () where P and the position operator X satisfy [X,P] = i. This Hamiltonian admits the factorization H a = AB +ǫ, A = B = i m a Pm b +β, () with ǫ a constant (in energy units) to be fixed. In the position representation (X = x and P = i d ) the function β is a root of the Riccati equation: dx V ǫ = ( [ a+ )( ] m )β β +β, d m 4 m dx, (3) and the factorization operators satisfy the commutation rule ( [A,B] = a+ ] )[mm 3(m ) β. (4) m 3 4 m The PDM problem can be simplified using the map ψ e g ϕ and the change of the independent variable x, ruled as follows: x y = s(x), y x = s (y). (5) To define s as a bijection we shall assume that J (D a ) 0, with D a Dom(H a ) and J = s (see e.g. []). The combination of transformations (5) produces y = [ ] / m(x) J(x)dx+y 0, J(x) = e g(x) = (6) m 0 and ψ(x) = J / (x)ϕ (y) (7)

3 where m 0 and y 0 are integration constants in proper units (hereafter we take y 0 = 0). The function ϕ is the representation of ϕ in the y-space: ϕ(x) = ϕ(s (y)) = [ϕ s ](y) ϕ (y) and viceversa ϕ s = ϕ. A further simplification (see [5] for details) leads to the Schrödinger equation H (a) ϕ (y) := [ ( ) ] d m 0 dy +V (y) ϕ (y) = Eϕ (y), (8) with D (a) = Dom(H (a) ) and Sp(H a ) = Sp(H (a) ). Hence two general cases are distinguishable: MDNT (Mass-dependent null terms) The mass-function m(x) is a solution of the non-linear, second order differential equation ( ) 7 mm 4 +a (m ) = 0. (9) A simple analysis shows that the roots of (9) leading to appropriate bijections s have the form m(x;a) = m 0 (x 0 +λx) 4/(3+4a) (0) with a in the set { A = a 0 = /4,a n = n }, n N, () 4n and x 0 a dimensionless real constant while m 0 and λ are real constants expressed in mass and inverse of distance units respectively. Then the functions m (0) (x) m(x;a 0 ) = m 0 (x 0 +λx), s (0)(x) = ln(x 0 +λx), x t 0 = x 0 λ λ, () define an invertible mapping from D a0 [t 0,+ ) to D (a 0) m (n) (x) m(x;a n ) = m 0 (x 0 +λx) 4n (n+), corresponds to a bijection between D an R and D (an) R. R and the pair (n+) s(n) (x) = (x 0 +λx), n N (3) λ(n+) MINT(Mass-independent null terms) Given a properly defined mass-function m(x), the ordering of P and m(x) in () is a priori fixed as a = b = /4. For instance, the regular functions m R (x) = m 0 +(λx), define a bijection connecting D /4 R to D ( /4) (see Section 5) is given by and connects D /4 R to D ( /4) s R(x) = arcsinh(λx), λ R (4) λ R. Another well behaved bijection m e (x) = m 0 e λx, s e (x) = eλx x 0, (5) λ [t 0,+ ). 3

4 3 Position-dependent mass linear oscillators Let the commutator (4) be a constant, namely [A,B] = ω 0, with ω 0 in frequency units. The β-function is easily found to be β = ω x 0 m / dr ( a+ )( ) m +β 4 m 3/ 0 (6) with β 0 an integration constant which will be omitted in the sequel. The identification ǫ = ω 0 /, after introducing (6) in the Riccati equation (3), leads to a very simple form of the potential in the y-representation V(x) = ω 0 [ x m / dr ] = m 0ω 0 [ x ] Jdr = ( ) m0 ω0 y V (y). (7) Thereby the transformation(6)-(7) tunes to the potentials V(x) exhibiting the equidistant energies ω 0 (n+/) of a constant mass quantum oscillator. Since s is a bijection, (7) admits another lecture: given V (y) = m 0ω0 y, the mapping (6)-(7) leads to a potential V(x) = m 0ω0 (s(x)) which is isospectral to the harmonic oscillator for the masses m(x) allowed by the rule s. It is convenient to introduce a dimensionless notation by taking [E] = ω 0 and [L] = /(m0 ω 0 ) as the units of energy and distance respectively. Hence, the potential (7) reads V (y) = y [E], with y a real number such that y = y[l]. From now on we drop the dimensions and use the same symbol for the physical and the dimensionless variables. We shall return to the expressions with units only if necessary. In each case, the notation will be self-consistent. The same holds for the subscript * labelling the representation of functions and operators in the y-space. The introduction of (6) in () cancels the explicit dependence of the factorization operators A and B on the ordering label a. We have A = a + + ( ) d dy lnj, B = a ( d dy lnj ), [A] = [B] = [E]/ (8) with a (a + ) the conventional annihilation (creation) operator of the quantum oscillator in the y-space a + = a := d dy +y, [a,a + ] =, a + a = N, [N,a ± ] = ±a ±. (9) Here N is the Fock s number operator. A most convenient relationship between A, B and a ± is easily calculated to read AJ / = J / a +, BJ / = J / a. (0) Hence, the action of A and B in the ϕ-space can be established as Aψ = ( J / a + ) ϕ, Bψ = ( J / a ) ϕ. () 4

5 The Hamiltonian H a is clearly isospectral to the one-dimensional quantum oscillator H a ψ = (AB +)ψ = A ( ) J / a ϕ+ψ = J / (N +)ϕ [ = J / d +y ]ϕ = J / H (a) ϕ = (n+)ψ, dy () where [H a ] = [H (a) ] = [E]/. 4 Position-dependent mass singular oscillators Some consequences of the commutation relations (9) are that [a,f] = f y, [a +,f] = f y, [N,f] = f yy f y d dy, f y df dy (3) with f a differentiable function of the position. Moreover, since [a,a + ] = 8H(a) we can introduce the operators c + := a + +f, c := a +f (4) to get [ ( y ] [c,c + ] = 8 H (a) + y := 8h )f (a). (5) The straightforward calculation shows that given f(y) = g 0 y, with g 0 a real constant, the operators k 0 = h (a) /4, k ± = c ± /4 close the su(,) algebra [k,k + ] = k 0, [k 0,k ± ] = ±k ±. (6) The operator h (a) in (5) is the Hamiltonian of the singular oscillator: h (a) = H (a) + g 0 y = d dy +y + g 0 y, [g 0] = [L] 4. (7) Remark that the potential V(y) = y + g 0 admits an infinite point spectrum if g y 0 > / (see e.g. Chs. III.8 and V.35 of []). Moreover, the presence of the centrifugal-like term g 0 constrainsthedomainofdefinitionofh (a) tobed (a) y = [0, + ), and the wave-functions of h (a) are necessarily equal to zero at the origin. In this way, if g 0 = 0, only the odd linear-oscillator functions are recovered (see Section 4.). On the other hand, from (0) and () one obtains the expressions for the generators in the PDM case: C + = A g 0 s (x), C = B g 0 s (x), h a = H a + g 0 s (x). (8) Notice that the relationships C ± J / = J / c ±, h a J / = J / h (a), (9) imply that the operators (8) also close the su(,) algebra provided that K 0 = h a /4 and K ± = C ± /4. 5

6 4. Physical Solutions Let us introduce the mappings ϕ y l e y / u(y) and y z so that the eigenvalue problem h (a) ϕ = Eϕ is mapped to the Kummer equation [3]: zu zz + (l+ ) z u z 4 (l+ E)u = 0, l l g 0 = 0. (30) ( ) Then, for each l ± = α ± = ± +g0, we have a general expression of the form ϕ ± = λ () ± y α ± e y / F (a,c,y )+λ () ± y α ± e y / F (ã, c,y ), (3) with 4a = α ± + E, c = α ± +, 4ã = 3 α ± E and c = 3 α ±. The straightforward calculation shows that both of the general expressions (3) lead to the same physical solution if < g 0. In such a case we omit ϕ and take λ () + = 0, a = n and α + = α to get ( ϕ n (y) = c n y α e y / F n,α+,y) ( = n! Γ(α+n+/) ) /y α e y / L (α /) n (y ), n = 0,,,... (3) where L (γ) n (x) are the Generalized Laguerre Polynomials [4]. The set of energies is then defined by E n = 4n++ +g 0 = 4n+α+ = 4(κ+n), n = 0,,,... (33) The change y y/ in h (a) makes clear that the set (33) corresponds to the quantum oscillator spectrum E n = n+, shifted by α /. This last case has been discussed in [5] for g 0 = /. On the other hand, if g 0 = 0 then α = and we have ϕ n (y) = π 4 (n+)! e y / Hen+( y), n = 0,,,... (34) with Hen+(x) the odd Hermite Polynomials [4] and E n = 4n+3. That is, the wavefunctions (3) are reduced to the odd-oscillator ones as g Applications Each of the pairs (V(x),m(x)), with V(x) = s (x)+g 0 /s (x) describe a PDM quantum system sharing its spectrum with a particle of mass m 0 subject to the singular oscillator interaction. We give some examples below. MDNT For a = a 0 A, we have V (x) = ln (x 0 +λx) g 0 λ + λ ln (x 0 +λx)), (35) 6

7 with () the mapping from D a0 to D (a 0) = [0,+ ) and viceversa. On the other hand, if a = a n A, n N, the potential is given by ( n+ V (x) = λ ) (x 0 +λx) /(n+) + g 0 with D an = [t 0,+ ) and D (an) = [0,+ ) connected by (3). ( ) λ (x 0 +λx) /(n+) (36) n+ MINT If we fix a = /4, the transformation (4) connects D /4 = [0,+ ) with = [0,+ ) and leads to the potential D ( /4) V 3 (x) = arcsinh (λx) g 0 λ + λ arcsinh (λx). (37) In turn, for x 0 =, transformation (5) defines the mapping from D /4 = [0,+ ) to D ( /4) = [0,+ ) and the potential ( ) ( ) ( ) V 4 (x) = 4eλx λx sinh g0 λ λx + e λx csch, x > 0. (38) λ 8 In Fig. the global behavior of the above defined potentials is shown for specific values of the parameters. Fig. shows the eigenfunctions of potentials V (x) and V 4 (x). In the former case the functions are expanded towards infinity by preserving the shape of the constant-mass ones at short distances. In the second case they are squeezed into the vicinity of the origin of coordinates. As expected, no change in the normalization is found. Details can be appreciated in Fig. 3 where one of the squeezed probability densities of V 4 (x) is contrasted with its constant-mass equivalent. As a final remark, the above results make clear that κ = ( + α) determines the representation of SU(,) we are dealing with (see e.g. [5] and [6]). E3 E E V V V3 En E3 E E Sing V5 En E0 E Figure : The PDM potentials (35) (38), together with the singular oscillator of constant mass (Sing). Potential V (x) has been depicted with n =. The first four energy levels (En) are included as reference. In all cases g 0 =, λ = and x 0 is fixed to give D a = [0,+ ). Vertical and horizontal axis are respectively in ω 0 / and dimensionless units. 5. New SU(, ) coherent states Similar to the case of c +, c and h (a), the operators C ± intertwine the Hamiltonian h a with itself, shifted by 4 units of the energy. In other words, C ± work as ladder operators 7

8 n=0 n= n= Figure : The first three wave-functions of the PDM potentials V 4 (x) left and V (x) right for the same values of the parameters as in Figure. Vertical and horizontal axis are in dimensionless units Figure 3: The second excited state probability density of the PDM potential (38) and the singular oscillator of constant mass (dashed) for the parameters indicated in Figure. Vertical and horizontal axis are in dimensionless units. when acting on the wave-functions of h a. We get C ψ 0 = 0 and C ± ψ n = γ ± (n)ψ n±, γ ± (n) = E n± E n +3 g 0. (39) It is convenient to rewrite the coefficients γ ± as follows γ + (n) = 4 (n+)(n+κ), γ (n) = 4 n(n +κ). (40) The su(,) CS for the PDM singular oscillator are now constructed as solutions of the equation C Φ z = zφ z, z C. We get ( ) κ / [ ] / z J(x) (z/4) l Φ z (x) = ϕ l (s(x)) (4) 4 I κ ( z /) l!γ(l+κ) where I ν (z) is the ν-order modified Bessel function of the first kind [4]. The constant mass case is recovered by taking m(x) = m 0, then J = and s(x) = x. A final change z 4z leads from (4) to the well known generalized coherent states of Barut and Girardello [3]. In a similar form we get the involved Perelomov SU(, ) coherent states: [ ] / Γ(l+κ) Φ z (x) = J / (x)( 4z ) κ (4z) l ϕ l (s(x)) (4) l!γ(κ) l=0 with the same recipe to recover the classical results [4]. In the description of Perelomov, each SU(,) CS is connected with a point in the coset space SU(,)/U(). 8 l=0

9 Thereby, the construction (4) corresponds to the applying of the unitary operators Ω(ξ) SU(,)/U() to the lowest state ψ 0. Here 4z = (ξ/ ξ )tanh ξ so that 4z <. The details can be consulted in [4] (see also [6]). Both of the above derived sets of CS preserve the form of the constant-mass case. Important properties like the resolution of the identity are found to be similar to the conventional case. These states also evolve in time without dispersion because their energy eigenvalues are equally spaced. Moreover, the quadratures X(X,P) = (K + +K )/ and P(X,P) = i(k + K )/ appear to satisfy the commutator [X,P] = ik 0. In this sense, the Barut-Girardello CS can be considered as the quadrature states minimizing the inequality relation X P 8 h a. (43) Some other PDM coherent states minimizing similar inequality relations are constructed in [5 9]. The relevant aspect is that the point transformations analyzed in Section simplify the construction of the PDM coherent states to be practically the same method as in the constant-mass case. 6 Concluding remarks The Darboux transformation of the singular oscillator with a varying frequency and constant mass V(y,t) = ω (t)y + g 0 has been studied in [7]. As a result, it was shown that y the CS belonging to V(y,t) are essentially unchanged, just as the results reported here for the stationary case and masses varying with the position m(x). A combination of these approaches would be applicable in the case of a time-dependent frequency ω(t) and a position-dependent mass m(x). It is then expected a similar result: the Barut-Girardello and Perelomov CS will preserve their global properties after the transformations. In this context, it is important to remark that the β-function defined in (6) obeys the fact that the operators A and B are taken to fulfill the oscillator algebra [A,B] = ω 0. That is, the function in (6) corresponds to a particular solution of the Riccati equation (3). As it is well known, general solutions give rise to different algebras (see e.g. [8] and [9]). These algebras have been applied in the construction of a new kind of CS connected with the linear oscillator and its Susy-partners [0, ]. Quite recently, it has been shown that non-linear Susy algebras can be linearized to exhibit the Heisenberg- Weyl structure. In particular, the SU(, ) algebra as connected with the infinite well was analyzed []. Then, the higher order Susy transformations can be also studied for the position-dependent mass systems we have presented in this work. The classical models of the harmonic oscillator and the Pöschl-Teller potentials have been also useful in the solving of mass-dependent systems [] and in the construction of CS [3]. Of particular interest, the study of the Wigner function gives rise to a better understanding of the PDM coherent states [6, 4]. Further insights are in progress. Acknownledgementes The support of CONACyT project F and IPN grant COFAA is acknowledged. 9

10 References [] Schrödinger E., Naturwissenschaften (96) [] Glauber R.J., Chem. Phys. Chem (006) [3] Barut A.O. and Girardello L., Commun. Math. Phys. 4 (97) [4] Perelomov A.M., Commun. Math. Phys. 6 (97); Perelomov A.M., Usp. Fiz. Nauk 3 3 (977); Perelomov A., Generalized Coherent States and Their Applications (Springer-Verlag, Heidelberg, 986) [5] Cruz y Cruz S. and Rosas-Ortiz O., J. Phys. A: Math. Theor (009) [6] Biswas A. and Roy B., Mod. Phys. Lett. A (009) [7] Guo-Xing J., Chang-Ying C., and Zhong-Zhou, R., Commun. Theor. Phys (009) [8] Rubi V.C., and Senthilvelan M., J. Math. Phys (00) [9] Midya B., Roy B. and Biswas A., Phys. Scr (009) [0] Roy B., Europhys. Lett. 7 (005) [] Choquet-Bruhat Y., DeWitt-Morette C. and Dillard-Bleick M., Analysis, Manifolds and Physics (North-Holland, Amstersam, 977) [] Landau L.D. and Lifshitz E.M., Quantum Mechanics (Pergamon, London, 960) [3] Negro J., Nieto L.M. and Rosas-Ortiz O., J. Math. Phys (000); Rosas-Ortiz O., Negro J. and Nieto L.M., Rev. Mex. Fis. 49 S 88 (003) [4] Abramowitz M. and Stegun I. (Eds), Handbook of Mathematical Functions (Dover Pub. Inc., New York, 97) [5] Barut A.O. and Raczka R., Theory of Group Representations And Applications (Polish Scientific Publishers, Poland, 977) [6] Brif C. Vourdas A. and Mann A., J. Phys. A: Math. Gen (996) [7] Samsonov B.F., J. Phys. A: Math. Gen (004) [8] Mielnik B. and Rosas-Ortiz O., J Phys A: Math Gen (004) [9] Andrianov A.A. and Cannata F., J Phys A: Math Gen (004) 0

11 [0] Fernández D., Hussin V. and Nieto L.M., J. Phys. A: Math. Gen (994); Fernández D., Nieto L.M. and Rosas-Ortiz O., J. Phys. A: Math. Gen (995); Rosas-Ortiz O., J. Phys. A: Math. Gen (996); Fernández D. and Hussin V., J. Phys. A: Math. Gen (999) [] Fernández D., Hussin V. and Rosas-Ortiz O., J. Phys. A: Math. Gen (007); Fernández D., Rosas-Ortiz O. and Hussin V., J. Phys. Conf. Ser (008) [] Cruz y Cruz S., Negro J. and Nieto L.M., Phys. Lett. A (007); Cruz y Cruz S., Negro J. and Nieto L.M., J. Phys.: Conf. Ser (008) [3] Cruz y Cruz S., Kuru S. and Negro J., Phys. Lett. A (008) [4] de Souza Dutra A and de Oliveira J.A., Phys. Scr (008)