title: Trigonometric Rosen–Morse potential
text: # Trigonometric Rosen–Morse potential
  
  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
  
      $V_{tRM}^{(a,b,\lambda )}(\chi )={\frac {a(a-\lambda )}{\sin ^{2}\lambda \chi }}-2b\cot \lambda \chi ,\quad \chi ={\frac {r}{R}},\quad \lambda \chi \in \left[0,\pi \right].$         (1)
  
  where $r$ is a relative distance, $\lambda$ is an angle rescaling parameter, and $R$ is so far a matching length parameter. Another parametrization of same potential is
  
      $V_{tRM}^{(a,b,\lambda )}(\chi )={\frac {a(a-\lambda )}{\cos ^{2}\lambda \chi }}-2b\tan \lambda \chi ,\quad \chi ={\frac {r}{R}},\quad \lambda \chi \in \left[-{\frac {\pi }{2}},+{\frac {\pi }{2}}\right],$         (2)
  
  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,
  
      $V(x/d)=B\operatorname {tanh} (x/d)-C\operatorname {sech} ^{2}(x/d),$         (3)
  
  a parallelism that explains the potential's name. The most prominent application concerns the $V_{tRM}^{(\ell +1,b,1)}(\chi )$ parametrization, with $\ell$ non-negative integer, and is due to Schrödinger who intended to formulate the hydrogen atom problem on Albert Einstein's closed universe, $R^{1}\otimes S^{3}$, the direct product of a time line with a three-dimensional closed space of positive constant curvature, the hypersphere $S^{3}$, and introduced it on this geometry in his celebrated equation as the counterpart to the Coulomb potential, a mathematical problem briefly highlighted below.
  
  ## The $V_{tRM}^{(\ell +1,0,1)}(\chi )$ case: Four-dimensional rigid rotator in inertial quantum motion on the three dimensional hypersphere $S^{3}$
  
  The hypersphere is a surface in a four-dimensional Euclidean space, $E_{4}$, and is defined as,
  
      $x_{1}^{2}+x_{2}^{2}+x_{3}^{2}+x_{4}^{2}=R^{2},$         (4)
  
  where $x_{1}$, $x_{2}$, $x_{3}$, and $x_{4}$ are the Cartesian coordinates of a vector in $E_{4}$, and $R$ is termed to as hyper-radius. Correspondingly, Laplace operator in $E_{4}$ is given by,
  
      $\Delta _{4}(x_{1},x_{2},x_{3},x_{4})={\frac {\partial ^{2}}{\partial x_{4}^{2}}}+{\frac {\partial ^{2}}{\partial x_{1}^{2}}}+{\frac {\partial ^{2}}{\partial x_{2}^{2}}}+{\frac {\partial ^{2}}{\partial x_{3}^{2}}}.$         (5)
  
  In now switching to polar coordinates,
  
      ${\begin{aligned}x_{4}=R\cos \chi ,\quad x_{3}=R\sin \chi \cos \theta ,\quad \chi ,\theta &\in \left[0,\pi \right],\\x_{2}=R\sin \chi \sin \theta \sin \varphi ,\quad x_{1}=R\sin \chi \sin \theta \cos \varphi ,\quad \varphi &\in \left[0,2\pi \right],\end{aligned}}$         (6)
  
  one finds the Laplace operator expressed as
  
      $\Delta _{4}(x_{1},x_{2},x_{3},x_{4})\longrightarrow \Delta _{4}(R,\chi ,\theta ,\varphi )={\frac {1}{R^{3}}}{\frac {\partial }{\partial R}}R^{3}{\frac {\partial }{\partial R}}-{\frac {1}{R^{2}}}{\mathcal {K}}^{2}(\chi ,\theta ,\varphi ),$         (7)
      ${\mathcal {K}}^{2}(\chi ,\theta ,\varphi )=-{\frac {1}{\sin ^{2}\chi }}{\frac {\partial }{\partial \chi }}\sin ^{2}\chi {\frac {\partial }{\partial \chi }}+{\frac {{\mathbf {L} }^{2}(\theta ,\varphi )}{\sin ^{2}\chi }}.$         (8)
  
  Here, ${\mathcal {K}}^{2}(\chi ,\theta ,\varphi )$ stands for the squared angular momentum operator in four dimensions, while ${\mathbf {L} }^{2}(\theta ,\varphi )$ is the standard three-dimensional squared angular momentum operator. Considering now the hyper-spherical radius $R$ as a constant, one encounters the Laplace-Beltrami operator on $S^{3}$ as
  
      $\Delta _{S^{3}}(R,\chi ,\theta ,\varphi )=-{\frac {1}{R^{2}}}{\mathcal {K}}^{2}(\chi ,\theta ,\varphi ),\quad R={\mbox{const}}.$         (9)
  
  With that the free wave equation on $S^{3}$ takes the form
  
      ${\begin{aligned}-{\frac {\hbar ^{2}c^{2}}{R^{2}}}\Delta _{S^{3}}(R,\chi ,\theta ,\varphi )Y_{K\ell m}(\chi ,\theta ,\varphi )&={\frac {\hbar ^{2}c^{2}}{R^{2}}}{\mathcal {K}}(\chi ,\theta ,\varphi )Y_{K\ell m}(\chi ,\theta ,\varphi )\\&={\frac {\hbar ^{2}c^{2}}{R^{2}}}K(K+2)Y_{K\ell m}(\chi ,\theta ,\varphi ),\end{aligned}}$         (10)
      $K=0,1,2,...,\quad \ell =0,1,2,...K,\quad |m|=0,1,2,...,\ell .$         (11)
  
  The solutions, $Y_{K\ell m}(\chi ,\theta ,\varphi )$, to this equation are the so-called four-dimensional hyper-spherical harmonics defined as
  
      $Y_{K\ell m}(\chi ,\theta ,\varphi )=\sin ^{\ell }\chi {\mathcal {G}}_{K-\ell }^{\ell +1}(\cos \chi )Y_{\ell }^{m}(\theta ,\varphi ),$         (12)
  
  where ${\mathcal {G}}_{n}^{\ell +1}(\cos \chi )$ are the Gegenbauer polynomials. Changing in (10) variables as
  
      $Y_{K\ell m}(\chi ,\theta ,\varphi )={\frac {\psi _{K\ell }(\chi )}{\sin \chi }}Y_{\ell }(\theta ,\varphi ),$         (13)
  
  one observes that the $\psi _{K\ell }(\chi )$ function satisfies the one-dimensional Schrödinger equation with the $\csc ^{2}\chi$ potential according to
  
      $-{\frac {\hbar ^{2}c^{2}}{R^{2}}}\left[{\frac {{\mathrm {d} }^{2}}{{\mathrm {d} }\chi ^{2}}}+{\frac {\ell (\ell +1)}{\sin ^{2}\chi }}\right]\psi _{K\ell }(\chi )={\frac {\hbar ^{2}c^{2}}{R^{2}}}(K+1)^{2}\psi _{K\ell }(\chi ).$         (14)
  
  The one-dimensional potential in the latter equation, in coinciding with the Rosen–Morse potential in (1) for $a=\ell +1$ and $b=0$, clearly reveals that for integer $a$ values, the first term of this potential takes its origin from the centrifugal barrier on $S^{3}$. Stated differently, the equation (10), and its version (14) describe inertial (free) quantum motion of a rigid rotator in the four-dimensional Euclidean space, $E_{4}$, such as the H Atom, the positronium, etc. whose "ends" trace the large "circles" (i.e. $S^{2}$ spheres) on $S^{3}$.
  
  Now the question arises whether the second term in (1) could also be related in some way to the $S^{3}$ geometry.
  
  ## The $V_{tRM}^{(\ell +1,b,1)}(\chi )$ case: Electric charge confinement on $S^{3}$ and a dipole potential shaped after $\alpha Z\cot \chi$
  
  To the amount the cotangent function solves the Laplace–Beltrami equation on $S^{3}$,
  
      $\Delta _{S^{3}}(R,\chi ,\theta ,\varphi )\cot \chi =0,$         (15)
  
  it represents a fundamental solution on $S^{3}$, 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 $E_{3}$ Laplacian. Due to this analogy, the cotangent function is frequently termed to as "curved Coulomb" potential. Such an interpretation ascribes the cotangent potential to a single charge source, and here lies a severe problem. Namely, while open spaces, as is $E_{3}$, support single charges, in closed spaces single charge can not be defined in a consistent way. 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
  
      ${\begin{aligned}-{\frac {\hbar ^{2}c^{2}}{R^{2}}}\left[\Delta _{S^{3}}(R,\chi ,\theta ,\varphi )-2b\cot \chi \right]\Psi _{K\ell m}(\chi ,\theta ,\varphi )&={\frac {\hbar ^{2}c^{2}}{R^{2}}}\left[{\mathcal {K}}^{2}(\chi ,\theta ,\varphi )-2b\cot \chi \right]\Psi _{K\ell m}(\chi ,\theta ,\varphi )\\&={\frac {\hbar ^{2}c^{2}}{R^{2}}}\left[K(K+2)-{\frac {b^{2}}{(K+1)^{2}}}\right]\Psi _{K\ell m}(\chi ,\theta ,\varphi ),\end{aligned}}$         (16)
  
  which transforms upon the variable change, $\Psi _{K\ell m}(\chi ,\theta ,\varphi )={\frac {U_{K\ell }^{(b)}(\chi )}{\sin \chi }}Y_{\ell }^{m}(\theta ,\varphi )$, into the familiar one-dimensional Schrödinger equation with the $V_{tRM}^{(\ell +1,b,1)}(\chi )$ trigonometric Rosen–Morse potential,
  
      $-{\frac {\hbar ^{2}c^{2}}{R^{2}}}\left[{\frac {{\mathrm {d} }^{2}}{{\mathrm {d} }\chi ^{2}}}+{\frac {\ell (\ell +1)}{\sin ^{2}\chi }}\right]U_{K\ell }^{(b)}(\chi )-2b\cot \chi U_{K\ell }^{(b)}(\chi )={\frac {\hbar ^{2}c^{2}}{R^{2}}}\left[(K+1)^{2}-{\frac {b^{2}}{(K+1)^{2}}}\right]U_{K\ell }^{(b)}(\chi ),$         (17)
  
  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 (16) and (17) do not describe strictly speaking a Hydrogen Atom on $S_{3}$, but rather quantum motion on $S^{3}$ of a light $(e^{+},e^{-})$ dipole perturbed by the dipole potential of another very heavy dipole, like the H Atom, so that the reduced mass, $\mu$, 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 $S^{3}$ of both the Gauss law and the superposition principle for the sake of being capable to formulate electrostatic there. With the cotangent function in (15) as a single-source potential, such can not be achieved. Rather, it is necessary to prove that the cotangent function represents a dipole potential. Such a proof has been delivered in. To understand the line of arguing of it is necessary to go back to the expression for the Laplace operator in (5) and before considering the hyper-radius as a constant, factorize this space into a time line and $S^{3}$. For this purpose, a "time" variable is introduced via the logarithm of the $S^{3}$ radius. Introducing this variable change in (7) amounts to the following Laplacian,
  
      $-\left[{\frac {\partial ^{2}}{\partial x_{4}^{2}}}+{\frac {\partial ^{2}}{\partial x_{1}^{2}}}+{\frac {\partial ^{2}}{\partial x_{2}^{2}}}+{\frac {\partial ^{2}}{\partial x_{3}^{2}}}\right]\longrightarrow -\left[{\frac {1}{R^{3}}}{\frac {\partial }{\partial R}}R^{3}{\frac {\partial }{\partial R}}-{\frac {1}{R^{2}}}{\mathcal {K}}^{2}(\chi ,\theta ,\varphi )\right]$         (18)
      $\qquad \qquad \qquad \qquad \qquad \qquad \quad \;\;\;{\stackrel {R=e^{\tau }}{\longrightarrow }}-e^{2\tau }\left[{\frac {\partial ^{2}}{\partial \tau ^{2}}}-\left({\mathcal {K}}^{2}(\chi ,\theta ,\varphi )+1\right)\right].$         (19)
  
  The $\tau$ parameter is known as "conformal time", and the whole procedure is referred to as "radial quantization". Charge-static is now built up in setting $\tau$=const in (19) and calculating the harmonic function to the remaining piece, the so-called conformal Laplacian, $\Delta _{S^{3}}^{1}(\chi ,\theta ,\varphi )$, on $S^{3}$, which is read off from (19) as
  
      $\Delta _{S^{3}}^{1}(\tau ,\chi ,\theta ,\varphi )|_{\tau =0}=\Delta _{S^{3}}(R,\chi ,\theta ,\varphi )|_{R=1}+1={\mathcal {K}}^{2}(\chi ,\theta ,\varphi )+1,$         (20)
  
  where we have chosen $\tau =0$, equivalently, $R=1$.
  
  Then the correct equation to be employed in the calculation of the fundamental solution is
  
  $\left(\Delta _{S^{3}}^{1}(\tau ,\chi ,\theta ,\varphi )-1\right)|_{\tau =0}{\mathcal {G}}=\Delta _{S^{3}}(R,\chi ,\theta ,\varphi )|_{R=1}{\mathcal {G}}=\delta -{\frac {1}{2\pi ^{2}}}$.
  
  This Green function to $\Delta _{S^{3}}(\tau ,\chi ,\theta ,\varphi )$ has been calculated for example in.
  
  Its values at the respective South and North poles, in turn denoted by ${\mathcal {G}}_{\pi }(\chi )$, and ${\mathcal {G}}_{0}(\chi )$, are reported as
  
      ${\mathcal {G}}_{\pi }(\chi )=-{\frac {1}{4\pi ^{2}}}\chi \cot \chi ,$         (21)
  
  and
  
      ${\mathcal {G}}_{0}(\chi )={\frac {1}{4\pi ^{2}}}\left(\pi -\chi \right)\cot \chi ,$         (22)
  
  From them one can now construct the dipole potential for a fundamental charge ${\mathcal {Q}}$ placed, say, on the North pole, and a fundamental charge of opposite sign, $-{\mathcal {Q}}$, placed on the antipodal South pole of $S^{3}$. The associated potentials, $V_{\pi }(\chi )$ and $V_{0}(\chi )$, are then constructed through multiplication of the respective Green function values by the relevant charges as
  
      $V_{\pi }(\chi )=-{\mathcal {Q}}{\mathcal {G}}_{\pi }(\chi )=-{\frac {\left(-{\mathcal {Q}}\right)}{4\pi ^{2}}}\chi \cot \chi ,$         (23)
      $V_{0}(\chi )={\mathcal {Q}}{\mathcal {G}}_{0}(\chi )=-{\frac {\mathcal {Q}}{4\pi ^{2}}}\left(\chi -\pi \right)\cot \chi .$         (24)
  
  In now assuming validity of the superposition principle, one encounters a Charge Dipole (CD) potential to emerge at a point $\chi$ on $S^{3}$ according to
  
      $V_{CD}(\chi )=V_{\pi }(\chi )+V_{0}(\chi )=-{\mathcal {Q}}{\mathcal {G}}_{\pi }(\chi )+{\mathcal {Q}}{\mathcal {G}}_{0}(\chi )=-{\frac {\mathcal {Q}}{4\pi }}\cot \chi .$         (25)
  
  The electric field to this dipole is obtained in the standard way through differentiation as
  
      $E(\chi )=-{\frac {\partial }{\partial \chi }}V_{CD}(\chi )={\frac {\mathcal {Q}}{4\pi }}{\frac {1}{\sin ^{2}\chi }},$         (26)
  
  and coincides with the precise expression prescribed by the Gauss theorem on $S^{3}$, as explained in. Notice that ${\mathcal {Q}}$ stands for dimension-less charges. In terms of dimensional charges, $q$, related to ${\mathcal {Q}}$ via
  
      ${\mathcal {Q}}={\frac {q}{\sqrt {\hbar c}}},$         (27)
  
  the potential perceived by another charge $(Zq)/{\sqrt {\hbar c}}$, is
  
      $V_{CD}(\chi )=-{\frac {q^{2}Z}{4\pi \hbar c}}\cot \chi .$         (28)
  
  For example, in the case of electrostatic, the fundamental charge $q$ is taken the electron charge, $e$, in which case the special notation of
  
      $\alpha ={\frac {e^{2}}{4\pi \hbar c}}$         (29)
  
  is introduced for the so-called fundamental coupling constant of electrodynamics. In effect, one finds
  
      $V_{CD}(\chi )=-\alpha Z\cot \chi .$         (30)
  
  In Fig. 2 we display the dipole potential $V_{CD}(\chi )$ in (30).
  
  With that, the one-dimensional Schrödinger equation that describes on $S^{3}$ 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
  
      $\left(-{\frac {\hbar ^{2}c^{2}}{R^{2}}}{\frac {{\mathrm {d} }^{2}}{{\mathrm {d} }\chi ^{2}}}+V_{\mbox{tRM}}^{(\ell +1,\alpha Z/2,1)}(\chi )\right)U_{K\ell }^{(\alpha Z)}(\chi )={\frac {\hbar ^{2}c^{2}}{R^{2}}}\left({\epsilon }_{\ell n}^{(\alpha Z)}\right)^{2}U_{K\ell }^{(\alpha Z)}(\chi ),$         (31)
      $\qquad \qquad \qquad \qquad \qquad \quad \;\;\,V_{\mbox{tRM}}^{(\ell +1,\alpha Z/2,1)}(\chi )={\frac {\hbar ^{2}c^{2}}{R^{2}}}{\frac {\ell (\ell +1)}{\sin ^{2}\chi }}-2{\frac {\hbar ^{2}c^{2}}{R^{2}}}\alpha Z\cot \chi ,$         (32)
      $\qquad \qquad \qquad \qquad \qquad \quad \;\;{\frac {\hbar ^{2}c^{2}}{R^{2}}}\left(\epsilon _{\ell n}^{(\alpha Z)}\right)^{2}={\frac {\hbar ^{2}c^{2}}{R^{2}}}(K+1)^{2}-{\frac {\hbar ^{2}c^{2}\alpha ^{2}Z^{2}}{4R^{2}(K+1)^{2}}}.$         (33)
  
  Because of the relationship, $K-\ell =n$, with $n$ being the node number of the wave function, one could change labeling of the wave functions, $U_{K\ell }^{(b)}(\chi )$, to the more familiar in the literature, $U_{\ell n}^{(b)}(\chi )$.
  
  In eqs. (31)-(32) one recognizes the one-dimensional wave equation with the trigonometric Rosen–Morse potential in (1) for $a=\ell +1$ and $2b=\alpha Z$.
  
     In this way, the cotangent term of the trigonometric Rosen–Morse potential could be derived from the Gauss law on $S^{3}$ 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 $\csc ^{2}\chi$ term of this potential has been generated by the kinetic energy operator on $S^{3}$. In this manner, the complete trigonometric Rosen–Morse potential could be derived from first principles.
  
  Back to Schrödinger's work, the $S^{3}$ hyper-radius for the H Atom has turned out to be very big indeed, and of the order of $10^{-3}cm$. 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} 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 (31)-(32), given by
  
      $U_{K\ell }^{(\alpha Z)}(\chi )=N_{K\ell }\sin ^{K+1}\chi e^{-{\frac {\alpha Z}{2(K+1)}}}R_{K-\ell }^{{\frac {\alpha Z}{K+1)}},-(K+1)}(\cot \chi ),\quad K-\ell =n,$         (34)
  
  where $R_{n}^{\alpha \beta }(x)$ stand for the Romanovski polynomials.
  
  ## 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, one employs as a simulation cell the hyper spherical surface $S^{3}$ in (4). As already mentioned above, the basic object on $S^{3}$ 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, $+q$ and $-q$. The potential of a bi-charge is calculated by solving on $S_{3}$ the Poisson equation,
  
      $\Delta _{S^{3}}(\chi )V_{\mbox{b-ch}}(\chi )=-{\frac {4\pi q}{R^{3}}}\left[\delta _{S^{3}}(\chi ,\chi _{0})-\delta _{S_{3}}(\chi ,{\bar {\chi }}_{0})\right].$         (35)
  
  Here, $\chi _{0}$ is the angular coordinate of a charge $q$ placed at angular position $\chi _{0}$, read off from the North pole, while ${\bar {\chi }}_{0}$ stands for the anti-podal to $\chi _{0}$ angular coordinate of the position, at which the charge of opposite signs is placed in the Southern hemisphere. The solution found,
  
      $V_{\mbox{b-ch}}(\chi )={\frac {q}{R}}\cot \chi ,$         (36)
  
  equals the potential in (30), modulo conventions regarding the charge signs and units. It provides an alternative proof to that delivered by the equations (19)-(30) of the fact that the cotangent function on $S^{3}$ has to be associated with the potential generated by a charge dipole. In contrast, the potentials in the above equations (23), and (24), have been interpreted in as due to so called single "pseudo-charge" sources, where a "pseudo-charge" is understood as the association of a point charge $q$ with a uniform neutralizing background of a total charge, $-q$. The pseudo-charge potential, $V_{\mbox{psd-ch}}$, solves $\Delta _{S^{3}}V_{\mbox{psd-ch}}=\delta -{\frac {1}{2\pi ^{2}}}$. Therefore, the bi-charge potential is the difference between the potentials of two antipodal pseudo-charges of opposite signs.
  
  ## Application to color confinement and the physics of quarks
  
  The confining nature of the cotangent potential in (28) 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 $(r)$, blue $(b)$, and green $(g)$, while anti-quarks carry the corresponding anti-colors, anti-red $({\bar {r}})$, anti-blue $({\bar {b}})$, or anti-green $({\bar {g}})$, 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 $g_{s}(Q^{2})$, are not constants, but depend on the values, $Q^{2}$, of the transferred momentum, giving rise to the so-called, running of the strong coupling constant, $\alpha _{s}(Q^{2})$, in which case the Gauss law becomes more involved. However, at low momentum transfer, near the so-called infrared regime, the momentum dependence of the color charge significantly weakens, and in starting approaching a constant value,
  
      $g_{s}(Q^{2}){\stackrel {Q^{2}\to 0}{\longrightarrow }}{\mbox{const}},$         (37)
  
  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 $S^{3}$, it has been shown in, that a potential, there denoted by $V_{CCD}(\chi )$, and obtained from the one in (28) through the replacement,
  
      $\alpha Z\longrightarrow \alpha _{s}N_{c},\quad \alpha _{s}={\frac {g_{s}^{2}}{4\pi }},$         (38)
  
  i.e. the potential
  
      $V_{CCD}(\chi )=-\alpha _{s}N_{c}\cot \chi ,$         (39)
  
  where $N_{c}$ is the number of colors, is the adequate one for the description of the spectra of the light mesons with masses up to $\sim 2500MeV$. Especially, the hydrogen like degeneracies have been well captured. This because the potential, in being a harmonic function to the Laplacian on $S^{3}$, has same symmetry as the Laplacian by itself, a symmetry that is defined by the isometry group of $S^{3}$, i.e. by $SO(4)$, the maximal compact group of the conformal group $SO(2,4)$. For this reason, the potential in (39), as part of $V_{tRM}^{(\ell +1,\alpha _{s}N_{c}/2,1)}(\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 $(q)$-anti-quark $({\bar {q}})$ color dipole in quantum motion on an $S^{3}$ geometry, and gets perturbed by the dipole potential in (39), generated by and other color dipole, such as a gluon $(g)$-anti-gluon $({\bar {g}})$, as visualized in Fig. 3.
  
  The $S^{3}$ geometry could be viewed as the unique closed space-like geodesic of a four-dimensional hyperboloid of one sheet, ${\mathbf {H} }_{1}^{4}$, 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, $dS_{4}$. 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 ${\mathbf {H} }_{1}^{4}$ geodesics can give rise to barriers describing resonances transmitted through them. An illustrative example for the application of the color confining dipole potential in (39) to meson spectroscopy is given in Fig. 4. It should be pointed out that the potentials in the above equations (23) and (24) have been alternatively derived in, from Wilson loops with cusps, predicting their magnitude as $\alpha _{s}N_{c}/(4\pi ^{2})$, and in accord with (38).
  
  The potential in (39) has furthermore been used in in the Dirac equation on $S^{3}$, 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.
  
  ## Applicability of $V_{tRM}^{\left(\ell +1,b,1\right)}$ to phase transitions
  
  The property of the trigonometric Rosen-Morse potential, be it in the parametrization with $b=\alpha Z/2$ in eq. (32) which is of interest to electrodynamics, or in the $b=\alpha _{s}N_{c}/2$ parametrization of interest to QCD from the previous section, qualifies it to studies of phase transitions in systems with electromagnetic or strong interactions on hyperspherical "boxes" of finite volumes . The virtue of such studies lies in the possibility to express the temperature, $T$, as the inverse, $T=1/R$, to the radius $R$ of the hypersphere. For this purpose, knowledge on the partition function (statistical mechanics), here denoted by ${\mathcal {Z}}(R,b)$, of the potential under consideration is needed. In the following we evaluate ${\mathcal {Z}}(R,b)$ for the case of the Schrödinger equation on $S^{3}$ with linear energy (here in units of MeV),
  
      ${\begin{aligned}-{\Big{\frac {\hbar ^{2}c^{2}}{2\mu c^{2}R^{2}}}\Delta _{S^{3}}(R,\chi ,\theta ,\varphi )|_{\theta =\pi /2,\varphi =0}&-2b{\frac {\hbar c}{R}}\cot \chi {\Big )}\psi (\chi )=E_{K}\psi (\chi ),\\E_{K}={\frac {\hbar ^{2}c^{2}}{2\mu c^{2}R^{2}}}(K+1)^{2}&-{\frac {\hbar c}{R}}{\frac {b^{2}}{(K+1)^{2}}},\end{aligned}}$         (40)
  
  where $\mu c^{2}$ 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,
  
      ${\begin{aligned}{\mathcal {Z}}(R,b)&=\sum _{K=0}^{K=\infty }(K+1)^{2}e^{-\beta \left(E_{K}-E_{0}\right)}\\&=1+e^{\beta E_{0}}\sum _{K=1}^{K=\infty }(K+1)^{2}e^{-\beta E_{K}}.\end{aligned}}$         (41)
  
  Here, the thermodynamic beta is defined as $\beta =(k_{B}T)^{-1}$ with $k_{B}$ standing for the Boltzmann constant. In evaluating ${\mathcal {Z}}(R,b)$ it is useful to recall that with the increase of $K$ the second term on the right hand side in (40) becomes negligible compared to the term proportional $(K+1)^{2}$, a behavior which becomes even more pronounced for the choices, $2b=\alpha Z$, and $2b=\alpha _{s}N_{c}$. In both cases $b$ is much smaller compared to the corresponding dimensionless factor, $(\hbar c)/(2\mu c^{2}R)$, multiplying $(\hbar c/R)(K+1)^{2}$. For this reason the partition function under investigation might be well approximated by,
  
      ${\begin{aligned}{\mathcal {Z}}(R,b)&=\sum _{K=0}^{K=\infty }(K+1)^{2}e^{-\beta \left(E_{K}-E_{0}\right)}\\&=1+e^{-\beta \left(E_{1}-E_{0}\right)}+e^{\beta E_{0}}\sum _{K=2}^{K=\infty }(K+1)^{2}e^{-\beta E_{K}}\\&\approx 1+e^{-\beta \left(E_{1}-E_{0}\right)}+e^{\beta E_{0}}\int _{0}^{\infty }x^{2}e^{-ax^{2}}{\mathrm {d} }x\\&=1+e^{-\beta \left(E_{1}-E_{0}\right)}+e^{\beta E_{0}}{\frac {\sqrt {\pi }}{4(a)^{\frac {3}{2}}}},\quad a=\beta {\frac {\hbar ^{2}c^{2}}{2\mu c^{2}R^{2}}}.\end{aligned}}$         (42)
  
  Along same lines, the partition function for the $b=\alpha Z/2$ parametrization corresponding to the Hydrogen atom on $S^{3}$ has been calculated in, where a more sophisticated approximation has been employed. When transcribed to the current notations and units, the partition function in presents itself as,
  
      ${\begin{aligned}{\mathcal {Z}}\left(R,{\frac {\alpha Z}{2}}\right)\approx 1+e^{\beta E_{0}}\int _{0}^{\infty }x^{2}e^{-ax^{2}+{\frac {p}{x^{2}}}},\quad p=\beta {\frac {\hbar c}{R}}\alpha ^{2}Z^{2}.\end{aligned}}$         (43)
  
  The infinite integral has first been treated by means of partial integration giving,
  
      ${\begin{aligned}\int _{0}^{\infty }x^{2}e^{-ax^{2}+{\frac {p}{x^{2}}}}&=\int _{0}^{\infty }x^{2}e^{-ax^{2}}e^{\frac {p}{x^{2}}}{\mathrm {d} }x\\&=-{\frac {1}{2a}}\int _{0}^{\infty }xe^{\frac {p}{x^{2}}}{\mathrm {d} }e^{-ax^{2}}\\&=-{\frac {1}{2a}}xe^{\frac {p}{x^{2}}}e^{-ax^{2}}|_{0}^{\infty }+{\frac {1}{2a}}\int _{0}^{\infty }e^{-ax^{2}}{\mathrm {d} }\,\left(xe^{\frac {p}{x^{2}}}\right)\\&={\frac {1}{2a}}\int _{0}^{\infty }e^{-ax^{2}+{\frac {p}{x^{2}}}}{\mathrm {d} }x+{\frac {2p}{2a}}\int _{0}^{\infty }e^{-ax^{2}+{\frac {p}{x^{2}}}}{\mathrm {d} }\left({\frac {1}{x}}\right).\end{aligned}}$         (44)
  
  Then the argument of the exponential under the sign of the integral has been cast as,
  
      ${\begin{aligned}-ax^{2}+{\frac {p}{x^{2}}}=-z^{2}+2i{\sqrt {ap}},\quad z={\sqrt {a}}x+i{\frac {\sqrt {p}}{x}},\end{aligned}}$         (45)
  
  thus reaching the following intermediate result,
  
      ${\begin{aligned}{\mathcal {Z}}\left(R,{\frac {\alpha Z}{2}}\right)\approx 1+{\frac {1}{2a}}e^{\beta E_{0}}e^{2i{\sqrt {ap}}}\int _{0}^{\infty }e^{-z^{2}}{\mathrm {d} }\left(x+{\frac {2p}{x}}\right).\end{aligned}}$         (46)
  
  As a next step the differential has been represented as
  
      ${\begin{aligned}x+{\frac {2p}{x}}={\frac {1}{2}}\left({\frac {1}{\sqrt {a}}}-2i{\sqrt {p}}\right)z+{\frac {1}{2}}\left({\frac {1}{\sqrt {a}}}+2i{\sqrt {p}}\right)z^{\ast },\end{aligned}}$         (47)
  
  an algebraic manipulation which allows to express the partition function in (46) in terms of the ${\mbox{erf}}(u)$ function of complex argument according to,
  
      ${\begin{aligned}{\mathcal {Z}}\left(R,{\frac {\alpha Z}{2}}\right)\approx 1+{\frac {1}{4a}}e^{\beta E_{0}}\times \left[\left({\frac {1}{\sqrt {a}}}-2i{\sqrt {p}}\right)e^{2i{\sqrt {ap}}}\int _{\Gamma }e^{-z^{2}}{\mathrm {d} }z+\left({\frac {1}{\sqrt {a}}}+2i{\sqrt {p}}\right)e^{-2i{\sqrt {ap}}}\int _{\Gamma }e^{-\left(z^{\ast }\right)^{2}}{\mathrm {d} }z^{\ast }\right],\\\end{aligned}}$         (48)
  
  where $\Gamma$ is an arbitrary path on the complex plane starting in zero and ending in $u\to \infty$. For more details and physical interpretations, see.
  
  