Abstracts
We analyze the quantum-mechanical behavior of a system described by a one-dimensional asymmetric potential constituted by a step plus (i) a linear barrier or (ii) an exponential barrier. We solve the energy eigenvalue equation by means of the integral representation method, classifying the independent solutions as equivalence classes of homotopic paths in the complex plane. We discuss the structure of the bound states as function of the height of the step and we study the propagation of a sharp-peaked wave packet reflected by the barrier. For both the linear and the exponential barrier we provide an explicit formula for the delay time as a function of the peak energy E. We display the resonant behavior of at energies close to . By analyzing the asymptotic behavior for large energies of the eigenfunctions of the continuous spectrum we also show that, as expected, approaches the classical value for , thus diverging for the step-linear case and vanishing for the step-exponential one.
Keywords:
integral transforms; special functions; solutions of wave equations: bound states; scattering theory
Analisamos o comportamento quântico de um sistema descrito por um potencial assimétrico unidimensional constituído por um degrau junto a (i) uma barreira linear ou (ii) uma barreira exponencial. Resolvemos a equação dos valores próprios da energia por meio do método da representação integral, classificando as soluções independentes como classes de equivalência de caminhos homotópicos no plano complexo. Discutimos a estrutura dos estados ligados como função da altura do degrau e estudamos a propagação de um pacote de onda pontiagudo refletido pela barreira. Para ambos os casos linear e exponencial fornecemos uma fórmula explícita para para o retardo como função da energia de pico E. Exibimos o comportamento ressonante de para energias próximas de . Ao analisar o comportamento assintótico para grandes energias das funções próprias do espectro contínuo, mostramos também que, como esperado, se aproxima do valor clássico para , divergindo assim para o caso linear e zerando para o caso exponencial.
Palavras-chave:
transformações integraais; funçoes especiais; solução das equeções de onda: estados ligados; teoria do espalhamento
1. Introduction
In Ref. [1[1] L. Rizzi, O.F. Piattella, S.L. Cacciatori and V. Gorini, Am. J. Phys. 78, 842 (2010).], we have analyzed the quantum-mechanical behavior of a system described by a one-dimensional asymmetric potential formed by a step plus a harmonic barrier (the “step-harmonic” potential), by using the integral representation method [2[2] H. Hochstadt, The Functions of Mathematical Physics (Dover Publications, New York, 1976), p. 100-105.]. We investigated the behavior of the discrete energy levels (as a function of the height of the step) and of the delay time of a wave packet coming from infinity and bouncing back on the harmonic barrier, as a function of the packet's peak energy and of the height of the step.
Among the convex or concave locally bounded symmetric and confining potentials the harmonic oscillator is a threshold, in that it gives rise to classical isochronous oscillations and evenly spaced quantum energy levels [3[3] J.F. Cariñena, A.M. Perelomov and M.F. Rañada, J. Phys. Conf. Ser. 87, 012007 (2007).,4[4] M. Asorey, J.F. Cariñena, G. Marmo and A. Perelomov, Annals Phys. 322, 1444 (2007).].
In our quantum mechanical step variant of the problem we recover both these features in the limit in which , and the potential reduces to the half-space harmonic oscillator. Then, it is conceivable that the harmonic one is the only confining barrier which displays a constant nonvanishing delay in the limit of high energies. For steeper barriers we expect to vanish at high energies, while for milder ones we expect the delay to become infinite in this limit, in accordance with the corresponding classical situations. Similarly, we expect that, as , the spacing between two neighboring discrete levels tends to infinity in the former case and to zero in the latter.
Here we analyze the “step-linear” and the “step-exponential” potentials. Both these problems can be solved exactly using the integral representation method, which is interesting per se, as it can be applied to more general problems.
2. The step-linear potential
Let M, be positive parameters, and consider the “step-linear” potential
If E denotes the energy of the particle and m its mass, the time-independent Schrödinger equation is
This, for , can be rewritten as
It is convenient to define
which allows us to recast Eq. (3) as
which is the Airy equation (see Refs. [5[5] M. Abramowitz and I.A. Stegun, Handbook of Mathematical Functions (Dover Publications, New York, 1972), 10th. ed.,6[6] O. Vallée and M. Soares, Airy Functions and Applications to Physics (Imperial College Press, London, 2004).]). The general solution of Eq. (5) is
where C and D are arbitrary integration constants, and the two linearly independent solutions and are expressed in section 5.1 in terms of the integral representation method.
Since for , in order for Eq. (6) to be an eigenfunction, we must set D = 0. Hence, we obtain
For , the solution of Eq. (2) has the following from
where A, B and F are arbitrary integration constants and
2.1. The case : Bound states and level spacing
If we obtain
The requirement of continuity of u(x) and of its first derivative in x = 0 is expressed by
System (11) has a non trivial solution iff
where . The energy levels are determined graphically by the intersections of the curves at the two sides of Eq. (12). An example is depicted in Fig. 1 for .
Graphical solutions of Eq. (12) for the energy levels. The red solid line represents the graph of the function while the blue dash-dotted line represents the graph of . In this plot, we have chosen .
In the limit the step of Eq. (1) becomes an infinite barrier. In this case, the energy levels correspond to the zeros of , the denominator of Eq. (12). As expected, these energy levels are the ones of the symmetric confining potential corresponding to the odd eigenfunctions of the latter (Appendix A A. The confining symmetric linear potential Consider the confining symmetric potential U(x)=M|x|. The eigenfunctions can be written as (79) u ( x ) = C 1 Ai ( - α x - β ) x < 0 , u ( x ) = C 2 Ai ( α x - β ) x > 0 , where C1 and C2 are constants fixed by the junction conditions in x = 0. If Ai(-β)≠0, then the continuity of u(x) in x = 0 implies C1=C2. Moreover, the continuity of the derivative implies Ai′(-β)=0. This condition determines the even eigenfunctions and their eigenvalues. If Ai(-β)=0, then the continuity of the derivative implies C1 = - C2. This condition determines the odd eigenfunctions and their eigenvalues. ).
To study the level separation for large energies, consider the asymptotics of for large values of (see Ref. [7[7] B.R. Fabijonas and F.W.J. Olver, SIAM Review 41, 762 (1999).])
where and the coefficients of the series expansions are given by
The inversion of the asymptotic expansion (13) allows us to find, for large values of , the following approximate solution of the equation (see Ref. [7[7] B.R. Fabijonas and F.W.J. Olver, SIAM Review 41, 762 (1999).])
Thus, at the leading order of Eq. (15) the approximate zeros have the form
This is an excellent approximation to the zeros of (see Table 1).
The first three zeros of the Airy function compared with the corresponding approximate zeros from Eq. (16).
For we get,
The spacing behavior n-1/3 is the threshold between concave and convex potentials.
2.2. The case : Scattering and delay
In this case, the (improper) eigenfunctions have the form
The junction conditions in x = 0 are
Solving for the constants, the normalized (with respect to k) improper eigenfunctions are given by
where
As expected, the continuous part of the spectrum () is simple. Note that
From Eq. (20) a generic wave packet
has the form
Then, writing , and take the following form
where
If c(k) is sufficiently regular and non-vanishing only in a small neighborhood of some , then and represent wave packets which move according to the following equations of motion [8[8] P. Caldirola, R. Cirelli and G.M. Prosperi, Introduction to Theoretical Physics (UTET, Milano, 1982) 1st. ed. (in Italian).,9[9] D.J. Griffiths, Introduction to Quantum Mechanics (Benjamin Cummings, San Francisco, 2004), 2nd. ed.]
for the “incoming” wave packet, and
for the reflected “outgoing” one.
The solution thus built represents a particle of well defined momentum which approaches the origin from the right, interacts with the linear potential (at t = ), and is totally reflected. Note that the argument of the complex valued function c(k) determines . The phase shift results in a delay in the rebound, caused by the interaction with the confining linear barrier. The delay is calculated with respect to the case of instantaneous reflection, which takes place in presence of an infinite barrier and for which . From Eqs. (4) and (9) it follows that
where . We compute from Eq. (22). Using the Airy equation , we obtain
where we have suppressed the tilde on the packet peak energy .
We are interested in the behavior of the interaction time for large values of , i.e. for incoming packets with high energy. To this end, we need the asymptotic expansion of for (see Refs. [5[5] M. Abramowitz and I.A. Stegun, Handbook of Mathematical Functions (Dover Publications, New York, 1972), 10th. ed.,7[7] B.R. Fabijonas and F.W.J. Olver, SIAM Review 41, 762 (1999).])
where and the coefficients are
with given in Eq. (14). Dividing the two asymptotic expansions of and we obtain to leading order in
Thus, Eq. (31) becomes
Hence, reintroducing the physical variables, the high-energy behavior of the interaction time is
which is exactly the time a classical particle arriving from infinity with energy E would spend in the region.
In Fig. 2 we plot for a choice of different values of . Note the resonances located at the points (), zeros of , corresponding to the formation of metastable states at the respective energies . The values are the energies of the excited states of the confining potential (see Appendix A A. The confining symmetric linear potential Consider the confining symmetric potential U(x)=M|x|. The eigenfunctions can be written as (79) u ( x ) = C 1 Ai ( - α x - β ) x < 0 , u ( x ) = C 2 Ai ( α x - β ) x > 0 , where C1 and C2 are constants fixed by the junction conditions in x = 0. If Ai(-β)≠0, then the continuity of u(x) in x = 0 implies C1=C2. Moreover, the continuity of the derivative implies Ai′(-β)=0. This condition determines the even eigenfunctions and their eigenvalues. If Ai(-β)=0, then the continuity of the derivative implies C1 = - C2. This condition determines the odd eigenfunctions and their eigenvalues. ), corresponding to even eigenfunctions. The resonances have lifetimes which decrease as the corresponding energies increase and move farther away from the threshold energy . Conversely, as increases, the lifetime of the resonance closest to the height of the step becomes progressively longer and then infinite when the resonance turns into the next bound state. This behavior is evident in Fig. 2, in which the first three plots correspond to values of for which there is only one bound state. In the successive three plots the resonance at has disappeared, having turned into the second bound state.
Plots of the delay (in units of ) versus the energy of the incoming wave packet (red solid lines) for . The parabolic blue dotted lines represent the classical delay, given in Eq. (36).
Comparing Fig. 2 with Fig. 6 of Ref. [1[1] L. Rizzi, O.F. Piattella, S.L. Cacciatori and V. Gorini, Am. J. Phys. 78, 842 (2010).] we note that, whereas in the step-harmonic case the graph of oscillates with decreasing amplitude about the straight line (the half period of the oscillator), in the step-linear case the corresponding graph similarly oscillates about the parabolic line , corresponding to the delay of the classical particle. Furthermore, whereas in the step-harmonic case the resonances are evenly spaced, in the step-linear case their spacing decreases with the energy, corresponding to the behavior as a function of the energy of the eigenvalues of the corresponding (symmetric) potentials and .
3. The step-exponential potential
Let , and be positive parameters, and consider the “step-exponential” potential
For , introduce the following dimensionless quantities
in terms of which the time-independent Schrödinger equation writes as
Setting , Eq. (39) can be cast in the form of a modified Bessel equation (see Ref. [5[5] M. Abramowitz and I.A. Stegun, Handbook of Mathematical Functions (Dover Publications, New York, 1972), 10th. ed.,10[10] G.N. Watson, A Treatise on the Theory of Bessel Functions (Cambridge University Press, Cambridge, 1966), 2nd. ed.])
whose general solution is
where C and D are arbitrary integration constants and and are the modified Bessel functions of order .
The function diverges exponentially for [5[5] M. Abramowitz and I.A. Stegun, Handbook of Mathematical Functions (Dover Publications, New York, 1972), 10th. ed.]. For this reason, in order for u(x) to be a proper (or improper) eigenfunction, we must set D = 0. Therefore Eq. (41) reduces to
Also the solutions of Eq. (40) can be studied with the integral representation method (see section 5.2.).
3.1. The case : Bound states
If we obtain
The junction conditions in give
which, setting can be recast in the form
Graphical solutions of Eq. (45) are shown in Fig. 3 .
Graphical solutions of (45). The red solid line is the graph of , whereas the blue dash-dotted line represents for the choice and .
Plots of the delay (in units of ) versus the energy of the incoming wave packet (solid red lines) for six values of the step height: (upper panels) (lower panels). The blue dotted lines represent the classical delay.
Analogously to what happens in the step-linear case (12), in the limit of an infinite barrier () the energy levels are specified by the zeros of , the denominator of Eq. (44), as a function of and they are the ones of the symmetric confining potential , corresponding to the odd eigenfunctions of the latter (Appendix B B. The confining symmetric exponential potential Consider the confining symmetric potential U(x)=κe|x|∕σ-1. The eigenfunctions can be written as (80) u ( x ) = C 1 K i β ( α e - x ∕ 2 σ ) x < 0 , u ( x ) = C 2 K i β ( α e x ∕ 2 σ ) x > 0 , where C1 and C2 are constants fixed by the junction conditions in x = 0. If Kiβ(α)≠0, then the continuity of u(x) in x = 0 implies C1 = C2. Moreover, the continuity of the derivative implies Kiβ′(α)=0. This condition determines the even eigenfunctions and their eigenvalues. If Kiβ(α)=0, then the continuity of the derivative implies C1 = - C2. This condition determines the odd eigenfunctions and their eigenvalues. ).
To study how the energy levels behave for large energies, we employ the following formula for the asymptotic behavior of the function for large (see Ref. [5[5] M. Abramowitz and I.A. Stegun, Handbook of Mathematical Functions (Dover Publications, New York, 1972), 10th. ed.])
for . Note that expansion (46) can be proved starting from Eq. (76). Therefore, the zeros of , as a function of , are asymptotically the solutions of the following equation
Solving for we obtain
where W (x) is the Lambert function [5[5] M. Abramowitz and I.A. Stegun, Handbook of Mathematical Functions (Dover Publications, New York, 1972), 10th. ed.]. Since for we have for large n that
We see from Eq. (49) that the potential behaves for large x as an infinite square well whose width, up to inessential factors, grows as , an intuitive fact. Moreover,
for , proving thus that the level spacing diverges.
3.2. The case : Scattering and delay
The unbound eigenstates have the form
and the junction conditions are
Therefore, the normalized (with respect to k) improper eigenfunctions are given by
where and
Hence,
Then, following the same argument adopted for the step-linear case, we obtain the following formula for the delay of the rebound of an incoming wavepacket with peak energy
Using Eq. (46), we obtain for large values of
from which
A comment is here in order. In general, taking the derivative of an asymptotic expansion with respect to the variable or a parameter may lead to wrong results. However, in our case this procedure can be justified using the integral representation of Eq. (77) (we leave this as an exercise for the interested reader).
Thus, plugging Eq. (58) into Eq. (56) we obtain the asymptotic behavior of the delay time for large 's, namely
or, in terms of the energy of the particle
As expected, Eq. (60) coincides with the large energy value of the half period of the classical particle subjected to the confining potential .
4. Conclusions
Regarding the structure of the discrete energy spectrum as a function of the height of the barrier, in the two potentials treated in this paper, the same considerations apply as those of the concluding section of Ref. [1[1] L. Rizzi, O.F. Piattella, S.L. Cacciatori and V. Gorini, Am. J. Phys. 78, 842 (2010).]. The only difference is that the energy levels , , of the harmonic oscillator have to be replaced here by the corresponding levels of the confining linear and exponential potentials, respectively (see Appendices A A. The confining symmetric linear potential Consider the confining symmetric potential U(x)=M|x|. The eigenfunctions can be written as (79) u ( x ) = C 1 Ai ( - α x - β ) x < 0 , u ( x ) = C 2 Ai ( α x - β ) x > 0 , where C1 and C2 are constants fixed by the junction conditions in x = 0. If Ai(-β)≠0, then the continuity of u(x) in x = 0 implies C1=C2. Moreover, the continuity of the derivative implies Ai′(-β)=0. This condition determines the even eigenfunctions and their eigenvalues. If Ai(-β)=0, then the continuity of the derivative implies C1 = - C2. This condition determines the odd eigenfunctions and their eigenvalues. and B B. The confining symmetric exponential potential Consider the confining symmetric potential U(x)=κe|x|∕σ-1. The eigenfunctions can be written as (80) u ( x ) = C 1 K i β ( α e - x ∕ 2 σ ) x < 0 , u ( x ) = C 2 K i β ( α e x ∕ 2 σ ) x > 0 , where C1 and C2 are constants fixed by the junction conditions in x = 0. If Kiβ(α)≠0, then the continuity of u(x) in x = 0 implies C1 = C2. Moreover, the continuity of the derivative implies Kiβ′(α)=0. This condition determines the even eigenfunctions and their eigenvalues. If Kiβ(α)=0, then the continuity of the derivative implies C1 = - C2. This condition determines the odd eigenfunctions and their eigenvalues. ). In the case of the step-linear potential, the level spacing Eq. (17) goes to zero as the energy increases, while in the case of the step-exponential one (50) it approaches infinity. As regards the continuous spectrum, we provide in both cases exact expressions for the delay of a wavepacket reflected from the barrier, as a function of the peak packet energy (see Eqs. (31) and (56)). As expected, in both cases these delays exhibit a series of resonances for energies not much larger than , while for large energies, they approach the classical values. The step-harmonic potential is a threshold separating the potential barriers for which the delay time goes to infinity at large energies from those for which it vanishes.
An entirely similar discussion can be applied to the step variant
of any symmetric potential V (x) (V (x) = V (-x)) such that . Indeed the energy eigenvalue equation for V (x) has two linearly independent solutions (x) and (x) the first of which approaches zero very rapidly as whereas the second one diverges steadily without oscillating, and two linearly independent solutions (x) and (x) having a corresponding behavior for (see Refs. [8[8] P. Caldirola, R. Cirelli and G.M. Prosperi, Introduction to Theoretical Physics (UTET, Milano, 1982) 1st. ed. (in Italian).,11[11] F.G. Tricomi, Differential Equations (Einaudi, Torino, 1953), 2d ed. (in Italian).]). Since , the energy eigenvalues are the roots of the equation b(E) = 0. Since the potential is symmetric, these roots correspond to even and odd eigenfunctions alternatively, the ground state being even. However, in the general case the eigenvectors cannot be found explicitly. Therefore, for example, no explicit formula is available in general for the delay time of the reflected packet in the corresponding step variant potential (61).
5. Airy and modified Bessel functions through the integral representations
In this section we solve the energy eigenvalue equations by means of the integral representation method, classifying the independent solutions as equivalence classes of homotopic paths in the complex plane. For the step-linear case we obtain Airy function, while for the step-exponential case we get modified Bessel functions. This technique is interesting per se, as it can be applied to more general cases, provided one is able to guess the correct integral kernel. The Airy case is somehow classical, while the Bessel case is more interesting. We present them both for completeness.
5.1. The step-linear case: Airy functions
We look for a solution of Eq. (5) of the form
where is a path in the complex plane and f is an holomorphic function. Plugging Eq. (62) into Eq. (5) we find
Integrating by parts, we obtain:
Therefore, E(y) is a solution of Eq. (5) if
Hence, a class of solutions of the Airy equation is of the form
where is a suitable path for which the contour term vanishes.
The integrand of Eq. (66) entire. Thus, by Cauchy theorem, every closed path represents the trivial solution E(y) = 0.
Consider an unbounded path. In order for to vanish, we require the leading term in the exponent of f (t) (i.e. t3) to have a negative real part. Therefore, the acceptable unbounded paths are those whose phase is confined to the regions (n = 0, 1, 2). These possible paths are showed in Fig. 5, where the allowed sectors (n = 0, 1, 2) are shaded.
Paths with both endpoints in the same sector (e.g. in Fig. 5) can be closed at infinity using Jordan's Lemma; therefore, they correspond to the trivial solution. The only non-trivial paths are those which link different sectors. There are only 3 non-equivalent classes of such paths which we dub , and respectively (see Fig. 5).
Taking into account Cauchy theorem, these paths satisfy the relation in the sense that the corresponding solutions are not independent. The conventional Airy functions and are the independent solutions of such that (see Ref. [5[5] M. Abramowitz and I.A. Stegun, Handbook of Mathematical Functions (Dover Publications, New York, 1972), 10th. ed.])
and
Denoting by the solutions in Eq. (66) corresponding to the paths (i = 1, 2), it is not difficult to show that
We leave the details to the interested reader (hint: Check the above expressions and their first derivatives in Eq. (66) for . In this case the integrals correspond to Euler Gamma functions).
5.2. The step-exponential case: Modified Bessel functions
Consider the modified Bessel equation
with . A convenient kernel for the integral representation is the following:
We look for solutions of the form
Plugging Eq. (71) into Eq. (69), we get
Integrating by parts, Eq. (72) gives
Then, the integral on the right hand side of Eq. (71) is a solution of Eq. (69) if
Up to a normalization, the solution is . If is integer then f (t) is either entire () or meromorphic (), otherwise it has infinite branch points located at (), see Fig. 6. In the latter instance, the usual procedure is to define a domain in which the function is holomorphic by cutting the t-plane and thus forbidding loops around the branch points. In Fig. 6 a convenient choice for the cuts is also shown.
Cut of the complex t-plane. The dashed red and the solid blue thick paths represent the only classes of paths which provide independent solutions to the modified Bessel equation.
Recalling that , the contour condition is
There are 4 different classes of paths for which Eq. (75) is satisfied and
is well defined. The “paths zoology” is more complicated and rich than in the linear and in the harmonic case [1[1] L. Rizzi, O.F. Piattella, S.L. Cacciatori and V. Gorini, Am. J. Phys. 78, 842 (2010).].
-
Closed paths. For any closed path, the contour condition is trivially satisfied and any integral along a path enclosing a region where the integrand function is holomorphic (i.e. the path does not cross the cuts) vanishes. An example is shown by the thin dashed black line in Fig. 6.
-
Infinite paths. For paths whose endpoints are both at infinity, the function diverges or oscillates. On the other hand, the exponential vanishes for (recall that ). Since then must stretch at infinity in one of the sectors defined by (), which are represented by the shaded regions in Fig. 6. Incidentally, in these bands, when there are no cuts, one can “close” the paths at infinity, by virtue of Jordan's Lemma. Examples of this class of paths are the black solid thin lines of Fig. 6.
-
Semi-infinite paths. By “semi-infinite” paths we mean paths starting from a point, say , and ending at infinity. These paths must go to infinity in the shaded bands . For the starting point , the contour condition demands that . This means that these paths must start from one of the points , which are the zeroes of the hyperbolic sine function. It is easy to prove that the integral of Eq. (76) performed along any two such paths lying in the same band gives the same result (indeed, recall that is periodic). Two examples of this class of paths are the blue solid thick lines in Fig. 6.
-
Finite paths. These are the paths starting and ending in two points say and , with . The contour condition can be satisfied in two different ways: either the values of the contour part are equal at the endpoints, or the contour part vanishes at the endpoints. The former case accounts for paths which do not cross any cut and start from any point , ending at (). Examples of this class of paths are represented by the red dash-dotted thick line in Fig. 6. The latter case is realized by paths connecting the branch points and is represented by the red dashed thick line in Fig. 6.
Taking into account the periodicity of the integrand function (the period is in the domain where it is holomorphic), and the Cauchy theorem, it is easy to show that the integrals along the two kinds of finite paths are proportional to the integral along the “fundamental” path . Moreover, the integrals along any one of the infinite paths are linear combinations of the ones performed along the finite and semi-infinite paths.
In conclusion, two linear independent solutions of Eq. (69) are
and
where and correspond, respectively, to the integrals performed along the solid blue and the dashed red thick lines in Fig. 6 (a semi-infinite path and a finite one). It can be shown that both solutions (derived here for ) can be analytically continued throughout the whole z-plane cut along the negative real axis (see Ref. [2[2] H. Hochstadt, The Functions of Mathematical Physics (Dover Publications, New York, 1976), p. 100-105.]).
A. The confining symmetric linear potential
Consider the confining symmetric potential . The eigenfunctions can be written as
where and are constants fixed by the junction conditions in x = 0. If , then the continuity of u(x) in x = 0 implies . Moreover, the continuity of the derivative implies . This condition determines the even eigenfunctions and their eigenvalues. If , then the continuity of the derivative implies = - . This condition determines the odd eigenfunctions and their eigenvalues.
B. The confining symmetric exponential potential
Consider the confining symmetric potential . The eigenfunctions can be written as
where and are constants fixed by the junction conditions in x = 0. If , then the continuity of u(x) in x = 0 implies = . Moreover, the continuity of the derivative implies . This condition determines the even eigenfunctions and their eigenvalues. If , then the continuity of the derivative implies = - . This condition determines the odd eigenfunctions and their eigenvalues.
References
- [1] L. Rizzi, O.F. Piattella, S.L. Cacciatori and V. Gorini, Am. J. Phys. 78, 842 (2010).
- [2] H. Hochstadt, The Functions of Mathematical Physics (Dover Publications, New York, 1976), p. 100-105.
- [3] J.F. Cariñena, A.M. Perelomov and M.F. Rañada, J. Phys. Conf. Ser. 87, 012007 (2007).
- [4] M. Asorey, J.F. Cariñena, G. Marmo and A. Perelomov, Annals Phys. 322, 1444 (2007).
- [5] M. Abramowitz and I.A. Stegun, Handbook of Mathematical Functions (Dover Publications, New York, 1972), 10th. ed.
- [6] O. Vallée and M. Soares, Airy Functions and Applications to Physics (Imperial College Press, London, 2004).
- [7] B.R. Fabijonas and F.W.J. Olver, SIAM Review 41, 762 (1999).
- [8] P. Caldirola, R. Cirelli and G.M. Prosperi, Introduction to Theoretical Physics (UTET, Milano, 1982) 1st. ed. (in Italian).
- [9] D.J. Griffiths, Introduction to Quantum Mechanics (Benjamin Cummings, San Francisco, 2004), 2nd. ed.
- [10] G.N. Watson, A Treatise on the Theory of Bessel Functions (Cambridge University Press, Cambridge, 1966), 2nd. ed.
- [11] F.G. Tricomi, Differential Equations (Einaudi, Torino, 1953), 2d ed. (in Italian).
Publication Dates
-
Publication in this collection
June 2016
History
-
Received
12 Sept 2015 -
Accepted
24 Oct 2015