Trigonometric Rosen–Morse potential explained
The trigonometric Rosen–Morse potential, named after the physicists Nathan Rosen and Philip M. Morse, is among the exactly solvable quantum mechanical potentials.
Definition
In dimensionless units and modulo additive constants, it is defined as [1]
where
is a relative distance,
is an angle rescaling parameter, and
is so far a matching length parameter. Another parametrization of same potential is
which is the trigonometric version of a one-dimensional hyperbolic potential introduced in molecular physics by Nathan Rosen and Philip M. Morse and given by,[2]
a parallelism that explains the potential's name. The most prominent application concerns the
parametrization, with
non-negative integer, and is due to Schrödinger
[3] who intended to formulate the
hydrogen atom problem on
Albert Einstein's closed universe,
, the
direct product of a time line with a three-dimensional closed space of positive constant curvature, the
hypersphere
, and introduced it on this geometry in his celebrated equation as the counterpart to the Coulomb potential, a mathematical problem briefly highlighted below.
The hypersphere is a surface in a four-dimensional Euclidean space,
, and is defined as,
where
,
,
, and
are the
Cartesian coordinates of a vector in
, and
is termed to as hyper-radius. Correspondingly,
Laplace operator in
is given by,
In now switching to polar coordinates,
one finds the Laplace operator expressed as
Here,
{lK}2(\chi,\theta,\varphi)
stands for the squared angular momentum operator in four dimensions, while
is the standard three-dimensional
squared angular momentum operator. Considering now the hyper-spherical radius
as a constant, one encounters the
Laplace-Beltrami operator on
as
With that the free wave equation on
takes the form
The solutions,
YK\ell(\chi,\theta,\varphi)
, to this equation are the so-called four-dimensional hyper-spherical harmonics defined as
where
are the
Gegenbauer polynomials. Changing in variables as
one observes that the
function satisfies the one-dimensional
Schrödinger equation with the
potential according to
The one-dimensional potential in the latter equation, in coinciding with the Rosen–Morse potential in for
and
, clearly reveals that for integer
values, the first term of this potential takes its origin from the centrifugal barrier on
. Stated differently, the equation, and its version describe inertial (free) quantum motion of a rigid rotator in the four-dimensional
Euclidean space,
, such as the H Atom, the
positronium, etc. whose "ends" trace the large "circles" (i.e.
spheres) on
.
Now the question arises whether the second term in could also be related in some way to the
geometry.
To the amount the cotangent function solves the Laplace–Beltrami equation on
,
it represents a fundamental solution on
, a reason for which Schrödinger considered it as the counterpart to the Coulomb potential in flat space, by itself a fundamental solution to the
Laplacian. Due to this analogy, the
cotangent function is frequently termed to as "curved Coulomb" potential.
[4] Such an interpretation ascribes the cotangent potential to a single charge source, and here lies a severe problem. Namely, while open spaces, as is
, support single charges, in closed spaces single charge can not be defined in a consistent way.
[5] Closed spaces are necessarily and inevitably charge neutral meaning that the minimal fundamental
degrees of freedom allowed on them are charge dipoles (see Fig. 1).
For this reason, the wave equation
which transforms upon the variable change,
\PsiK\ell(\chi,\theta,\varphi)=
, into the familiar one-dimensional
Schrödinger equation with the
trigonometric Rosen–Morse potential,
in reality describes quantum motion of a charge dipole perturbed by the field due to another charge dipole, and not the motion of a single charge within the field produced by another charge. Stated differently, the two equations and do not describe strictly speaking a Hydrogen Atom on
, but rather quantum motion on
of a light
dipole perturbed by the dipole potential of another very heavy dipole, like the H Atom, so that the reduced mass,
, would be of the order of the electron mass and could be neglected in comparison with the energy.
In order to understand this decisive issue, one needs to focus attention to the necessity of ensuring validity on
of both the Gauss law and the
superposition principle for the sake of being capable to formulate electrostatic there. With the cotangent function in as a single-source potential, such can not be achieved.
[6] Rather, it is necessary to prove that the cotangent function represents a dipole potential. Such a proof has been delivered in.
[7] To understand the line of arguing of
[7] it is necessary to go back to the expression for the Laplace operator in and before considering the hyper-radius as a constant, factorize this space into a time line and
. For this purpose, a "time" variable is introduced via the logarithm of the
radius.
[8] Introducing this variable change in amounts to the following Laplacian,
The
parameter is known as "conformal time", and the whole procedure is referred to as "radial
quantization".
[8] Charge-static is now built up in setting
=const in and calculating the harmonic function to the remaining piece, the so-called conformal Laplacian,
, on
, which is read off from as
where we have chosen
, equivalently,
.
Then the correct equation to be employed in the calculation of the fundamental solution is
\left(
(\tau,\chi,\theta,\varphi)-1\right)|\tau=0
(R,\chi,\theta,\varphi)|R=1{lG}=\delta-
.
This Green function to
(\tau,\chi,\theta,\varphi)
has been calculated for example in.
[9] Its values at the respective South and North poles, in turn denoted by
, and
, are reported as
and
From them one can now construct the dipole potential for a fundamental charge
placed, say, on the North pole, and a fundamental charge of opposite sign,
, placed on the antipodal South pole of
. The associated potentials,
and
, are then constructed through multiplication of the respective Green function values by the relevant charges
[10] as
In now assuming validity of the superposition principle, one encounters a Charge Dipole (CD) potential to emerge at a point
on
according to
The electric field to this dipole is obtained in the standard way through differentiation as
and coincides with the precise expression prescribed by the Gauss theorem on
, as explained in.
[6] Notice that
stands for dimension-less charges. In terms of dimensional charges,
, related to
via
the potential perceived by another charge
, is
For example, in the case of electrostatic, the fundamental charge
is taken the electron charge,
, in which case the special notation of
is introduced for the so-called fundamental coupling constant of electrodynamics. In effect, one finds
In Fig. 2 we display the dipole potential
in .
With that, the one-dimensional Schrödinger equation that describes on
the quantum motion of an
electric charge dipole perturbed by the trigonometric Rosen–Morse potential, produced by another electric charge dipole, takes the form of
Because of the relationship,
, with
being the node number of the wave function, one could change labeling of the
wave functions,
, to the more familiar in the literature,
.
In eqs. - one recognizes the one-dimensional wave equation with the trigonometric Rosen–Morse potential in for
and
.
In this way, the cotangent term of the trigonometric Rosen–Morse potential could be derived from the Gauss law on
in combination with the superposition principle, and could be interpreted as a dipole potential generated by a system consisting of two opposite fundamental charges. The centrifugal
term of this potential has been generated by the kinetic energy operator on
. In this manner, the complete trigonometric Rosen–Morse potential could be derived from first principles.
Back to Schrödinger's work,[3] the
hyper-radius for the H Atom has turned out to be very big indeed, and of the order of
. This is by eight orders of magnitudes larger than the H Atom size. The result has been concluded from fitting magnetic dipole elements to hydrogen hyper-fine structure effects (see
[11] } and reference therein). The aforementioned radius is sufficiently large to allow approximating the hyper-sphere locally by plane space in which case the existence of single charge still could be justified. In cases in which the hyper spherical radius becomes comparable to the size of the system, the charge neutrality takes over. Such an example will be presented in section 6 below.
Before closing this section, it is in order to bring the exact solutions to the equations -, given by
where
stand for the
Romanovski polynomials.
[12] [13] [14] Application to Coulomb fluids
Coulomb fluids consist of dipolar particles and are modelled by means of direct numerical simulations. It is commonly used to choose cubic cells with periodic boundary conditions in conjunction with Ewald summation techniques. In a more efficient alternative method pursued by,[15] [16] one employs as a simulation cell the hyper spherical surface
in . As already mentioned above, the basic object on
is the electric charge dipole, termed to as "bi-charge" in
fluid dynamics, which can be visualized classically as a rigid "dumbbell" (rigid rotator) of two antipodal charges of opposite signs,
and
. The potential of a bi-charge is calculated by solving on
the
Poisson equation,
Here,
is the angular coordinate of a charge
placed at angular position
, read off from the North pole, while
stands for the anti-podal to
angular coordinate of the position, at which the charge of opposite signs is placed in the Southern hemisphere. The solution found,
equals the potential in, modulo conventions regarding the charge signs and units. It provides an alternative proof to that delivered by the equations - of the fact that the cotangent function on
has to be associated with the potential generated by a charge dipole. In contrast, the potentials in the above equations, and, have been interpreted in
[15] as due to so called single "pseudo-charge" sources, where a "pseudo-charge" is understood as the association of a point charge
with a uniform neutralizing background of a total charge,
.The pseudo-charge potential,
, solves
.Therefore, the bi-charge potential is the difference between the potentials of two antipodal pseudo-charges of opposite signs.
[17] Application to color confinement and the physics of quarks
The confining nature of the cotangent potential in finds an application in a phenomenon known from the physics of strong interaction which refers to the non-observability of free quarks, the constituents of the hadrons. Quarks are considered to possess three fundamental internal degree of freedom, conditionally termed to as "colors", red
, blue
, and green
, while anti-quarks carry the corresponding anti-colors, anti-red
, anti-blue
, or anti-green
, meaning that the non-observability of free quarks is equivalent to the non-observability of free color-charges, and thereby to the "color neutrality" of the
hadrons. Quark "colors" are the fundamental degrees of freedom of the
Quantum Chromodynamics (QCD), the
gauge theory of strong interaction. In contrast to the
Quantum Electrodynamics, the
gauge theory of the electromagnetic interactions, QCD is a
non-Abelian theory which roughly means that the "color" charges, denoted by
, are not constants, but depend on the values,
, of the transferred momentum, giving rise to the so-called, running of the strong coupling constant,
, in which case the
Gauss law becomes more involved.
[18] However, at low momentum transfer, near the so-called infrared regime, the momentum dependence of the color charge significantly weakens,
[19] and in starting approaching a constant value,
drives the Gauss law back to the standard form known from Abelian theories. For this reason, under the condition of color charge constancy, one can attempt to model the color neutrality of hadrons in parallel to the neutrality of Coulomb fluids, namely, by considering quantum color motions on closed surfaces. In particular for the case of the hyper-sphere
, it has been shown in,
[20] that a potential, there denoted by
, and obtained from the one in through the replacement,
i.e. the potential
where
is the number of colors, is the adequate one for the description of the spectra of the light mesons with masses up to
. Especially, the hydrogen like degeneracies have been well captured. This because the potential, in being a
harmonic function to the
Laplacian on
, has same symmetry as the Laplacian by itself, a symmetry that is defined by the isometry group of
, i.e. by
, the maximal compact group of the conformal group
. For this reason, the potential in, as part of
| (\ell+1,\alphasNc/2,1) |
V | |
| tRM |
(\chi)
, accounts not only for
color confinement, but also for
conformal symmetry in the infrared regime of QCD. Within such a picture, a
meson is constituted by a quark
-anti-quark
color dipole in quantum motion on an
geometry, and gets perturbed by the dipole potential in, generated by and other color dipole, such as a gluon
-anti-gluon
, as visualized in Fig. 3.
The
geometry could be viewed as the unique closed space-like geodesic of a four-dimensional
hyperboloid of one sheet,
, foliating outside of the causal Minkowski light-cone the space-like region, assumed to have one more spatial dimension, this in accord with the so-called
de Sitter Special Relativity,
.
[21] Indeed, potentials, in being instantaneous and not allowing for time orderings, represent virtual, i.e. acausal processes and as such can be generated in one-dimensional
wave equations upon proper transformations of virtual quantum motions on surfaces located outside the causal region marked by the
Light Cone. Such surfaces can be viewed as
geodesics of the surfaces foliating the space like region. Quantum motions on open
geodesics can give rise to barriers describing resonances transmitted through them.
[7] An illustrative example for the application of the color confining dipole potential in to
meson spectroscopy is given in Fig. 4. It should be pointed out that the potentials in the above equations and have been alternatively derived in,
[22] [23] from Wilson loops with cusps, predicting their magnitude as
, and in accord with .
The potential in has furthermore been used in [24] in the Dirac equation on
, and has been shown to predict realistic electromagnetic nucleon form-factors and related constants such as mean square electric-charge and magnetic-dipole radii, proton and nucleon magnetic dipole moments and their ratio, etc.
The property of the trigonometric Rosen-Morse potential, be it in the parametrization with
in eq. (32)which is of interest to electrodynamics, or in the
parametrization of interest to QCD from theprevious section, qualifies it to studies of phase transitions in systems with electromagnetic or strong interactions onhyperspherical "boxes" of finite volumes
[25] .
[26] The virtue of such studies lies in thepossibility to express the temperature,
, as the inverse,
, to the radius
of the hypersphere. For this purpose, knowledge on the
partition function (statistical mechanics), here denoted by
, of the potential under consideration is needed. In the following we evaluate
for the case of the Schrödinger equation on
with linear energy (here in units of MeV),
where
is the reduced mass of the two-body system under consideration. The
partition function (statistical mechanics) for this energy spectrum is defined in the standard way as,
Here, the thermodynamic beta is defined as
with
standing for the
Boltzmann constant. In evaluating
it is useful to recall that with the increase of
the second term on the right hand side in becomes negligible compared to the term proportional
, abehavior which becomes even more pronounced for the choices,
, and
. In both cases
ismuch smaller compared to the corresponding dimensionless factor,
, multiplying
. For this reason the partition function under investigation might be well approximated by,
Along same lines, the partition function for the
parametrization corresponding to the
Hydrogen atom on
has been calculated in,
[27] where a more sophisticated approximation has been employed. When transcribed tothe current notations and units, the partition function in
[27] presents itself as,
The infinite integral has first been treated bymeans of partial integration giving,
Then the argument of the exponential under the sign of the integral has been cast as,
thus reaching the following intermediate result,
As a next step the differential has been represented as
an algebraic manipulation which allows to express the partition function in in terms of the
function of complex argument according to,
where
is an arbitrary path on the complex plane starting in zero and ending in
. For more details and physical interpretations, see.
[27] [28] See also
Notes and References
- Book: F. . Cooper . A. . Khare . U. P. . Sukhatme . 2001 . Supersymmetry in Quantum Mechanics . World Scientific . Singapore . 978-9-81-024612-9.
- N. . Rosen . P. M. . Morse . On the Vibrations of Polyatomic Molecules . Phys. Rev. . 42 . 2 . 210 . 1932 . 10.1103/PhysRev.42.210. 1932PhRv...42..210R .
- E. . Schrödinger . The Factorization of the Hypergeometric Equation . Proc. Roy. Irish Acad. A . 47 . 53–54 . 1941 . 20488434 .
- A. O. . Barut . R. . Wilson . On the dynamical group of the Kepler problem in a curved space of constant curvature . Phys. Lett. A . 110 . 7–8 . 1985 . 351 . 10.1016/0375-9601(85)90052-0. 1985PhLA..110..351B .
- Book: L. D. . Landau . E. M. . Lifschitz . The Classical Theory of Fields . Vol. 2 of Course of Theoretical Physics . 3rd. . Pergamon Press . 1971 . 335 . 978-0-08-016019-1.
- P. . Pouria . Modification of Coulomb's law in closed spaces . Am. J. Phys. . 78 . 4 . 2010 . 403 . 10.1119/1.3272020 . 0912.0225 . 2010AmJPh..78..403P . 115166698 .
- M. . Kirchbach . C. B. . Compean . Modelling duality between bound and resonant meson spectra by means of free quantum motions on the de Sitter space-time
. Eur. Phys. J. A . 52 . 7 . 2016 . 210 . 10.1140/epja/i2016-16210-3 . 1608.05041 . 2016EPJA...52..210K . 118370681 .
- S. . Fubini . A. J. . Hanson . R. . Jackiw . New Approach to Field Theory . Phys. Rev. D . 7 . 6 . 1973 . 1732 . 10.1103/PhysRevD.7.1732 . 1973PhRvD...7.1732F .
- B. . Alertz . Electrodynamics in Robertson-Walker spacetimes . Annales de l'Institut Henri Poincaré A . 53 . 319 . 3 . 1990.
- Book: Kellogg, O. D. . Foundations of Potential Theory . Dover . New York . 1953 . 978-0-48-660144-1.
- N. . Bessis . G. . Bessis . R. . Shamseddine . Atomic fine structure in a space of constant curvature . J. Phys. A: Math. Gen. . 15 . 10 . 3131 . 1982 . 10.1088/0305-4470/15/10/017. 1982JPhA...15.3131B .
- V. . Romanovski . Sur quelques classes nouvelles de polynomes orthogonaux . C. R. Acad. Sci. Paris . 188 . 1929 . 1023 . French.
- E. J. . Routh . On some properties of certain solutions of a differential equation of second order . Proc. London Math. Soc. . 16 . 1884 . 245 . 10.1112/plms/s1-16.1.245.
- A. P. . Raposo . H. J. . Weber . D. E. . Álvarez Castillo . M. . Kirchbach . Romanovski polynomials in selected physics problems . Cent. Eur. J. Phys. . 5 . 3 . 253–284 . 2007 . 10.2478/s11534-007-0018-5. 0706.3897 . 2007CEJPh...5..253R . 119120266 .
- J. M. . Caillol . A new potential for the numerical simulations of electrolyte solutions on a hypersphere . J. Chem. Phys. . 99 . 11 . 8953 . 1993 . 10.1063/1.465565. 1993JChPh..99.8953C .
- J. M. . Caillol . M. . Trulsson . A new dipolar potential for numerical simulations of polar fluids on the 4D hypersphere . J. Chem. Phys. . 141 . 12 . 124111 . 2014 . 10.1063/1.4896181. 25273416 . 1407.7739 . 2014JChPh.141l4111C . 5818844 .
- J. M. . Caillol. D. . Levesque . Numerical simulations of homogeneous and inhomogeneous ionic systems: An efficient alternative to the Ewald method . J. Chem. Phys. . 94 . 1 . 8953 . 1991 . 10.1063/1.460326. 1991JChPh..94..597C .
- M. . Serna . K. . Cahill . Riemannian gauge theory and charge quantization . Journal of High Energy Physics . 2003 . 10 . 054 . 2003 . 10.1088/1126-6708/2003/10/054. hep-th/0205250 . 2003JHEP...10..054S . 6573051 .
- A. . Deur . V. . Burkert . J. P. . Chen . W. . Korsch . Determination of the effective strong coupling constant
from CLAS spin structure function data . Phys. Lett. B . 665 . 5 . 349–351 . 2008 . 10.1016/j.physletb.2008.06.049. 2008PhLB..665..349D . 0803.4119 . 118442505 .
- M. . Kirchbach . C. B. . Compean . Addendum to: Modelling duality between bound and resonant meson spectra by means of free quantum motions on the de Sitter space-time
. Eur. Phys. J. A . 53 . 4 . 65 . 2017 . 10.1140/epja/i2017-12269-6. 2017EPJA...53...65K . 195307166 .
- R. . Aldrovandi . J. P. . Beltrán Almeida . J. G. . Pereira . de Sitter special relativity . Class. Quantum Grav. . 24 . 1385–1404 . 6 . 2007 . 10.1088/0264-9381/24/6/002. gr-qc/0606122 . 2007CQGra..24.1385A . 11703342 .
- A. V. . Belitsky . A. S. . Gorsky . G. P. . Korchemsky . Gauge/string duality for QCD conformal operators . Nucl. Phys. B . 667 . 1–2 . 3–54 . 2003 . 10.1016/S0550-3213(03)00542-X. hep-th/0304028 . 2003NuPhB.667....3B . 15273221 .
- A. S. . Gorsky . Spin chains and gauge-string duality . Theor. Math. Phys. . 142 . 2 . 153 . 2005 . 10.1007/s11232-005-0042-9. hep-th/0308182 . 446917 .
- M. . Kirchbach . C. B. . Compean . Proton's electromagnetic form factors from a non-power confinement potential . Nucl. Phys. A . 980 . 32 . 2018 . 10.1016/j.nuclphysa.2018.09.083. 2018NuPhA.980...32K . 1810.03665 . 119494972 .
- O. . Aharony . J. . Marsano . S. . Minwalla . K. . Papadodimas . M. . Van Raamsdonk . The Hagedorn deconfinement phase transition in weakly coupled large N gauge theories. Adv. Theor. Math. Phys.. 8 . 4 . 603–696 . 2004. 10.4310/ATMP.2004.v8.n4.a1. hep-th/0310285.
- S. . Hands . T. J. . Hollowood . J. C. . Myers. QCD with Chemical Potential in a Small Hyperspherical Box. Journal of High Energy Physics. 2010 . 7 . 86 . 2010. 10.1007/JHEP07(2010)086. 1003.5813. 2010JHEP...07..086H . 119258554 .
- S. M. . Blinder . Canonical partition function for the hydrogen atom in curved space. J. Math. Chem.. 19. 43 . 1996. 10.1007/BF01165129. 2027.42/43064 . 25491310 . free.
- Aram . Bahroz Brzo . David . Alvarez-Castillo. Thermodynamic Properties of a Quantum Gas of Mesons Considered as Color Charge Dipoles in a Three-Dimensional Space of Constant Curvature. Modern Physics Letters A. 36. 15. 2150095-1. 2021. 10.1142/S0217732321500954. 2010.14642. 225094121 .