Nuclear Physics
Step potential and Ramsauer-Townsend effect in Wigner-Dunkl quantum mechanics
A. Askaria
H. Hassanabadia
W. S. Chungb
a Faculty of physics, Shahrood University of Technology, Shahrood, Iran. E-mail: aliaskari@shahroodut.ac.ir; h.hasanabadi@shahroodut.ac.ir.
b Department of Physics and Research Institute of Natural Science, College of Natural Science, Gyeongsang National University, Jinju 660-701, Korea. E-mail:; mimip44@naver.com
Abstract
In this paper the continuity equation for Wigner-Dunkl-Schrödinger equation is studied. Some properties of v-deformed functions related to Dunkl derivative are also studied. Based on these, the step potential and Ramsauer-Townsend effect are discussed in Wigner-Dunkl quantum mechanics.
Keywords: Wigner-Dunkl-Schrödinger equation; Ramsauer-Townsend effect
1. Introduction
In 1950 Wigner [1] proposed a deformed Heisenberg algebra including reflection operator in the form,
x^,p^=i(1+2νP),
(1)
where P denotes the reflection operator Pf(x) = f(-x). The realization of this algebra [2,3] was given by
p^=1iDx, x^=x,
(2)
where we set ℏ = 1 and the Dunkl derivative [3] is defined as
Dx=∂x+νx(1-P).
(3)
In fact, Wigner was the first who discussed the question on supreme level of the Heisenberg-Lie equations on the commutation relation between the momentum and position operator. In 1951 Yang simply developed the problem treated by Wigner and obtained the well-known non-canonical description of the non-relativistic momentum operator [2]. In 1980 N. Mukunda et al investigated Energy position and momentum eigenstates of para-Bose oscillator operators. They found that the two apparently different solutions obtained by Ohnuki and Kamefuchi in this context are actually unitarily equivalent [3]. Also in next year coherent states and the minimum uncertainty states of para-Bose oscillator operators investigate by J. K. Sharma et al. [4]. In 1989 Dunkl constructed a commutative set of first-order differential difference operators associated to the second-order operator [5]. The canonical approach to the non-relativistic quantum mechanics with v = 0 is a special case of the non-canonical approach for arbitrary ν that was proposed by Wigner. Reference [6] shows the Wigner function of the ground state for the para-Bose oscillator (it is an oscillator that was discussed in Refs. [3,4]) that generalizes Gaussian distribution almost overlaps with the Wigner function of the canonical non-relativistic bosonic harmonic oscillator. If there is overlap of the Wigner function of non-canonical ground state with the any excited canonical state, then it means that not only Ramsauer-Townsend effect, but also a lot of other physical effects can be obtained theoretically if to replace canonical momentum operator with non-canonical one. Some studies related to the Wigner-Dunkl quantum mechanics have been accomplished in [7-12].
In one-dimensional Wigner-Dunkl quantum mechanics, the inner product is given by [7,12]
〈f|g〉=∫-∞∞g*(x)f(x)|x|2νdx,
(4)
where |x|
2v
is a weight function. The expectation value of a physical operator 𝒪 with respect to the state ψ(x, t) is defined by
〈O〉=〈ψ|Oψ〉=∫-∞∞ψ*(x,t)Oψ(x,t)|x|2νdx,
(5)
and 𝒪 is a Hermitian operator if it obeys
〈Oψ|ψ〉=〈ψ|Oψ〉.
(6)
For the weight function (6) the momentum operator p^=1/iDx is a Hermitian operator.
In this paper we study the continuity Equation for Wigner-Dunkl-Schrödinger equation. We discuss some properties of v-deformed functions related to Dunkl derivative. Using these we discuss the step potential and Ramsauer-Townsend effect in Wigner-Dunkl quantum mechanics. This paper is organized as follows: In 2 we discuss continuity equation for Wigner-Dunkl-Schrödinger equation. In 3 we discuss the v-deformed functions. In 4 we discuss Step potential. In 5 we discuss Ramsauer-Townsend effect.
2. Continuity equation for Wigner-Dunkl-Schrödinger equation
Now let us consider the time-dependent Wigner-Dunkl-Schrödinger equation,
i∂ψ(x,t)∂t=-12mDx2+V(x)ψ(x,t).
(7)
Now let us derive the continuity equation for the Wigner-Dunkl-Schrödinger equation. From Eq.(7) we have
i∂|ψ(x,t)|2∂t=12mψ(x,t)Dx2ψ*(x,t)-ψ*(x,t)Dx2ψ(x,t).
(8)
The wave function can be split into the even part and odd part,
ψe=12(1+P)ψ=12ψ(x)+ψ(-x),
(9)
ψo=12(1-P)ψ=12ψ(x)-ψ(-x).
(10)
Now let us multiply the weight function K(x)=|x|2ν by Eq.(8),
∂K(x)|ψ(x,t)|2∂t=12miK(x)ψ(x,t)Dx2ψ*(x,t)-ψ*(x,t)Dx2ψ(x,t),
(11)
Let us set
ρ=K(x)|ψ(x,t)|2,
(12)
where we set
ψ=ψe+ψo.
(13)
Then we have
ρ=ρe+ρo,
(14)
where
ρe=K(x)(|ψe|2+|ψo|2),
(15)
ρo=K(x)(ψeψo*+ψoψe*).
(16)
Thus we have the continuity equation of the form,
∂tρ=-∂xJ+f(x),
(17)
where the flux is
J=12mi|x|2ν(ψ∂xψ*-ψ*∂xψ),
(18)
and the source is
f(x)=-νmix2|x|2ν(ψe*ψo-ψeψo*).
(19)
The derivation of Eq.(17) is given in Appendix A. Then we have
ddt∫-∞∞|x|2ν|ψ|2dx=0,
(20)
where we have
∫-∞∞f(x)dx=0,
(21)
because f is odd. Thus, |x|2ν|ψ|2 can be interpreted as the probability density function.
Now let us consider the time-independent Wigner-Dunkl-Schrödinger equation,
-12mDx2+V(x)ψ(x)=Eψ(x),
(22)
or
Dx2ψ(x)=2m(V(x)-E)ψ(x).
(23)
This means that Dx2ψ(x) is not continuous in general because the potential can be discontinuous.
To find the continuity condition in the Wigner-Dunkl-Schrödinger equation, we should find the inverse of Dunkl derivative. The Dunkl derivative can be written as
Dx=1+νx(1-P)∂x-1∂x.
(24)
The inverse of Dunkl derivative is
Dx-1=∂x-11+νx(1-P)∂x-1-1=∂x-1∑n=0∞(-1)nνx(1-P)∂x-1n,
(25)
where ∂x-1 denotes the ordinary integration. Thus we have
Dx-1xN=xN+1N+1ν,
(26)
where Dunkl number is
nν=n+ν(1-(-1)n).
(27)
Because the inverse of Dunkl derivative is expressed in terms of multiple of the ordinary integration, from Eq. (23) we know that both Dxψ and ψ are continuous.
3. The v-deformed functions
Now let us find the v-exponential function obeying
Dxeν(ax)=aeν(ax), eν(0)=1.
(28)
We consider the ν-deformed differential equation
Dxy(x)=ay(x), y(0)=1.
(29)
Let us set
y(x)=ye(x)+yo(x),
(30)
where y
e
(x) is the even function obeying Pye(x)=ye(x) while y
o
(x) is the odd function obeying Pyo(x)=-yo(x). Inserting Eq. (30) into Eq. (28) and splitting into the even part and odd part we get
dye(x)dx=ayo(x),
(31)
dyo(x)dx+2νxyo(x)=aye(x).
(32)
Let us set
ye(x)=∑n=0∞anx2n,
(33)
yo(x)=∑n=0∞bnx2n+1.
(34)
Inserting Eqs. (33,34) into Eq. (31,32), we get
2(n+1)an+1=abn,
(35)
(2n+1+2ν)bn=aan.
(36)
From the above equations we have
an+1=a22(n+1)(2n+1+2ν)an.
(37)
Thus, we have
an=1n!ν+12na22n,
(38)
bn=1n!ν+12n+1a22n+1.
(39)
Thus, we have
y(x)=eν(ax)=coshν(ax)+sinhν(ax),
(40)
where
coshν(ax)=∑n=0∞1n!ν+12nax22n=0F1;ν+12;a2x24,
(41)
sinhν(ax)=∑n=0∞1n!ν+12n+1ax22n+1=ax2ν+10F1;ν+32;a2x24,
(42)
and
0F1(;a;x)=∑n=0∞1n!(a)nxn,
(43)
and
(a)0=1, (a)n=a(a+1)(a+2)⋯(a+n-1).
(44)
One can also express the ν-deformed hyperbolic functions as
coshν(ax)=|ax|2-ν+12×Γν+12Iν-1/2(|ax|),
(45)
sinhν(ax)=ax2|ax|2-ν-12×Γν+12Iν+1/2(|ax|),
(46)
where I
α
(x) denotes the modified Bessel function. These deformed hyperbolic functions reduce to cosh(ax) and sinh(ax) in the limit v → 0. The v-deformed hyperbolic functions obey
P coshν(ax)=coshν(ax),
(47)
P sinhν(ax)=-sinhν(ax).
(48)
Action of the v-derivative gives
Dxeν(ax)=aeν(ax),
(49)
Dxcoshν(ax)=asinhν(ax),
(50)
Dxsinhν(ax)=acoshν(ax).
(51)
If we replace x → ix, , we have the v-deformed Euler relation
eν(iax)=cosν(ax)+isinν(ax),
(52)
where
cosν(ax)=∑n=0∞(-1)nn!ν+12nax22n=0F1;ν+12;-a2x24,
(53)
sinν(ax)=∑n=0∞(-1)nn!ν+12n+1ax22n+1=ax2ν+10F1;ν+32;-a2x24.
(54)
One can also express the v-deformed trigonometric functions as
cosν(ax)=|ax|2-ν+12Γν+12Jν-1/2(|ax|),
(55)
sinν(ax)=ax2|ax|2-ν-12Γν+12Jν+1/2(|ax|),
(56)
where Jα(x) denotes the Bessel function. The v-deformed trigonometric functions obey the following relations
Dxcosν(ax)=-asinν(ax),
(57)
Dxsinν(ax)=acosν(ax).
(58)
Figure 1 shows the plot of y=cosν(x) for v = 0 (Pink), v = 0.1 (Brown) and for v = -0.1 (Gray). Figure 2 shows the plot of y=sinν(x) for v = 0 (Pink), v = 0.1 (Brown) and for v = -0.1 (Gray).
4. Step potential
Now let us consider the step potential problem whose potential is given by
V(x)=0 (x<0)V0 (x>0).
(59)
Now let us consider the case of 0 < V
0 < E. Then, Wigner-Dunkl-Schrödinger equation reads
-12mDx2ψI(x)=EψI(x),
(60)
for x < 0, while it reads
-12mDx2+V0ψII(x)=EψII(x),
(61)
for x > 0. The solution is given by
ψI(x)=eν(ik0x)+reν(-ik0x),
(62)
and
ψII(x)=teν(iqx),
(63)
where
k0=2mE, q=2m(E-V0).
(64)
From the continuity of ψ and Dxψ we have
1+r=t, 1-r=qk0t.
(65)
Thus we have
r=k0-qk0+q, t=2k0k0+q,
(66)
which is the same as the case of v = 0. In this case the transmission and reflection flux are
R=r2, T=qk0t2.
(67)
We see that transmission and reflection flux are independent of the v-deformed parameter. Figure 3 shows the plot of R and T versus E with V
0 = 1.
5. Ramsauer-Townsend effect
Let us consider the quantum well whose potential is given by
V(x)=0 (x<0, Region I)-V0 (0<x<a, Region II)0 x>a, Region III,
(68)
where V
0 is a positive constant and we assume E > 0.
Now let us consider the Wigner-Dunkl-Schrödinger equation for three cases:
-12mDx2ψI=EψI,
(69)
-12mDx2-V0ψII=EψII,
(70)
-12mDx2ψIII=EψIII.
(71)
Solving three equations, we get
ψI=eν(ikx)+Aeν(-ikx),
(72)
ψII=Beν(iqx)+Ceν(-iqx),
(73)
ψIII=Deν(ikx),
(74)
where
k=2mE, q=2m(E+V0).
(75)
Now the boundary conditions are the continuity of the wave functions and their first Dunkl derivatives at the boundaries, which are
A-B-C=-1,
(76)
kA+qB-qC=k,
(77)
Beν(iqa)+Ceν(-iqa)=Deν(ika),
(78)
qBeν(iqa)-qCeν(-iqa)=kDeν(ika).
(79)
Solving Eqs. (76-79) for A we get
A=i(k2-q2)sinν(qa)-2qkcosν(qa)+i(k2+q2)sinν(qa).
(80)
The reflection probability density is given by
R=|A|2eν(-ika)eν(ika)=|A|2(sinν2(ka)+cosν2(ka)),
(81)
where
|A|2=(k2-q2)2sinν2(qa)(k2+q2)2sinν2(qa)+4k2q2cosν2(qa).
(82)
In Fig. 4 shows the behavior of the reflection probability density versus energy. Now let us investigate the Ramsauer-Townsend effect. This effect is a physical phenomenon involving the scattering of low-energy electrons by atoms of a noble gas. Ramsauer-Townsend effect is no reflection condition. From |A|2 = 0, we have
sinν(qa)=0.
(83)
Thus, we have
q=2αν+1/2,pa, p=1,2,⋯,
(84)
where αν+1/2,p denotes p-th zero of Jν+1/2(x).
6. Conclusion
In this paper we derived the continuity equation for Wigner-Dunkl-Schrödinger equation. We found the flux and probability density. We found that the probability conserved in time. From the Dunkl integral ( inverse of Dunkl derivative ), we found that the wave function and first order Dunkl derivative are continuous although the potential is not continuous. We introduced the v-deformed functions related to Dunkl derivative and investigated their mathematical properties. Using the continuity condition in the Wigner-Dunkl-Schrödinger equation we discussed two examples; the step potential problem and Ramsauer-Townsend effect.
References
1. E. Wigner, Do the Equations of Motion Determine the Quantum Mechanical
Commutation Relations? Phys. Rev.
77 (1950) 711.
[ Links ]
2. L. Yang, A Note on the Quantum Rule of the Harmonic Oscillator Phys.
Rev.
84 (1951) 788.
[ Links ]
3. N. Mukunda, E.C.G. Sudarshan, J.K. Sharma and C.L. Mehta, J. Math.
Phys.
21 (1980).
[ Links ]
4. J.K. Sharma, C.L. Mehta, N. Mukunda and E.C.G. Sudarshan, J. Math.
Phys.
22 (1981).
[ Links ]
5. C. Dunkl, Trans. Amer. Math. Soc.
311 (1989) 167.
[ Links ]
6. E.I. Jafarov, S. Lievens and J. Van der Jeugt, J. Phys.
A41 (2008) 235201.
[ Links ]
7. V. Genest, M. Ismail, L. Vinet and A. Zhedanov, J. Phys. A
46 (2013) 145201.
[ Links ]
8. V. Genest, M. Ismail, L. Vinet and A. Zhedanov, The Dunkl Oscillator in the
Plane II: Representations of the Symmetry Algebra Comm. Math.
Phys.
329 (2014) 999.
[ Links ]
9. V. Genest , L. Vinet andA. Zhedanov, J. Phys: Conf. Ser.
512 (2014) 012010.
[ Links ]
10. H. Carrion and R. Rodrigues, Mod. Phys. Lett. A
25 (2010) 2507.
[ Links ]
11. S. Sargolzaeipor, H. Hassanabadi and W. Chung, Mod. Phys.
Lett. A 33 (2018) 1850146.
[ Links ]
12. W. Chung and H. Hassanabadi, Mod. Phys. Lett. A
34 (2019) 1950190.
[ Links ]
Appendix A
Appendix A
In right side of Eq. (11), let us set
I=K(x)(ψ(x,t)Dx2ψ*(x,t)-ψ*(x,t)Dx2ψ(x,t)).
(1)
Then we have
I=Kψ*∂2+2νx∂-νx2(1-P)ψ-Kψ∂2+2νx∂-νx2(1-P)ψ*=K(ψe*+ψo*)×∂2+2νx∂-νx2(1-P)(ψe+ψo)-K(ψe+ψo)∂2+2νx∂-νx2(1-P)(ψe*+ψo*)=Kψe*∂2+2νx∂ψe+Kψo*∂2+2νx∂ψe+Kψe*∂2+2νx∂-2νx2ψo+Kψo*∂2+2νx∂-2νx2ψo-Kψe∂2+2νx∂ψe*-Kψo∂2+2νx∂ψe*-Kψe∂2+2νx∂-2νx2ψo*-Kψo∂2+2νx∂-2νx2ψo*=-∂xK(x)(ψ∂xψ*-ψ*∂xψ)-2νx2K(x)(ψe*ψo-ψeψo*),
(A.2)
where we used
K'(x)=2νxK(x).
(A.3)