In this paper, we present a novel theoretical analytical perform further investigation for the exact solvability of non-relativistic quantum spectrum systems for modified inverse-square potential (m.i.s.) potential is discussed by means Boopp’s shift method instead to solving deformed Schrödinger equations with star product, in the framework of both noncommutativity (two -three) dimensional real space and phase (NC: 2D-RSP) and (NC: 3D-RSP). The exact corrections for excited states are found straightforwardly for interactions in one-electron atoms by means of the standard perturbation theory. Furthermore, the obtained corrections of energies are depended on four infinitesimals parameters ( , ) and ( , ), which are induced by position-position and momentum-momentum noncommutativity, (NC: 2D-RSP) and (NC: 3D-RSP), respectively, in addition to the discreet atomic quantum numbers: and (the angular momentum quantum number) and we have also shown that, the usual states in ordinary two and three dimensional spaces are cancelled and has been replaced by new degenerated sub-states in the new quantum symmetries of (NC: 2D-RSP) and (NC: 3D-RSP).

Keywords: The inverse-square potential; Noncommutative space; Phase; Star product; Boopp’s shift method.

MIS: Modified Inverse Square potential; NC: 2D-3D-RSP: Noncommutativity (two-three) Dimensional Real Space Phase; CCRs: Canonical Commutations Relations; NNCCRs: New Noncommutative Canonical Commutations Relations; MSE: Modified Schrödinger Equations.

It is well-known, that, the modern quantum mechanics, satisfied a big successful in the last few years, for describing atoms, nuclei, and molecules and their spectral behaviors based on three fundamental equations: Schrödinger, Klein-Gordon and Dirac. Schrödinger equation rest the first and the latest in terms of interest, it is playing a crucial role in devising well-behaved physical models in different fields of physics and chemists, many potentials are treated within the framework of nonrelativistic quantum mechanics based on this equation in two, three and D generalized spaces [1-32], the quantum structure based to the ordinary canonical commutations relations (CCRs) in both Schrödinger and Heisenberg (the operators are depended on time) pictures (CCRs), respectively, as:
$\left[x_i,p_j\right]=i{\delta }_{ij}\text{ and }\left[x_i,x_j\right] = \left[p_i,p_j\right] = 0$ (1.1)

$\left[x_i\left(t\right),p_j\left(t\right)\right] = i{\delta }_{ij}\text{ and }\left[x_i\left(t\right),x_j\left(t\right)\right] = \left[p_i\left(t\right),p_j\left(t\right)\right] = 0$ (1.2)

Where the two operators $\left(x_i\left(t\right),p_i\left(t\right)\right)$ in Heisenberg picture are related to the corresponding two operators $\left(x_i,p_i\right)$ in Schrödinger picture from the two projections relations:

$x_i\left(t\right)=\exp (iH\left(t-t_0\right))x_i\exp (-iH\left(t-t_0\right))\text{ and }{p}_i\left(t\right)=\exp (iH\left(t-t_0\right))p_i\exp (-iH\left(t-t_0\right))$ (1.3)

Here denote to the ordinary quantum Hamiltonian operator, recently, much considerable effort has been expanded on the solutions of Schrödinger, Dirac and Klein-Gordon equations to noncommutative quantum mechanics, the present paper investigates first the present new quantum structure which based to new noncommutative canonical commutations relations (NNCCRs) in both Schrödinger and Heisenberg pictures, respectively, as follows [33-60]:

$\begin{array}{l}\left[{\widehat{x}}_i\stackrel{\ast }{,}{\widehat{p}}_j\right]=i{\delta }_{ij},\left[{\widehat{x}}_i\stackrel{\ast }{,}{\widehat{x}}_j\right]=i{\theta }_{ij}\text{ and }\left[{\widehat{p}}_i\stackrel{\ast }{,}{\widehat{p}}_j\right]=i{\overline{\theta }}_{ij}\\ \left[{\widehat{x}}_i\left(t\right)\stackrel{\ast }{,}{\widehat{p}}_j\left(t\right)\right]=i{\delta }_{ij},\left[{\widehat{x}}_i\left(t\right)\stackrel{\ast }{,}{\widehat{x}}_j\left(t\right)\right]=i{\theta }_{ij}\text{ and }\left[{\widehat{p}}_i\left(t\right)\stackrel{\ast }{,}{\widehat{p}}_j\left(t\right)\right]=i{\overline{\theta }}_{ij}\end{array}$ (1.4)

Where the two new operators $\left({\widehat{x}}_i\left(t\right),{\widehat{p}}_i\left(t\right)\right)$ in Heisenberg picture are related to the corresponding two new operators $\left({\widehat{x}}_i,{\widehat{p}}_i\right)$ in Schrödinger picture from the two projections relations:

${\widehat{x}}_i\left(t\right)=\exp (i{H}_{nc}\left(t-t_0\right))*{\widehat{x}}_i*\exp (-i{H}_{nc}\left(t-t_0\right))\text{ and }{\widehat{p}}_i\left(t\right)=\exp (i{H}_{nc}\left(t-t_0\right))*{\widehat{p}}_i*\exp (-i{H}_{nc}\left(t-t_0\right))$ (1.5)

Here $H_{nc}$ denote to the new quantum Hamiltonian operator in the symmetries of (NC: 2D-RSP) and (NC: 3D-RSP). The very small two parameters ${\theta }^{\mu \nu }$ and ${\overline{\theta }}^{\mu \nu }$ (compared to the energy) are elements of two ant symmetric real matrixes and (*) denote to the new star product, which is generalized between two arbitrary functions $f\left(x,p\right)$ and $g\left(x,p\right)$ to $\left(f\ast g\right)\left(x,p\right)$ instead of the usual product $\left(fg\right)\left(x,p\right)$ in ordinary (two-three) dimensional spaces [39-63]:

$\left(f\ast g\right)\left(x,p\right) \equiv \exp (\frac{i}{2}{\theta }^{\mu \nu }{\partial }_{\mu }^x{\partial }_{\nu }^x+\frac{i}{2}{\overline{\theta }}^{\mu \nu }{\partial }_{\mu }^p{\partial }_{\nu }^p)\left(fg\right)\left(x,p\right) ={ (fg-\frac{i}{2}{\theta }^{\mu \nu }{\partial }_{\mu }^{x}f{\partial }_{\nu }^{x}g-\frac{i}{2}{\overline{\theta }}^{\mu \nu }{\partial }_{\mu }^{p}f{\partial }_{\nu }^{p}g)\left(x,p\right)|}_{\left(x^{\mu }=x^{\nu },p^{\mu }=p^{\nu }\right)}+O\left({\theta }^2,{\overline{\theta }}^2\right)$ (2)

Where the two covariant derivatives $\left({\partial }_{\mu }^{x}f\left(x,p\right),{\partial }_{\mu }^{p}f\left(x,p\right)\right)$ are denotes to the $\left(\frac{\partial f\left(x,p\right)}{\partial x^{\mu }},\frac{\partial f\left(x,p\right)}{\partial p^{\mu }}\right)$ , respectively, and the two following terms [ $-\frac{i}{2}{\theta }^{\mu \nu }{\partial }_{\mu }^{x}f\left(x,p\right){\partial }_{\nu }^{x}g\left(x,p\right)$ , $-\frac{i}{2}{\overline{\theta }}^{\mu \nu }{\partial }_{\mu }^{p}f\left(x,p\right){\partial }_{\nu }^{p}g\left(x,p\right)$ ] are induced by (space-space) and (phase-phase) noncommutativity properties, respectively, a Boopp's shift method can be used, instead of solving any quantum systems by using directly star product procedure [39-66]:

$\left[{\widehat{x}}_i,{\widehat{x}}_j\right]=i{\theta }_{ij}\begin{array}{cc}\text{and}& \left[{\widehat{p}}_i,{\widehat{p}}_j\right]=i{\overline{\theta }}_{ij}\end{array}$ (3.1)

The, four generalized positions and momentum coordinates in the noncommutative quantum mechanics $\left(\widehat{x},\widehat{y}\right)$ and $\left({\widehat{p}}_x,{\widehat{p}}_y\right)$ are depended with corresponding four usual generalized positions and momentum coordinates in the usual quantum mechanics $\left(x,y\right)$ and $\left(p_x,p_y\right)$ by the following four relations [32-55]:

$\{\begin{array}{l}\widehat{x}=x-\frac{\theta }{2}{p}_y\begin{array}{cc},& \widehat{y}=y+\frac{\theta }{2}{p}_x\end{array}\\ {\widehat{p}}_x=p_x+\frac{\overline{\theta }}{2}y\begin{array}{cc}\text{ and}& {\widehat{p}}_y=p_x-\frac{\overline{\theta }}{2}x\end{array}\end{array}$ (3.2)

$\{\begin{array}{l}\widehat{x}=x-\frac{{\theta }_{12}}{2}{p}_y-\frac{{\theta }_{13}}{2}{p}_z,\widehat{y}=y-\frac{{\theta }_{21}}{2}{p}_x-\frac{{\theta }_{23}}{2}{p}_z\\ \text{and }\widehat{z}=z-\frac{{\theta }_{31}}{2}{p}_x-\frac{{\theta }_{32}}{2}{p}_y\end{array}$ (3.3)

and

$\{\begin{array}{l}{\widehat{p}}_x=p_x-\frac{{\overline{\theta }}_{12}}{2}y-\frac{{\overline{\theta }}_{13}}{2}z,{\widehat{p}}_y=p_y-\frac{{\overline{\theta }}_{{}_{21}}}{2}x-\frac{{\overline{\theta }}_{{}_{23}}}{2}z\\ \text{and }{\widehat{p}}_z=p_z-\frac{{\overline{\theta }}_{{}_{31}}}{2}x-\frac{{\overline{\theta }}_{{}_{32}}}{2}y\end{array}$ (3.4)

The non-vanish 9-commutators in (NC-2D: RSP) and (NC-3D: RSP) can be determined as follows:

$\begin{array}{l}\left[\widehat{x},{\widehat{p}}_x\right]=\left[\widehat{y},{\widehat{p}}_y\right]=i,\\ \left[\widehat{x},\widehat{y}\right]=i{\theta }_{12}\text{ and }\left[{\widehat{p}}_x,{\widehat{p}}_y\right]=i{\overline{\theta }}_{12}\end{array}$ (3.5)

and

$\begin{array}{l}\left[\widehat{x},{\widehat{p}}_x\right]=\left[\widehat{y},{\widehat{p}}_y\right]=\left[\widehat{z},{\widehat{p}}_z\right]=i,\\ \left[\widehat{x},\widehat{y}\right]=i{\theta }_{12},\left[\widehat{x},\widehat{z}\right]=i{\theta }_{13},\left[\widehat{y},\widehat{z}\right]=i{\theta }_{23}\\ \left[{\widehat{p}}_x,{\widehat{p}}_y\right]=i{\overline{\theta }}_{12},\left[{\widehat{p}}_y,{\widehat{p}}_z\right]=i{\overline{\theta }}_{23},\left[{\widehat{p}}_x,{\widehat{p}}_z\right]=i{\overline{\theta }}_{13}\end{array}$ (3.6)

Which allow us to getting the two operators ${\widehat{r}}^2$ and ${\widehat{p}}^2$ on a noncommutative two dimensional space-phase as follows [32-48]:

${\widehat{r}}^2=r^2-\theta L_z\text{ and }{\widehat{p}}^2=p^2+\overline{\theta }{L}_z$ (4.1)

${\widehat{r}}^2=r^2-\overrightarrow{L}\overrightarrow{\Theta }\text{ and }\frac{{\widehat{p}}^2}{\text{2}\mu }=\frac{p^2}{\text{2}\mu }+\frac{\overrightarrow{L}\overrightarrow{\overline{\theta }}}{\text{2}\mu }$ (4.2)

Where the two couplings $L\Theta$ and $\overrightarrow{L}\overrightarrow{\overline{\theta }}$ are given by, respectively:

$L\Theta \equiv L_x{\Theta }_{12}+L_y{\Theta }_{23}+L_z{\Theta }_{13}\text{ and }\overrightarrow{L}\overrightarrow{\overline{\theta }}\equiv L_x{\overline{\theta }}_{12}+L_y{\overline{\theta }}_{23}+L_z{\overline{\theta }}_{13}$ (5.1)

It is-well known, that, in quantum mechanics, the three components ( $L_x$ , $L_y$ and $L_z$ ) are determined, in Cartesian coordinates:

$L_x=y{p}_z-z{p}_y,L_y={\text{zp}}_{\text{x}}{\text{-xp}}_{\text{z}}\text{ and }{L}_z=x{p}_y-y{p}_x$ (5.2)

The study of inverse-square potential has now become a very interest field due to their applications in different fields [1], this work is aimed at obtaining an analytic expression for the eigenenergies of a inverse-square potential in (NC: 2D-RSP) and (NC: 3D-RSP) using the generalization Boopp’s shift method based on mentioned formalisms on above equations to discover the new symmetries and a possibility to obtain another applications to this potential in different fields, it is important to notice that, this potential was studied, in ordinary two dimensional spaces, by authors Shi-Hai Dong and Guo-Hua Sun of the Ref. the Schrödinger equation with a Coulomb plus inverse-square potential in D dimensions [1]. The organization scheme of the study is given as follows: In next section, we briefly review the Schrödinger equation with inverse-square potential on based to Ref. [1]. The Section 3, devoted to studying the (two-three) deformed Schrödinger equation by applying both Boopp's shift method to the inverse-square potential. In the fourth section and by applying standard perturbation theory we find the quantum spectrum of the excited states in (NC-2D: RSP) and (NC-3D: RSP) for spin-orbital interaction. In the next section, we derive the magnetic spectrum for studied potential. In the sixth section, we resume the global spectrum and corresponding noncommutative Hamiltonian for inverse-square potential. Finally, the important results and the conclusions are discussed in last section.

Here we will firstly describe the essential steps, which gives the solutions of time independent Schrödinger equation for a fermionic particle like electron of rest mass and its energy moving in inverse-square potential [1]:

$V\left(r\right)=\frac{A}{r^2}-\frac{B}{r}$ (6)

Where A and B are two positive constant coefficients. The above potential is the sum of Colombian $\left(-\frac{B}{r}\right)$ and inverse-square potential $\left(\frac{A}{r^2}\right)$ , if we insert this potential into the non-relativistic Schrödinger equation; we obtain the following equation, in two and three dimensional spaces, respectively, as follows:

$\left\{-\frac{{{\displaystyle \hslash }}^2}{2m_0}\left[\frac{1}{r}\frac{\partial }{\partial r}\left(r\frac{\partial }{\partial r}\right)+\frac{1}{r^2}\frac{{{\displaystyle \partial }}^2}{\partial {{\displaystyle \varphi }}^2}\right]-\frac{A}{r^2}+\frac{B}{r}\right\}\Psi \left(r,\varphi \right)={{\rm E}}_{2d}\Psi \left(r,\varphi \right)$ (7.1)

$\left\{-\frac{{{\displaystyle \hslash }}^2}{2m_0}\left[\frac{1}{{{\displaystyle r}}^2}\frac{\partial }{\partial r}\left({{\displaystyle r}}^2\frac{\partial }{\partial r}\right)+\frac{1}{{{\displaystyle r}}^2\sin \theta }\frac{\partial }{\partial \theta }\left(\sin \theta \frac{\partial }{\partial \theta }\right)+\frac{1}{{{\displaystyle r}}^2{{\displaystyle \left(\sin \theta \right)}}^2}\frac{{{\displaystyle \partial }}^2}{\partial {{\displaystyle \varphi }}^2}\right]-\frac{A}{r^2}+\frac{B}{r}\right\}\Psi \left(r,\theta ,\varphi \right)={{\rm E}}_{3d}\Psi \left(r,\theta ,\varphi \right)$ (7.2)

Here $\Psi \left(r,\varphi \right)$ and $\Psi \left(r,\theta ,\varphi \right)$ is the solution in the (2-3) dimensional in (polar and spherical) coordinates, the complete wave function ( $\Psi \left(r,\varphi \right)$ and $\Psi \left(r,\theta ,\varphi \right)$ separated as follows:

$\Psi \left(r,\varphi \right)={{\displaystyle R}}_l\left(r\right)e^{\pm i\varphi }$ (8.1)

and

$\Psi \left(x\right)={{\displaystyle R}}_l\left(r\right){{\displaystyle Y}}_l^l\left(\theta ,\varphi \right)$ (8.2)

Substituting eq. (8.1) and (8.2) into eq. (7.1) and (7.2), we obtain the radial function ${{\displaystyle R}}_l\left(r\right)$ satisfied the following equation, in (two-three) dimensional spaces [1]:

$\frac{{{\displaystyle d^{2}R}}_l\left(\rho \right)}{d{\rho }^2}+\frac{2}{\rho }\frac{{{\displaystyle dR}}_l\left(\rho \right)}{d\rho }+\left(-\frac{1}{4}+\frac{\tau }{\rho }-\frac{2A+l^2}{{\rho }^2}\right){{\displaystyle R}}_l\left(r\right) = 0$ (9.1)

$\frac{1}{{{\displaystyle r}}^2}\frac{\partial }{\partial r}\left({{\displaystyle r}}^2\frac{\partial }{\partial r}\right){{\displaystyle R}}_l\left(r\right)+\left[2\left({\rm E}-V\left(r\right)\right)-\frac{l\left(l+1\right)}{{{\displaystyle r}}^2}\right]{{\displaystyle R}}_l\left(r\right) = 0$ (9.2)

Here $\rho =r\sqrt{-8E}$

and $\tau =B\sqrt{-\frac{1}{2E}}$ . The proposed solutions of eqs. (9.1) and (9.2) are determined from the unifed relation:

${{\displaystyle R}}_l\left(\rho \right)={{\displaystyle \rho }}^{\lambda }{{\displaystyle e}}^{-\frac{\rho }{2}}F\left(\rho \right)$ (10)

where $\lambda =\frac{2-D+\sqrt{2A+k^2}}{2}$ and $k=2l+D-2$ . We Companie between eqs. (9.1), (9.2) and (10) to obtains [1]:

$\rho \frac{{{\displaystyle d}}^{2}F\left(\rho \right)}{d{{\displaystyle \rho }}^2}+\left(2\lambda +D-1-\rho \right)+\frac{dF\left(\rho \right)}{d\rho }+\left(\tau -\lambda -\frac{D-1}{2}\right)F\left(\rho \right)=0$ (11)

The confluent hypergeometric functions $\phi \left(\lambda -\tau +1/2,2\lambda +1;\rho \right)$ are present the solutions of eq. (11) [1] :

$R\left(\rho \right)={\rm N}{{\displaystyle \rho }}^{\lambda }{{\displaystyle e}}^{-\frac{\rho }{2}}\phi \left(\lambda -\tau +(D-1)/2,2\lambda +D-1;\rho \right)$ (12)

The constraint conditions on the potential parameters are determined from relations [1]:

$\begin{array}{l}\tau -\lambda -(D-1)/2=n^{\prime }=0,1,2,\mathrm{.......}\\ n=n^{\prime }+\frac{\kappa }{2}-D/2+2=n^{\prime }+l+1\\ \tau ={\rm B}\sqrt{\frac{1}{-2{\rm E}}}=n-l-1+\lambda +(D-1)/2\end{array}$ (13)

The normalized wave functions $\Psi \left(\rho ,\varphi \right)$ expressed in terms of the radial functions and spherical harmonic functions read as [1]:

$\Psi \left(\rho ,\varphi \right)=\left(\frac{4{\rm B}}{2n-2m+2s_2-1}\right){\left(\frac{\left(n-m-1\right)!}{\left(2n-2m+2s_2-1\right)\left(n-m+2s_2-1\right)!}\right)}^{1/2}{{\displaystyle \rho }}^{s_2}{{\displaystyle e}}^{\frac{-\rho }{2}}{{\displaystyle L}}_{n-m-1}^{2s_2}\left(\rho \right)\exp \left(\pm im\varphi \right)$ (14.1)

$\Psi \left(\rho \right)={{\displaystyle \left(\frac{2B}{n}\right)}}^{\frac{3}{2}}{{\displaystyle \left[\frac{\left(n-l-1\right)!}{2n\left(n+l\right)!}\right]}}^{\frac{1}{2}}{{\displaystyle \rho }}^l{{\displaystyle e}}^{-\frac{\rho }{2}}{{\displaystyle L}}_{n-l-1}^{2l+1}\left(\rho \right){{\displaystyle Y}}_l^l\left(\theta ,\varphi \right)$ (14.2)

And the corresponding eigenvalues ${\rm E}\left(n,l,D\right)$ is determined from relation [1]:

${\rm E}\left(n,l,D\right)=-\frac{2{\rm B}}{\left(2n-2l-1+\sqrt{8{\rm A}+{{\displaystyle \kappa }}^2}\right)}=\{\begin{array}{l}-\frac{{{\displaystyle 2{\rm B}}}^2}{{{\displaystyle \left(2n-2m-1+\sqrt{2A+m^2}\right)}}^2}\begin{array}{cc}\text{ for}& D=2\end{array}\\ -2{{\displaystyle B}}^2\left\{{{\displaystyle \left(2n\right)}}^{-2}-\frac{8A}{\kappa }{{\displaystyle \left(2n\right)}}^{-3}+\frac{16{{\displaystyle A}}^2}{{{\displaystyle \kappa }}^3}{{\displaystyle \left(2n\right)}}^{-3}+\frac{48{{\displaystyle A}}^2}{{{\displaystyle \kappa }}^2}{{\displaystyle \left(2n\right)}}^{-4}-\mathrm{...}\right\}\begin{array}{cc}\text{ for}& D=3\end{array}\end{array}$ (15)

The rest of this section is devoted to the reapply of some essential properties of generalized Laguerre polynomials ${{\displaystyle L}}_n^{\left(\beta \right)}\left(\rho \right)$ which are given by:

${{\displaystyle L}}_n^{\left(\beta \right)}\left(\rho \right)\equiv \frac{1}{2\prod i}{\displaystyle \oint \frac{\exp \left(-\frac{\rho t}{1-t}\right)}{{\left(1-t\right)}^{\beta +1}{t}^{n+1}}dt}$ (16)

Where $\beta$ is integer, this can be taking the exciplicitly mathematically forms [1,65,66,67]:

${{\displaystyle L}}_n^{\left(\beta \right)}\left(\rho \right)=\frac{\beta \left(\beta +n+1\right)}{n!\beta \left(\beta +1\right)}{}_1{{\displaystyle F}}_1\left(-n,\beta +1;\rho \right)$ (17)

The Laguerre polynomials may be defined in terms of hypergeometric functions ${}_1{{\displaystyle F}}_1\left(-n,\beta +1;\rho \right)$ , specifically the confluent hyper geometric functions, as:

${}_1{{\displaystyle F}}_1\left(-n,\beta +1;\rho \right)={\displaystyle \sum _{n=0}^{\infty }\frac{a^{\left(n\right)}{\rho }^n}{b^{\left(n\right)}n!}}$ (18.1)

Where $a^{\left(n\right)}$ is the Pochhammer symbol, which can be takes the particulars values $a^{\left(0\right)}=0$ and $a^{\left(n\right)}=a\left(a+1\right)\mathrm{.....}\left(a+n-1\right)$ , it is important to notice that, the hypergeometric functions have another common notation $\Phi \left(a,b,\rho \right)$ which considered as a function of a, $b=0,-1,-2,\mathrm{...}$ , and the variable . The generalized Laguerre polynomial can also be defined by the following equation:

${{\displaystyle L}}_n^{\left(\beta \right)}\left(\rho \right)={\displaystyle \sum _{i=0}^n\frac{{\left(-1\right)}^i\left(\begin{array}{l}n+\beta \\ n-i\end{array}\right)}{\left(\text{i}\right)\text{!}}}{\rho }^i$ (18.2)

This section is devoted to constructing of non relativistic modified Schrödinger equations (m.s.e) in both (NC-2D: RSP) and (NC-3D: RSP) for (m.i.s.) potential; to achieve this subject, we apply the essentials following steps [32-48]:

1. Ordinary two dimensional Hamiltonian operators ( ${\widehat{H}}_{is2}\left(p_i,x_i\right)$ , ${\widehat{H}}_{is3}\left(p_i,x_i\right)$ ) will be replaced by new two dimensional Hamiltonian operators ( ${\widehat{H}}_{nc2-is}\left({\widehat{p}}_i,{\widehat{x}}_i\right)$ , ${\widehat{H}}_{nc3-is}\left({\widehat{p}}_i,{\widehat{x}}_i\right)$ ),
2. Ordinary complex wave function $\Psi \left(\overrightarrow{r}\right)$ will be replacing by new complex wave function $\stackrel{⌢}{\Psi }\left(\overleftrightarrow{\stackrel{⌢}{r}}\right)$ ,
3. Ordinary energies $E\left(n,l,2\right)$ and $E\left(n,l,3\right)$ will be replaced by new values $E_{nc2-is}\left(n,l,2,\mathrm{...}\right)$ and $E_{nc3-is}\left(n,l,3,\mathrm{...}\right)$ , respectively.

And the last step corresponds to replace the ordinary old product by new star product , which allow us to constructing the modified two dimensional Schrödinger equation in both (NC-2D: RSP) and (NC-3D: RSP) as for (m.i.s.) potential:

${\widehat{H}}_{nc2-is}\left({\widehat{p}}_i,{\widehat{x}}_i\right)\ast \stackrel{⌢}{\Psi }\left(\overleftrightarrow{\stackrel{⌢}{r}}\right)=E_{nc2-is}\left(n,l,2,\mathrm{...}\right)\stackrel{⌢}{\Psi }\left(\overleftrightarrow{\stackrel{⌢}{r}}\right)$ (19.1)

and

${\widehat{H}}_{nc3-is}\left({\widehat{p}}_i,{\widehat{x}}_i\right)\ast \stackrel{⌢}{\Psi }\left(\overleftrightarrow{\stackrel{⌢}{r}}\right)=E_{nc3-is}\left(n,l,3,\mathrm{...}\right)\stackrel{⌢}{\Psi }\left(\overleftrightarrow{\stackrel{⌢}{r}}\right)$ (19.2)

In order to use the ordinary product without star product, with new vision, as mentioned before, we apply the Boopp’s shift method on the above eqs. (19.1) and (19.2) to obtain two reduced Schrödinger in both (NC-2D: RSP) and (NC-3D: RSP) for (m.i.s.) potential:

$H_{nc2-is}\left({\widehat{p}}_i,{\widehat{x}}_i\right)\psi \left(\overrightarrow{r}\right)=E_{nc2-is}\left(n,l,2,\mathrm{...}\right)\psi \left(\overrightarrow{r}\right)$ (20.1)

and

$H_{nc3-is}\left({\widehat{p}}_i,{\widehat{x}}_i\right)\psi \left(\overrightarrow{r}\right)=E_{nc3-is}\left(n,l,3,\mathrm{...}\right)\psi \left(\overrightarrow{r}\right)$ (20.2)

Where the new operators of Hamiltonian $H_{nc2-is}\left({\widehat{p}}_i,{\widehat{x}}_i\right)$ and $H_{nc3-is}\left({\widehat{p}}_i,{\widehat{x}}_i\right)$ can be expressed in three general varieties: both noncommutative space and noncommutative phase (NC-2D: RSP, NC-3D: RSP), only noncommutative space (NC-2D: RS, NC-3D: RS) and only noncommutative phase (NC: 2D-RP, NC: 3D-RP) as, respectively:

$H_{nc(2-3)-is}\left({\widehat{p}}_i,{\widehat{x}}_i\right)\equiv H\left(p_x+\frac{\overline{\theta }}{2}y,p_y-\frac{\overline{\theta }}{2}x,x-\frac{\theta }{2}{p}_y,y+\frac{\theta }{2}{p}_x\right)\begin{array}{cc}\text{ for}& \text{NC-2D: RSP and NC-3D: RSP}\end{array}$ (21.1)

$H_{nc(2-3)-is}\left({\widehat{p}}_i,{\widehat{x}}_i\right)\equiv H\left(p_x,p_y,x-\frac{\theta }{2}{p}_y,y+\frac{\theta }{2}{p}_x\right)\begin{array}{cc}\text{ for}& \text{NC-2D: RS}\end{array}\text{ and NC-3D: RS}$ (22.2)

$H_{nc(2-3)-is}\left({\widehat{p}}_i,{\widehat{x}}_i\right)\equiv H\left(p_x+\frac{\overline{\theta }}{2}y,p_x-\frac{\overline{\theta }}{2}x,x,y\right)\begin{array}{cc}\text{ for}& \text{NC-2D: RP}\end{array}\text{ and NC-3D: RP}$ (22.3)

In recently work, we are interest with the first variety (21.1), after straightforward calculations, we can obtain the five important terms, which will be use to determine the (m.i.s.) potential in (NC: 2D- RSP) and (NC: 3D-RSP), respectively, as:

$\frac{A}{{{\displaystyle \widehat{r}}}^2}=\frac{A}{{{\displaystyle r}}^4}+\frac{A\theta L_z}{{{\displaystyle r}}^4},\frac{B}{\widehat{r}}=\frac{B}{r}-\frac{B\theta L_z}{{{\displaystyle 2r}}^3}\text{ and }\frac{{\widehat{p}}^2}{2m_0}=\frac{p^2}{2m_0}+\frac{\overrightarrow{L}\overrightarrow{\overline{\theta }}}{2m_0}$ (23)

and

$\frac{A}{{{\displaystyle \widehat{r}}}^2}=\frac{A}{{{\displaystyle r}}^4}+\frac{A\overrightarrow{L}\overrightarrow{\Theta }}{{{\displaystyle r}}^4},\frac{B}{\widehat{r}}=\frac{B}{r}-\frac{B\overrightarrow{L}\overrightarrow{\Theta }}{{{\displaystyle 2r}}^3}\text{ and }\frac{{\widehat{p}}^2}{2m_0}=\frac{p^2}{2m_0}+\frac{\overline{\theta }{L}_z}{2m_0}$ (24)

Which allow us to obtaining the global potential operator $H_{nc2-is}\left({\widehat{p}}_i,{\widehat{x}}_i\right)$ and $H_{nc3-is}\left({\widehat{p}}_i,{\widehat{x}}_i\right)$ for (m.i.s) potential in both (NC: 2D-RSP) and (NC: 3D-RSP), respectively, as:

$H_{nc2-is}\left({\widehat{p}}_i,{\widehat{x}}_i\right)=\frac{A}{r^2}-\frac{B}{r}+\frac{p^2}{2m_0}+\frac{\overline{\theta }{L}_z}{2m_0}+\left(\frac{A}{r^4}-\frac{B}{2r^3}\right)\theta L_z$ (25.1)

and

$H_{nc3-is}\left({\widehat{p}}_i,{\widehat{x}}_i\right)=\frac{A}{r^2}-\frac{B}{r}+\frac{p^2}{2m_0}+\frac{\overrightarrow{L}\overrightarrow{\overline{\theta }}}{2m_0}+\left(\frac{A}{r^4}-\frac{B}{2r^3}\right)\overrightarrow{L}\overrightarrow{\Theta }$ (25.2)

It’s clearly, that the four first terms are given the ordinary inverse-square potential and kinetic energy in (2D-3D) spaces, while the rest terms are proportional’s with infinitesimals parameters ( , ) and ( , ), thus, we can considered as a perturbations terms, we noted by ${\widehat{H}}_{2-pert}\left(r,A,B,\theta ,\overline{\theta }\right)$ and ${\widehat{H}}_{3-pert}\left(r,A,B,\Theta ,\overline{\theta }\right)$ for (NC: 2D-RSP) and (NC: 3D-RSP) symmetries, respectively, as:

${\widehat{H}}_{2-pert}\left(r,A,B,\theta ,\overline{\theta }\right)=\frac{L_z\overline{\theta }}{2m_0}+\left(\frac{A}{r^4}-\frac{B}{2r^3}\right)\theta L_z$ (26.1)

and

${\widehat{H}}_{3-pert}\left(r,A,B,\Theta ,\overline{\theta }\right)=\frac{\overrightarrow{L}\overrightarrow{\overline{\theta }}}{2m_0}+\left(\frac{A}{r^4}-\frac{B}{2r^3}\right)\overrightarrow{L}\overrightarrow{\Theta }$ (26.2)

The exact spin-orbital hamiltonian for (m.i.s.) potential in both (NC: 2D- RSP) and (NC: 3D- RSP) symmetries for one-electron atoms

 Again, the perturbative two terms ${\widehat{H}}_{2-pert}\left(r,A,B,\theta ,\overline{\theta }\right)$  and ${\widehat{H}}_{3-pert}\left(r,A,B,\Theta ,\overline{\theta }\right)$  can be rewritten to the equivalent physical form for (m.i.p.) potential:
$${\widehat{H}}_{2-pert}\left(r,A,B,\theta ,\overline{\theta }\right)=\left\{\frac{\overline{\theta }}{2m_0}+\theta \left(\frac{A}{r^4}-\frac{B}{2r^3}\right)\right\}\overleftrightarrow{S}\overleftrightarrow{L}$$ (26.3)
${\widehat{H}}_{3-pert}\left(r,A,B,\Theta ,\overline{\theta }\right)=\left\{\frac{\overline{\theta }}{2m_0}+\Theta \left(\frac{A}{r^4}-\frac{B}{2r^3}\right)\right\}\overleftrightarrow{S}\overleftrightarrow{L}$ (26.4)
Furthermore, the above perturbative terms ${\widehat{H}}_{2-pert}\left(r,A,B,\theta ,\overline{\theta }\right)$  and ${\widehat{H}}_{3-pert}\left(r,A,B,\Theta ,\overline{\theta }\right)$  can be rewritten to the following new equivalent form for (m.i.p.) potential:
${\widehat{H}}_{2-pert}\left(r,A,B,\theta ,\overline{\theta }\right)=\frac{1}{2}\left\{\frac{\overline{\theta }}{2m_0}+\theta \left(\frac{A}{r^4}-\frac{B}{2r^3}\right)\right\}\left({\overleftrightarrow{J}}^2-{\overleftrightarrow{L}}^2-{\overleftrightarrow{S}}^2\right)$ (27.1)
${\widehat{H}}_{3-pert}\left(r,A,B,\Theta ,\overline{\theta }\right)=\frac{1}{2}\left\{\frac{\overline{\theta }}{2m_0}+\Theta \left(\frac{A}{r^4}-\frac{B}{2r^3}\right)\right\}\left({\overleftrightarrow{J}}^2-{\overleftrightarrow{L}}^2-{\overleftrightarrow{S}}^2\right)$ (27.2)
To the best of our knowledge, we just replace the coupling spin-orbital $\overleftrightarrow{S}\overleftrightarrow{L}$  by the expression $\frac{1}{2}\left({\overleftrightarrow{J}}^2-{\overleftrightarrow{L}}^2-{\overleftrightarrow{S}}^2\right)$ , in quantum mechanics. The set ( $H_{nc(2-3)-is}\left({\widehat{p}}_i,{\widehat{x}}_i\right)$ , ${\text{J}}^{\text{2}}$ , ${\text{L}}^2$ , ${\text{S}}^2$ and $J_z)$  forms a complete of conserved physics quantities and the eigenvalues of the spin orbital coupling operator are:

$p_{\pm }\left(j=l\pm 1/2,l,s=1/2\right)\equiv \frac{1}{2}\{\begin{array}{l}\left(l+\frac{1}{2}\right)(l+\frac{1}{2}+1)+l(l+1)-\frac{3}{4} \equiv p_{+}\begin{array}{cc}\text{ for}& j= l+\frac{1}{2}\Rightarrow \text{polarization}-\text{up}\end{array}\\ \left(l-\frac{1}{2}\right)(l-\frac{1}{2}+1)+l(l+1)-\frac{3}{4} \equiv p_{-}\begin{array}{cc}\text{ for}& j= l+\frac{1}{2}\Rightarrow \text{polarization}-\text{d}own\end{array}\end{array}$ (27.3)
Which allows us to form a diagonal $\left(2\times 2\right)$  and $\left(3\times 3\right)$ two matrixes, with non null elements are [ ${\left({\widehat{H}}_{so-is}\right)}_{11}$ and ${\left({\widehat{H}}_{so-is}\right)}_{22}$ ] and [ ${\left({\widehat{H}}_{so-is}\right)}_{11}$ , ${\left({\widehat{H}}_{so-is}\right)}_{22}$ , ${\left({\widehat{H}}_{so-is}\right)}_{33}$ ] for (m.i.s.) potential in (NC: 2D-RSP) and (NC: 3D-RSP), respectively, as:

$\begin{array}{l}{\left(H_{so-ip}\right)}_{11}=p_{+}\left(\frac{\overline{\theta }}{2m_0}+\theta \left(\frac{A}{r^4}-\frac{B}{2r^3}\right)\right)\text{if }j=l+\frac{1}{2}\Rightarrow \text{spin -up}\\ {\left(H_{so-ip}\right)}_{22}=p_{-}\left(\frac{\overline{\theta }}{2m_0}+\theta \left(\frac{A}{r^4}-\frac{B}{2r^3}\right)\right)\text{ if }j=l-\frac{1}{2}\Rightarrow \text{spin -down}\end{array}$ (28.1)

and

$\begin{array}{l}{\left({\widehat{H}}_{so-is}\right)}_{11}=p_{+}\left\{\frac{\overline{\theta }}{2m_0}+\Theta \left(\frac{A}{r^4}-\frac{B}{2r^3}\right)\right\}\text{if }j=l+\frac{1}{2}\Rightarrow \text{spin up}\\ {\left({\widehat{H}}_{so-is}\right)}_{22}= p_{-}\left\{\frac{\overline{\theta }}{2m_0}+\Theta \left(\frac{A}{r^4}-\frac{B}{2r^3}\right)\right\}\text{ if}j = l-\frac{1}{2}\Rightarrow \text{spin down}\\ {\left({\widehat{H}}_{so-is}\right)}_{33}=0\end{array}$ (28.2)

Substituting two equations (26.1) and (26.2) into two equations (20.1) and (20.12), respectively and then, the radial parts of the modified Schrödinger equations, satisfying the following important two equations:

$\frac{{{\displaystyle d^{2}R}}_l\left(\rho \right)}{d{\rho }^2}+\frac{2}{\rho }\frac{{{\displaystyle dR}}_l\left(\rho \right)}{d\rho }+\left(-\frac{1}{4}+\frac{\tau }{\rho }-\frac{2A+l^2}{{\rho }^2}-\left\{\frac{\overline{\theta }}{2m_0}+\theta \left(\frac{A}{r^4}-\frac{B}{2r^3}\right)\right\}\overleftrightarrow{S}\overleftrightarrow{L}\right){{\displaystyle R}}_l\left(r\right) = 0$ (29.1)

and

$\frac{1}{{{\displaystyle r}}^2}\frac{\partial }{\partial r}\left({{\displaystyle r}}^2\frac{\partial }{\partial r}\right){{\displaystyle R}}_l\left(r\right)+\left[2\left(E_{nc3-is}\left(n,l,3,\mathrm{...}\right)-V\left(r\right)\right)-\frac{l\left(l+1\right)}{{{\displaystyle r}}^2}-\left\{\frac{\overline{\theta }}{2m_0}+\Theta \left(\frac{A}{r^4}-\frac{B}{2r^3}\right)\right\}\overleftrightarrow{S}\overleftrightarrow{L}\right]{{\displaystyle R}}_l\left(r\right)=0$ (29.2)

for (m.i.s.) potential in (NC: 2D-RSP) and (NC: 3D-RSP), ii is clearly that the above equations including equations (26.1) and (26.2), the perturbative terms of Hamiltonian operator, which we are subject of discussion in next sub-section.

The exact spin-orbital spectrum for (m.i.s.) potential in both (NC: 2D- RSP) and (NC: 3D- RSP) symmetries for states for one-electron atoms

 In this sub section, we are going to study the modifications to the energy levels ( $E_{nc-per:u}\left(\theta ,\overline{\theta }\right)$ , $E_{nc-per:D}\left(\theta ,\overline{\theta }\right)$ ) and ( $E_{nc-per:u}\left(\Theta ,\overline{\theta }\right)$ , $E_{nc-per:D}\left(\Theta ,\overline{\theta }\right)$ ) for spin up and spin down, respectively, at first order of parameters ( $\theta$ , $\overline{\theta }$ ) and ( $\Theta$ , $\overline{\theta }$ ), for excited states $n^{th}$ , obtained by applying the standard perturbation theory, using eqs. (14.1) (14.2), (27.1) and (27.2) corresponding (NC-2D: RSP) and (NC-3D: RSP), respectively, as:
$\begin{array}{l}{E}_{nc-per:u}\left(\theta ,\overline{\theta }\right)\equiv 2\prod p_{+}{\displaystyle \int R^{*}\left(r\right)\left[\theta \left(\frac{A}{{{\displaystyle r}}^4}-\frac{B}{{{\displaystyle 2r}}^3}\right)+\frac{\overline{\theta }}{2m_0}\right]R\left(r\right)rdr}\begin{array}{cc}Si& j=l+\frac{1}{2}\end{array}\\ E_{nc-per:D}\left(\theta ,\overline{\theta }\right)\equiv 2\prod p_{-}{\displaystyle \int R^{*}\left(r\right)\left[\theta \left(\frac{A}{{{\displaystyle r}}^4}-\frac{B}{{{\displaystyle 2r}}^3}\right)+\frac{\overline{\theta }}{2m_0}\right]R\left(r\right)rdr}\begin{array}{cc}Si& j=l-\frac{1}{2}\end{array}\end{array}$  (30.1)
and

$\begin{array}{l}{E}_{nc-per:u}\left(\Theta ,\overline{\theta }\right)\equiv -\frac{\alpha p_{+}}{{\left(-8E\right)}^{3/2}}{\left(\frac{2B}{n}\right)}^3\frac{\left(n-l-1\right)!}{2n\left(n+l\right)!}{\displaystyle \int {{\displaystyle \rho }}^{2l+2}{{\displaystyle e}}^{-\rho }{\left[{{\displaystyle L}}_{n-l-1}^{2l+1}\left(\rho \right)\right]}^2\left[\Theta \left(\frac{A\text{'}}{{{\displaystyle \rho }}^4}-\frac{B\text{'}}{{{\displaystyle 2\rho }}^3}\right)+\frac{\overline{\theta }}{2m_0}\right]d\rho }\\ E_{nc-per:D}\left(\Theta ,\overline{\theta }\right)\equiv -\frac{\alpha p_{-}}{{\left(-8E\right)}^{3/2}}{\left(\frac{2B}{n}\right)}^3\frac{\left(n-l-1\right)!}{2n\left(n+l\right)!}{\displaystyle \int {{\displaystyle \rho }}^{2l+2}{{\displaystyle e}}^{-\rho }{\left[{{\displaystyle L}}_{n-l-1}^{2l+1}\left(\rho \right)\right]}^2\left[\Theta \left(\frac{A\text{'}}{{{\displaystyle \rho }}^4}-\frac{B\text{'}}{{{\displaystyle 2\rho }}^3}\right)+\frac{\overline{\theta }}{2m_0}\right]d\rho }\end{array}$ (30.2)

$\begin{array}{l}{E}_{nc-per:u}\left(\theta ,\overline{\theta }\right)\equiv \frac{2\prod p_{+}}{\left(-8E\right)}{\left(\frac{4{\rm B}}{2n-2m+2s_2-1}\right)}^2\left(\frac{\left(n-m-1\right)!}{\left(2n-2m+2s_2-1\right)\left(n-m+2s_2-1\right)!}\right)\left(\theta {\displaystyle \sum _{i=1}^2{T}_{i-2}}+\frac{\overline{\theta }}{2{{\displaystyle m}}_0}{T}_{3-2}\right)\\ E_{nc-per:D}\left(\theta ,\overline{\theta }\right)\equiv \frac{2\prod p_{-}}{\left(-8E\right)}{\left(\frac{4{\rm B}}{2n-2m+2s_2-1}\right)}^2\left(\frac{\left(n-m-1\right)!}{\left(2n-2m+2s_2-1\right)\left(n-m+2s_2-1\right)!}\right)\left(\theta {\displaystyle \sum _{i=1}^2{T}_{i-2}}+\frac{\overline{\theta }}{2{{\displaystyle m}}_0}{T}_{3-2}\right)\end{array}$ (31.1)

and

$\begin{array}{l}{E}_{nc-per:u}\left(\Theta ,\overline{\theta }\right)\equiv -\frac{\alpha p_{+}}{{\left(-8E\right)}^{3/2}}{\left(\frac{2B}{n}\right)}^3\frac{\left(n-l-1\right)!}{2n\left(n+l\right)!}\left(\Theta {\displaystyle \sum _{i=1}^2{T}_{i-3}}+\frac{\overline{\theta }}{2{{\displaystyle m}}_0}{T}_{3-3}\right)\\ E_{nc-per:D}\left(\Theta ,\overline{\theta }\right)\equiv -\frac{\alpha p_{-}}{{\left(-8E\right)}^{3/2}}{\left(\frac{2B}{n}\right)}^3\frac{\left(n-l-1\right)!}{2n\left(n+l\right)!}\left(\Theta {\displaystyle \sum _{i=1}^2{T}_{i-3}}+\frac{\overline{\theta }}{2{{\displaystyle m}}_0}{T}_{3-3}\right)\end{array}$ (32.2)

Where, the 6- terms: ( $T_{i-2}$ , $T_{i-3}$   $i=1,2$ ), $T_{3-2}$  and $T_{3-3}$  are given by:

$\begin{array}{l}{T}_{1-2=A^{\text{'}}}{\displaystyle \underset{0}{\overset{+\infty }{\int }}{e}^{-\rho }}{\rho }^{2S_{2^{-3}}}{\left[L_{n-m-1}^{2S_2}\left(\rho \right)\right]}^{2}d\rho \\ T_{2-2={\frac{B}{2}}^{\text{'}}}{\displaystyle \underset{0}{\overset{+\infty }{\int }}{e}^{-\rho }}{\rho }^{2S_{2^{-2}}}{\left[L_{n-m-1}^{2S_2}\left(\rho \right)\right]}^{2}d\rho \\ T_{3-2=}{\displaystyle \underset{0}{\overset{+\infty }{\int }}{e}^{-\rho }}{\rho }^{2S_{2^{+1}}}{\left[L_{n-m-1}^{2S_2}\left(\rho \right)\right]}^{2}d\rho \end{array}$ (33.1)