Advances in High Energy Physics

Volume 2013 (2013), Article ID 325431, 7 pages

http://dx.doi.org/10.1155/2013/325431

Research Article

Masses and Widths of Resonances for the Cornell Potential

M. N. Sergeenko

Stepanov Institute of Physics, Belarus National Academy of Sciences, 68 Nezavisimosti Avenue, 220072 Minsk, Belarus

Received 5 April 2013; Revised 10 July 2013; Accepted 1 September 2013

Academic Editor: Terry Sloan

Copyright © 2013 M. N. Sergeenko. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Physical properties of the Cornell potential in the complex-mass scheme are investigated. Two exact asymptotic solutions of relativistic wave equation for the coulombic and linear components of the potential are used to derive the resonance complex-mass formula. The centered masses and total widths of the leading -family resonances are calculated.

1. Introduction

Resonance is the tendency of a system to oscillate at a greater amplitude at some frequencies. Most particles listed in the Particle Data Group tables [1] are unstable. A thorough understanding of the physics summarized by the PDG is related to the concept of a resonance.

Operators in Quantum Mechanics (QM) are Hermitian, and the corresponding eigenvalues are real. However, in scattering experiment, the wave function requires different boundary condition that is why the complex energy is required [2, 3]. The concept of a purely outgoing wave belonging to the complex eigenvalue was introduced in 1939 by Siegert [4]. It is an appropriate tool in the studying of resonances. Resonance is present as transient oscillations associated with metastable states. Fundamentals of scattering theory and strict mathematical definition of resonances in QM were considered in [5, 6].

In elastic and inelastic scattering experiments, the resonances are associated with the scattering particle and the target [7]. Resonances can be thought of as unstable particles with the particle's complex energy poles in the scattering amplitude. In such case the energy decreases exponentially with the time, so that the damping constant is a measure of the lifetime, , of the oscillation.

In particle physics resonances arise as unstable intermediate states with complex masses. Resonances in Quantum Field Theory (QFT) are described by the complex-mass poles of the scattering matrix [3]. In potential theory the decay widths come from second- or higher-order perturbation theory. The complex eigenvalues correspond to a first-order pole of the -matrix [8]. The origin of mass and widths of resonances in QFT originate from the real and imaginary parts of Feynman diagrams, that is, from propagators appearing in loops. The masses of the states develop imaginary masses from loop corrections [9, 10].

In the pole approach the parameters of resonances are defined in terms of the pole position in the complex -plane as [5, 6, 11, 12]

where is the complex Mandelstam variable, is the resonance centered mass, and is the resonance width. In this approach, resonances are defined as the poles of the meromorphic continuation of the cut-off resolvent, which are shown to be the same as the poles of the meromorphically continued -matrix. This definition has the advantage of being quite universal regarding the pole position but can only be applied if the amplitude can be analytically continued in a reliable way.

The complex-mass scheme provides a consistent framework for dealing with unstable particles and has been successfully applied to various loop calculations [9]. The rigorous QM definition of a resonance requires determining the pole position in the second Riemann sheet of the analytically continued partial-wave scattering amplitude in the complex Mandelstam variable plane [13]. This definition has the advantage of being quite universal regarding the pole position but can only be applied if the amplitude in scattering experiment can be analytically continued in a reliable way.

However, there are two basic problems in quantum physics: the scattering problem and bound state problem. Resonances are usually studied in scattering experiment, for example, potential scattering. But, resonances are quasibound states in the -channel at that means one can use another approach to consider resonances, that is, in bound state region using the quark potential model.

In perturbative QCD, the essential interaction at small distances is instantaneous coulombic one-gluon exchange (OGE); in QCD, it is , , or Coulomb scattering [14]. Therefore, one expects from OGE a Coulomb-like contribution to the potential: at . For large distances, from lattice-gauge-theory computations [15], it follows that this rise is an approximately linear; that is, , where GeV² is the string tension. These two contributions by simple summation lead to the famous Cornell potential [15–17]:

and its parameters are directly related to basic physical quantities noted above. All phenomenologically acceptable QCD-inspired potentials are only variations around this potential.

The potential (2) is one of the most popular in hadron physics and incorporates in clear form the basic features of the strong interaction. In hadron physics, the nature of the potential is very important. There are normalizable solutions for scalar-like potentials, but not for vector-like. The effective interaction has to be Lorentz-scalar in order to confine quarks and gluons [18]; that is, we take the potential (2) to be Lorentz-scalar.

Our aim is to find in analytic form the energy/mass eigenvalues for the potential (2). This problem is not easy if one uses known relativistic wave equations. However, this aim can be achieved with the use of the semiclassical wave equation [19]. An important feature of this equation is that, for two and more turning-point problems, it can be solved exactly by the conventional WKB method [19]. Another way is to use the complex-mass scheme.

In this work, in contrast to the usual analysis of the scattering amplitude, we consider mesonic resonances in bound state region, that is, quasibound states of the meson constituents. In traditional approach to investigate resonances, one deals with the scattering theory, exploring the properties of -matrix and partial amplitudes. We consider mesonic resonances to be the transient excited states of a quark-antiquark system in bound state region. We use the potential approach and analyze the mass spectrum generated from solution of the relativistic wave equation for two components of the Cornell potential—coulombic and linear.

We use the advantage of analyzing the system in the complex plane (the complex-mass scheme) that has important features such as a simpler and more general framework [9, 10]. We show that two asymptotic components of the Cornell potential, the coulombic term and linear one, yield the complex masses of resonances that allows to calculate in unified way their centered masses and total widths.

2. The Complex Eigenmomenta

Operators of physical values in conventional QM are Hermitian. However, the nature demonstrates its nonexhausting complexity. And now, we must accept the fact that energy and mass are complex. This is clearly shown in particle physics [7, 9, 20].

It is possible to extend the Hamiltonian that describes a quantum system into the complex domain while still retaining the fundamental properties of a quantum theory. One of such approaches is Complex Quantum Mechanics (CQM) [21]. The CQM has proved to be so interesting that research activity in this area has motivated studies of complex classical mechanics [22]. Each unstable particle in CQM is associated with a well-defined object, which is a discrete generalized eigenstate of the Hamiltonian having its real and negative imaginary parts being the binding energy, , and half-width, , of the particle, respectively [21, 22]. Here corresponds to the inverse of its lifetime, , of a resonance state. The time dependence of the survival amplitude in the small-time region is strongly correlated with the asymptotic behavior of the energy spectrum of the system.

Complex energy means the complex momentum and mass. We observe an important hint studying the hydrogen atom problem. The total energy eigenvalues for the nonrelativistic Coulomb problem can be written with the use of complex quantities in the form of kinetic energy for a free particle:

Here is the electron's momentum eigenvalue with the imaginary discrete velocity, , . This means that the motion of the electron in a hydrogen atom is free but restricted by the “walls” of the potential [19, 23].

The Cornell potential (2) is unique in that sense, and it yields the complex masses of resonances. To show that, analyze the eigenvalues obtained separately for two components of the potential (2), that is, the coulombic term and the linear one.

Relativistic two-body Coulomb problem for particles of equal masses can be solved analytically [23]. The exact expression for the c.m. energy squared is well known and can be written in the form of two free relativistic particles as [24, 25]

where is given above and . Here we have introduced the imaginary momentum eigenvalues, .

The linear term of the Cornell potential (2) can be dealt with analogously. In this case the exact solution is also well known [23–25]:

This expression does not contain the mass term and can be written similar to (4) relativistic form:

where is the real momentum eigenvalue.

Thus, two asymptotic additive terms of the potential (2), and , separately, yield the imaginary (4) and the real (6) momentum eigenvalues. These terms of the potential represent two “different physics” (coulombic OGE and linear string tension) and therefore two different realms of the interaction. Each of these two expressions, (4) and (6), is exact and was obtained independently; therefore, we can consider the complex sum:

Thus, we accept the complex momentum eigenvalues (7) that means the total energy and mass should be complex as well.

3. The Complex Quark Mass

The real and imaginary parts of complex numbers/functions are independent. Particle decays are associated with complex mass. The centered mass, , and the total width, , are process-independent parameters of the resonance; they are the real and imaginary parts of the complex mass .

It is an experimental fact that the dependence is linear for light mesons [20]. However, at present, the best way to reproduce the experimental masses of particles is to rescale the entire spectrum given by (5) assuming that the masses of the mesons are expressed by the relation [18]

where is a constant energy (shift parameter). Relation (8) is used to shift the spectra and appears as a means to simulate the effects of unknown structure approximately. But, if we rewrite (8) in the usual relativistic form,

where Re is given by (6), we come to the concept of the imaginary mass, . Here in (9), we have introduced the notation, . What is the mass ?

The required shift of the spectra naturally follows from the asymptotic solution of the semiclassical wave equation for the potential (2) [23–26]. To show that, we account for the “weak coupling effect”; that is, together with the linear dependence in (6), we should include the contribution of the coulombic term, , of the potential.

These kinds of calculations result in the asymptotic expression similar to (8) [23, 24]:

The additional term, , arises from the interference of the coulombic and linear components of the Cornell potential (2). Comparing (10) with (9), we obtain

The interference term in (10) contains only the parameters of the potential (2) and is Lorentz-scalar, that is, additive to the particle masses. This is why, we accept the last term in (10) to be the mass term, that is, (11) is the imaginary-part mass generated by the interference term of the Cornell potential (2).

Thus, we have the particle real-part (constituent) mass, , and the imaginary-part mass, , originating from the potential. As in case of the eigenmomenta, we introduce the complex mass:

which we use to reconstruct the complex masses of resonances, that is, their centered masses and widths.

4. The Complex Masses of Resonances

Resonances are complex-mass values and can be described by complex numbers. These numbers are important even if one wants to find real solutions of a problem. Using complex numbers, we are getting more than what we insert. Separate consideration of two components of the Cornell potential (2) results in the complex momentum (7) and mass (12) of bound particles. These complex expressions can be combined in one: the mass squared of two free relativistic particles with the complex momenta (7) and masses (12) as ; that is,

The complex-mass expression (13), in which the real and imaginary parts are exact, is also exact. This formula is an Anzatz like the Cornell potential (2). It can be written in another form as

where

The real part (15) exactly coincides with the universal mass formula obtained independently by another method with the use of the two-point Padé approximant [24, 25, 27, 28]. The simple mass formula (15) describes equally well the mass spectra of all and mesons ranging from the (, , ) states up to the heaviest known systems [27, 28]. It has been used to calculate the glueball masses as well [24, 25].

The universal mass formula (15) is very transparent physically, as well as the potential (2) (OGE + linear). It allows us to get an analytic expression for Regge trajectories in the whole region. Transforming (15) we obtain the cubic equation for the angular momentum [24]:

Solution obtained from (17) results in the “saturating” Regge trajectories including the Pomeron [24, 25] and appears to be successful in many applications.

The “saturating” Regge trajectories (see Figure 1) [27, 28] were applied with success to the photoproduction of vector mesons that provide an excellent simultaneous description of the high and low behavior of the cross-sections [33, 34], given an appropriate choice of the relevant coupling constants (JML-model) [35, 36]. It was shown that the hard-scattering mechanism is incorporated in an effective way by using the “saturated” Regge trajectories that are independent of at large momentum transfers [27, 28].

Figure 1: The Chew-Frautschi plots of the leading , , and complex Regge trajectories calculated from (17) [24, 25]. Parameters are found from the fit of combined HERA and data (triangles) [29–31] and glueball candidate (cross) [32]. Solid curve is the trajectory; dotted and dashed lines show the leading and saturating Regge trajectories.

Saturating trajectories have a close phenomenological connection to the quark-antiquark interaction which governs the mesonic structure [25, 27, 28]. They provide an effective way to implement gluon exchange between the quarks forming the exchanged meson [33, 37] and lead to the asymptotic quark counting rules [38] that, model independently, determine the energy behavior of the cross section at large . This approach was successfully adopted to explain the large momentum transfer hadron-hadron interactions, as well as several photon-induced reactions [39]. The pion saturating trajectory ( when ) is in a form that reproduces the reaction around [39].

In general (mathematically), the -matrix is a meromorphic function of complex variable , where the complex -plane is replaced by the two-sheet Riemann surface, , made up of two sheets and , each cut along the positive real axis, , and with placed in front of [5, 13]. The square root of the complex expression (14) gives

where

where and .

The expressions (19) and (20) define the resonance positions which are symmetrically located in the Riemann -surface. Resonance masses arise in complex-conjugate pairs. Poles in the left half-plane correspond to either bound or antibound states [2, 40]. If is a pole in the fourth quadrant of the surface , then is also a pole, but in the third quadrant (antiparticle) [2].

As mentioned above, the masses of the states develop imaginary masses from loop corrections. In this case, the probability density comes from the particle's propagator, with the complex mass, . This formula is related to the particle's decay rate by the optical theorem [3, 40]. That means the total widths are

where the imaginary mass is given by (20). The centered mass, , and the total width, , are process-independent parameters of the resonance.

There is the lack of a precise definition of what is meant by mass and width of resonance. Comprehensive definitions of these resonance's parameters require further investigations. An alternate definition of the resonance's width can be obtained from (16). According to definition (21), the width is given by the negative imaginary-part mass of the resonance's complex mass, .

In the pole approach, the pole position in the complex -plane is given by (1). Comparing (14) with (1), we obtain the equality

The question arises: how should we treat (22)? One choice is

An alternate definition of the resonance's width can be obtained as follows. According to (21), the width is given by the imaginary-part mass of the resonance's complex mass, . The mass and width in (22) are in factorized form (system mass) × (width); that is, . Therefore, omitting the real-mass factor (), one can consider the following definition:

This value depends on the imaginary-part mass of complex mass (12) and quantum numbers. It has the maximum value given by the imaginary mass, ,

As an example, consider the -family resonances of the leading Regge trajectory, . Calculations are given in Table 1, where masses and widths are in MeV. Note that expressions (20) and (23) give very close results, for . Centered masses were calculated with the use of (19).

Table 1: Masses and total widths of the -family resonances.

Meson

(21)
(24)

776
775
149
150
75

1318
1323
107
108
93

1689
1689
161
170
188

1996
1985
255
194
249

—
2234
—
202
294

—
2462
—
205
328

Parameters values in these calculations are found from the best fit to the available data [1]: , GeV², and MeV. The widths (20) and (24) are calculated with the use of (20) and (24). The maximum imaginary mass is MeV. The spin-dependent corrections to the potential (2) have been considered in [28]. Location of the -family resonances in the complex -surface is shown in Figure 2.

Figure 2: The complex-conjugate poles in the Riemann -surface for the family resonances. Stars show location of the complex-mass resonances on the leading Regge trajectory; the imaginary-part component is the resonance half-width. Crosses show calculation results with the use of (20) and (24).

5. Conclusion

A thorough understanding of the physics summarized by the PDG is related to the concept of a resonance. Resonances in hadron physics are never observed directly but through their decay channels that causes problems with their definition. In addition to known definitions of what is meant by mass and width of resonance, we have considered the alternate definition (24).

In contrast to the usual analysis dealing with the scattering problem, we have modeled mesonic resonances to be the quasibound eigenstates of two interacting quarks. We have chosen the Cornell potential to describe the interquark interaction. Using the complex-mass scheme, we have analyzed the exact eigenvalues for the coulombic and linear terms of the potential, separately, and obtained the complex-mass formula for meson resonances. This approach has allowed us to simultaneously describe in the unified way the masses and total widths of resonances in a good agreement with data.

There exists a widely spread belief that the width of heavy hadronic resonances linearly depends on their mass. We have found this inconsistent with an existence of linear Regge trajectories. In our analysis the total width of resonances is finite or decreasing; according to (25), it may have the maximum value, which equals to imaginary-part mass, . The complex-mass expression (13) may have relation to some nonhermitian Hamiltonian.

Note a feature of the -family resonance data. There is a dip (MeV) for the resonance (see Figure 2). This dip is described in our scheme and may have the following explanation. It is known that complex poles of the -matrix always arise in conjugate pairs [40, 41]. Poles in the lower half-plane are complex-conjugated with zeros in upper half-plane. Zeros above the real axis correspond to bound states, and resonances are located in the lower half-plane.

According to (18), the imaginary-part mass can be positive and negative, poles can be located above and below the real -axis in the complex -surface; that is, in bound state or scattering region. This choice depends on boundary conditions. Analysis and numerical calculations using (18)–(24) show that the imaginary-part mass of the resonance is positive; that is, the -state is located in the first quadrant of the complex -surface (, ). All other states of the -family Regge trajectory are embedded in the fourth quadrant (, ). What does it mean?

Explanation to this phenomenon may be as follows. If the -state is located in the first quadrant, it is a quasibound state, which has the zero quantum numbers, that is, (ground state). If we define resonance to be an excited state with nonzero quantum numbers, either or , then -meson is not a resonance but the quasibound state with parallel spins of quarks. This is only difference between and -meson, that is, quark spin orientation. That means that in this sense all -states are mesons, that is, unstable bound states with positive imaginary-part mass.

Resonances represent a very economical way in theoretical description of hadronic reactions at high energies. Such a task is very important nowadays since a great significance of the width of heavy resonances. Our analysis may be important for further development of the string model of hadrons and for improvement of such transport codes as the hadron string dynamics by including the finite width of heavy resonances. The universal complex-mass formula (13) gives a possibility to calculate and predict masses and width for hadrons/resonances listed in the Particle Data Group tables. This is important for experiments at LHC and others.

Acknowledgments

The author thanks Yu. A. Kurochkin for support and constant interest in this work. This work was done in the framework of investigations for the experiment ATLAS (LHC), code 02-0-1081-2009/2013, “physical explorations at LHC” (JINR-ATLAS).

References

J. Beringer, “Review of particle physics,” Physical Review D, vol. 86, no. 1, Article ID 010001, 2012.

N. Moiseyev, “Quantum theory of resonances: calculating energies, widths and cross-sections by complex scaling,” Physics Reports, vol. 302, no. 5-6, pp. 211–293, 1998.

J. R. Taylor, The Quantum Theory of Nonrelativistic Collisions, Dover Publications, 2006.

A. J. F. Siegert, “On the derivation of the dispersion formula for nuclear reactions,” Physical Review, vol. 56, no. 8, p. 750, 1939.

P. D. Hislop and C. Villegas-Blas, “Semiclassical Szego limit of resonance clusters for the hydrogen atom Stark Hamiltonian,” http://arxiv.org/abs/1104.4466

R. Brummelhuis and A. Uribe, “A semi-classical trace formula for Schrödinger operators,” Communications in Mathematical Physics, vol. 136, no. 3, pp. 567–584, 1991.

G. Rupp, S. Coito, and E. Beveren, “Meson spectroscopy: too much excitement and too few excitations,” Acta Physica Polonica B, vol. 5, supplement 3, pp. 1007–1014, 2012.

W. Heitler and N. Hu, “Interpretation of complex roots of the s-matrix,” Nature, vol. 159, no. 4049, pp. 776–777, 1947.

T. Bauer, D. Djukanovic, J. Gegelia, S. Scherer, and L. Tiator, “Complex-mass scheme and resonances in EFT,” in Proceedings of the 8th International Workshop on the Physics of Excited Nucleons (NSTAR '11), vol. 1432 of AIP Conference Proceedings, Reading, Mass, USA, 2011.

A. Denner, S. Dittmaier, M. Roth, and D. Wackeroth, “Predictions for all processes e⁺e⁻ → fermions + γ,” Nuclear Physics B, vol. 560, no. 1–3, pp. 33–65, 1999.

A. Bernicha, G. L. Castro, and J. Pestieau, “Pion-proton scattering and isospin breaking in the Δ⁰ - Δ⁺⁺ system,” Nuclear Physics A, vol. 597, no. 4, pp. 623–635, 1996.

D. Morgan and M. R. Pennington, “Velocity dependence of absolute cross sections for charge capture by Ar⁷⁺ from ground-state and excited-state sodium,” Physical Review A, vol. 48, no. 2, pp. 1185–1188, 1993.

J. Nieves and E. R. Arriola, “Meson resonances at large NC: complex poles vs Breit-Wigner masses,” Physics Letters B, vol. 679, no. 5, pp. 449–453, 2009.

J. D. Bjorken and E. A. Paschos, “Inelastic electron-proton and γ-proton scattering and the structure of the nucleon,” Physical Review, vol. 185, no. 5, pp. 1975–1982, 1969.

G. S. Bali, “QCD forces and heavy quark bound states,” Physics Reports, vol. 343, no. 1-2, pp. 1–136, 2001.

N. Brambilla, A. Pineda, J. Soto, and A. Vairo, “Effective-field theories for heavy quarkonium,” Reviews of Modern Physics, vol. 77, no. 4, pp. 1423–1496, 2005.

H. M. E. Eichten, S. Godfrey, and J. L. Rosner, “Quarkonia and their transitions,” Reviews of Modern Physics, vol. 80, no. 3, pp. 1161–1193, 2008.

J. Sucher, “Confinement in relativistic potential models,” Physical Review D, vol. 51, no. 10, pp. 5965–5966, 1995.

M. N. Sergeenko, “Semiclassical wave equation and exactness of the WKB method,” Physical Review A, vol. 53, no. 6, pp. 3798–3804, 1996.

P. D. B. Collins, An Introduction to Regge Theory & High Energy Physics, Cambridge University Press, Cambridge, Mass, USA, 1977.

C. M. Bender, D. C. Brody, and D. W. Hook, “Quantum effects in classical systems having complex energy,” Journal of Physics A, vol. 41, no. 35, Article ID 352003, 15 pages, 2008.

C. D. P. Dorey and R. Tateo, “The ODE/IM correspondence,” Journal of Physics A, vol. 40, no. 32, p. R205, 2007.

M. N. Sergeenko, “Relativistic semiclassical wave equation and its solution,” Modern Physics Letters A, vol. 12, no. 37, pp. 2859–2871, 1997.

M. N. Sergeenko, “Glueball masses and Regge trajectories for the QCD-inspired potential,” The European Physical Journal C, vol. 72, article 2128, 2012.

M. N. Sergeenko, “Glueballs and the pomeron,” Europhysics Letters, vol. 89, no. 1, p. 11001, 2010.

M. N. Sergeenko, “Classical solution of the wave equation,” International Journal of Modern Physics A, vol. 18, no. 17, p. 1, 2003.

M. N. Sergeenko, “Some properties of Regge trajectories of heavy quarkonia,” Physics of Atomic Nuclei, vol. 56, pp. 365–371, 1993.

M. N. Sergeenko, “An interpolating mass formula and Regge trajectories for light and heavy quarkonia,” Zeitschrift für Physik C, vol. 64, no. 2, pp. 315–322, 1994.

D. Aston, M. Atkinson, R. Bailey, et al., “Photoproduction of ⁰ and ω on hydrogen at photon energies of 20 to 70 GeV,” Nuclear Physics B, vol. 209, no. 1, pp. 56–76, 1982.

J. Breitweg, M. Derrick, D. Krakauer, et al., “Measurement of the diffractive structure function F₂D(4) at HERA,” The European Physical Journal C, vol. 1, no. 1-2, pp. 81–96, 1998.

T. Ahmed, V. Andreev, B. Andrieu, et al., “Total photoproduction cross section measurement at HERA energies,” Physics Letters B, vol. 299, no. 3-4, pp. 374–384, 1993.

A. Kirk, “New effects observed in central production by experiment WA102 at the CERN Omega Spectrometer,” Nuclear Physics A, vol. 663, pp. 608–612, 2000.

J. M. Laget, “Space-time structure of hard-scattering processes,” Physical Review D, vol. 70, no. 5, Article ID 054023, 2004.

M. Battaglieri, M. Brunoldi, R. de Vita, et al., “Photoproduction of the omega meson on the proton at large momentum transfer,” Physical Review Letters, vol. 90, no. 2, Article ID 022002, 2003.

M. Battaglieri, E. Anciant, M. Anghinolfi, et al., “Photoproduction of the ⁰ meson on the proton at large momentum transfer,” Physical Review Letters, vol. 87, no. 17, Article ID 172002, 5 pages, 2001.

L. Morand, D. Doré, M. Garçon, et al., “Deeply virtual and exclusive electroproduction of ω-mesons,” The European Physical Journal A, vol. 24, no. 3, pp. 445–458, 2005.

F. Cano and J. M. Laget, “Compton scattering, vector-meson photoproduction, and the partonic structure of the nucleon,” Physical Review D, vol. 65, no. 7, Article ID 074022, 4 pages, 2002.

S. J. Brodsky and G. R. Farrar, “Scaling laws at large transverse momentum,” Physical Review Letters, vol. 31, no. 18, pp. 1153–1156, 1973.

M. Guidal, J.-M. Laget, and M. Vanderhaeghen, “Pion and kaon photoproduction at high energies: forward and intermediate angles,” Nuclear Physics A, vol. 627, no. 4, pp. 645–678, 1997.

L. D. Landau and E. M. Lifshitz, Quantum Mechanics, Pergamon, 1965.

N. Fernández-García and O. Rosas-Ortiz, “Gamow-Siegert functions and Darboux-deformed short range potentials,” Annals of Physics, vol. 323, no. 6, pp. 1397–1414, 2008.