\section{Introduction}
\label{sec:intro}
In this last few decades, non-commutative (NC) theories have been extensively studied. This happens because some results in string theory\cite{ST,STA} suggests that space-time coordinates may not commute\cite{snyder}, since it seems to be relevant to the quantization of D-branes in  background magnetic $(B_{\mu\nu})$ fields\cite{ST1,ST1A} and, also because NC space-time coordinates is an alternative mechanism for Lorentz-invariance breaking\cite{JC,JC1,JC2}. Despite of great interest today, the non-commutativity of coordinates is an older idea presented at the beginning of the quantum theory, indeed, some assumptions on its \textit{quantized} differential geometry can be found in Ref.\cite{dirac,dirac1}. Ever since Heisenberg's paper on uncertainty principle, there has been wide hesitancy in consider, simultaneously, classical-valued positions and momenta variables in any meaningful formula expressing quantum behavior, since these are incompatible observables. However, H. Groenewold\cite{HG} and J. Moyal\cite{MOYAL} provided an original technical solution to the problem above, \textit{i.e.}, they developed a special binary operation, the $\ast$-product, which it preserves the classical nature of positions$(x_i)$ and momenta$(p_i)$, but it also permits $(x_i)$ and $(p_i)$ to combine in a way  that it is equivalent to the familiar operator algebra of Hilbert space quantum theory. This original technical solution bases on the introduction of a Wigner phase space distribution function\cite{wigner} with deformed product and brackets (it is not the Poisson brackets), whose they correspond to the quantum commutators brackets. Due to this, essential aspects of quantum mechanics can be given a classical formalism in terms of the $\ast$-product. In the Ref.\cite{flato} the deformation structure in the symplectic space and the uniqueness of this formulation was systematically analyzed and mathematically consolidated\cite{fedosov} and, more recently, in Ref.\cite{kont} the author proved that any finite-dimensional Poisson manifold can be canonically quantized in the sense of deformation quantization. This alternative quantum framework is called \textit{Quantum Mechanics in Phase Space} (QMPS)\cite{QMPS1,QMPS2,QMPS3,QMPS4}. Inspired in these articles, we propose, in this present article, a generalization of the Bopp's shifts based on the NC symplectic induction procedure which it begins with the deformation of the brackets among the fields, embraced by the symplectic variables, which it is identified as being the inverse of the symplectic matrix (two-form tensor). From the deformed symplectic matrix, the transformation rule among the newer symplectic variables, which they describe the deformed system, and the older one, which they describe the original system, it is determined and it is interpreted as being the general Bopp's shifts and, due to this, the newer Lagrangian and Hamiltonian might to embrace the ``\textit{quantum}'' or NC features into them; at this point, the canonical quantized process can be applied, since the new brackets, now, might correspond to the quantum commutators. This will be done in section II. Applications will be explored in section III and, at the end, conclusions will be presented.

\section{The Non-Commutative Symplectic Induction Formalism and the General Bopp's Shifts Matrix Method}
\label{sec:2}
Considering a general system, represented in a $N$-dimensional configuration space $\Gamma(q_i,\dot q_i)$, with $i=1,2,\dots, N$, whose its dynamics is governed by a Lagrangian $L(q_i,\dot q_i,t)$, which in the reduced first-order description is given by
\be
\label{0000}
L=A^{(0)}_{\xi_\alpha}\cdot\dot\xi_\alpha-H(\xi_\alpha),~~~~~\alpha=(1,2,\dots ,2N),
\ee
where $\xi_\alpha=(q_1, p_1, q_2, p_2, \dots, q_N, p_N)$,  $A^{(0)}_{\xi_\alpha}$ and $H(\xi_\alpha)$ are, respectively, the symplectic variables, the one-form momenta and the symplectic potential; $p_i$ are the canonical momenta conjugated to $q_i$. At this point, the system is described in the symplectic or phase-space $\Sigma(q_i,p_i)\equiv\Sigma(\xi_\alpha)$, {\it i.e}, $\Gamma(q_i,\dot q_i)\rightarrow \Sigma(\xi_\alpha)$. From the Euler-Lagrange equation of motion, we get:
\ba
\label{0010}
\left(\frac{\partial\,A^{(0)}_{\xi_\beta}}{\partial\xi_\alpha}-\frac{\partial\,A^{(0)}_{\xi_\alpha}}{\partial\xi_\beta}\right)\dot\xi_\beta&=&\frac{\partial\,H}{\partial\xi_\alpha},\nonumber\\
f^{(0)}_{\alpha\beta}\dot\xi_\beta&=&\frac{\partial\,H}{\partial\xi_\alpha}, \\
\dot\xi_\beta&=&{f^{(0)}_{\beta\alpha}}^{-1}\frac{\partial\,H}{\partial\xi_\alpha},\nonumber
\ea
where $f^{(0)}_{\alpha\beta}$ is the ($2N\times2N$) symplectic matrix, namely,
\be
\label{0015}
f^{(0)}_{\alpha\beta}=\frac{\partial\,A^{(0)}_{\xi_\beta}}{\partial\xi_\alpha}-\frac{\partial\,A^{(0)}_{\xi_\alpha}}{\partial\xi_\beta}, \nonumber\\
\ee
\be
\label{0020}
{f^{(0)}}=\left(
\begin{array}{rcccl}
-\varepsilon  ~&~  0           ~&~ \cdots  ~&~ 0             ~&~ 0        \cr
  0           ~&~ -\varepsilon ~&~         ~&~ 0             ~&~ 0        \cr
	\vdots      ~&~              ~&~ \ddots  ~&~ \             ~&~ \vdots   \cr
	0           ~&~  0           ~&~         ~&~ -\varepsilon  ~&~ 0        \cr	
	0           ~&~  0           ~&~ \cdots  ~&~ 0             ~&~ -\varepsilon
\end{array}
\right),
\ee
with
\be
\label{0050}
\varepsilon=
\left(
\begin{array}{rl}
\{q_i\,,\,q_i\} ~&~ \{q_i\,,\,p_i\} \cr
\{p_i\,,\,q_i\} ~&~ \{p_i\,,\,p_i\} 
\end{array}
\right) = 
\left(
\begin{array}{rl}
0 ~&~ +1 \cr
-1 ~&~~~ 0 
\end{array}
\right),
\ee
where $A^{(0)}_{\xi_\alpha}=\frac 12\xi_\beta f^{(0)}_{\beta\alpha}$ and ${f^{(0)}_{\beta\alpha}}^{-1}$ is the inverse of the symplectic matrix. From the Hamilton's equation of motion

\be
\label{0025}
\dot\xi_\alpha=\left\{\xi_\alpha,H\right\}=\left\{\xi_\alpha,\xi_\beta\right\}\frac{\partial{H}}{\partial\xi_\beta},
\ee
so, with Eq.(\ref{0010}), the Poisson brackets among the symplectic variables are identified as being

\be
\label{0030}
\left\{\xi_\alpha,\xi_\beta\right\}={f^{(0)}_{\alpha\beta}}^{-1},
\ee
where

\be
\label{0040}
{f^{(0)}}^{-1}=\left(
\begin{array}{rcccl}
\varepsilon  ~&~  0          ~&~ \cdots ~&~ 0             ~&~ 0       \cr
  0          ~&~ \varepsilon ~&~        ~&~ 0             ~&~ 0       \cr
	\vdots     ~&~             ~&~ \ddots ~&~ \             ~&~ \vdots  \cr
	0          ~&~  0          ~&~        ~&~ \varepsilon   ~&~ 0       \cr	
	0          ~&~  0          ~&~ \cdots ~&~ 0             ~&~ \varepsilon
\end{array}
\right).
\ee
The Poisson brackets between two variables are defined as being
\be
\label{0055}
\{F,G\}=\frac{\partial{F}}{\partial\xi_\alpha}\,{f^{(0)}_{\alpha\beta}}^{-1}\,\frac{\partial{G}}{\partial\xi_\beta},
\ee
and due to the Eq.(\ref{0030}), we get

\be
\label{0057}
\{F,G\}=\frac{\partial{F}}{\partial\xi_\alpha}\,\{\xi_\alpha,\xi_\beta\}\,\frac{\partial{G}}{\partial\xi_\beta}.
\ee
This is a brief review of the symplectic formalism for nonconstrained system presented in \cite{FJ2,FJ3,FJ3,FJ4,FJ5}. In order to enable the symplectic formalism to handle with constrained systems, where the Dirac brackets are identified as being the inverse of symplectic matrix, {\it vide} Ref.\cite{FJ2}. The central concept behind of the symplectic formalism is the basic framework for the development of alternative applications: gauge symmetry induction\cite{ANO1,ANO1a,ANO1b,ANO1c}, Clebsch parametrization\cite{WOCN} and symplectic induction of non-commutativity \cite{ANO2,ANO2a,ANO2b}. The later method will be generalized in the present article and, from this generalization, we will obtain a general version for the Bopp's shifts\cite{bopp1,bopp2} and, consequently, a general version for the $\ast$-product. Realizing that the inverse of symplectic matrix, Eq.(\ref{0040}), encompasses the brackets among the variables, it is possible to proposal a transformation in the inverse of symplectic matrix ${f^{(0)}_{\alpha\beta}}^{-1}\rightarrow {\wt{f}}^{-1}_{\alpha\beta}$ which also changes the original brackets, namely,
\be
\label{0060}
\left\{\tilde{\xi}_\alpha,\tilde{\xi}_\beta\right\}={\wt{f}}^{-1}_{\alpha\beta},
\ee
where $\xi_\alpha\rightarrow \tilde{\xi}_\alpha\equiv R_{\alpha\beta}\xi_\beta$, where $R$ is, for a while, an unknown transformation matrix, and $\,\det(\wt{f}^{-1})\ne 0$. After that, the original system, described in the phase-space $\Sigma(\xi_\alpha)$, it is now described in the following non-commutative phase-space $\Sigma(\tilde{\xi}_\alpha)$. Now, the symplectic matrix can be computed as well, since
\be
\label{0070}
\wt{f}^{-1}_{\alpha\beta}\cdot \wt{f}_{\beta\gamma} = \delta_{\alpha\gamma}.
\ee
Therefore, $ \left({f^{(0)}_{\alpha\beta}}^{-1}\rightarrow \wt{f}^{-1}_{\alpha\beta}\right)\Rightarrow \left({f}^{(0)}_{\alpha\beta}\rightarrow \wt{f}_{\alpha\beta}\right) \Rightarrow \left[\left(\xi_\alpha,A^{(0)}_{\xi_\alpha}\right)\rightarrow \left(\tilde{\xi}_\alpha,A_{\wt\xi_\alpha}\right) \right]$.

In order to generalize the way of how to introduce the non-commutativity into a $N$-dimensional system with $4(N-1)$ non-commutative parameters, we consider the existence of $(N\!-\!1)$ \ $(2\times 2)-$matrices, $M_1,M_2,\dots,\text{ and } M_{N-1}$, where each one embraces $4$ independent non-commutative parameters and with their respective determinants represented by
\be
\label{0080}
\det{M_1}={b}_1^2 ~;~~~\det{M_2}={b}_2^2~;~\dots;~~~\det{M_{N-1}}={b}_{N-1}^2,
\ee
with 
\be
\label{0090}
{b}^2={b}_1^2+{b}_2^2+\cdots+{b}_{N-1}^2,
\ee
representing the sum of all of them. For $\,i=1,2,\dots,{N-1}$, 
\be
M_i = \left(
\begin{array}{cc}
\{\tilde{q}_1~,~\tilde{q}_{i+1}\} ~&~ \{\tilde{q}_1~,~\tilde{p}_{i+1}\} \cr
\{\tilde{p}_1~,~\tilde{q}_{i+1}\} ~&~ \{\tilde{p}_1~,~\tilde{p}_{i+1}\} 
\end{array}
\right).     
\ee
The inverse of symplectic matrix  $\wt{f}^{-1}$ in the non-commutative symplectic space $(\Sigma(\tilde{\xi}_\alpha))$ is
\be\arraycolsep=1.4pt\def\arraystretch{1.3}
\label{0100}
\wt{f}^{-1}=\left(
\begin{array}{cccccc}
\varepsilon ~&~  M_1         ~&~  M_2          ~&~ \cdots  ~&~  M_{N-2}     ~&~  M_{N-1} \cr
-M_1^T      ~&~  \varepsilon ~&~  0            ~&~ \cdots  ~&~  0           ~&~  0       \cr
-M_2^T      ~&~  0           ~&~  \varepsilon  ~&~         ~&~  0           ~&~  0       \cr
\vdots      ~&~  \vdots      ~&~               ~&~ \ddots  ~&~              ~&~  \vdots  \cr
-M_{N-2}^T  ~&~  0           ~&~  0            ~&~         ~&~  \varepsilon ~&~  0       \cr
-M_{N-1}^T  ~&~  0           ~&~  0            ~&~ \cdots  ~&~  0           ~&~  \varepsilon
\end{array}\right).  
\ee
This isn't the most possible general matrix for the inverse of the NC symplectic metric, but such choice leads to compact and manageable results. From Eq.(\ref{0070}) and Eq.(\ref{0100}) the symplectic matrix is computed as being
\be\arraycolsep=1.4pt\def\arraystretch{1.5}
\label{0110}
\wt{f}=\frac{1}{1-{b}^2}\left(
\begin{array}{cccc}
-\varepsilon                      ~~&~~ -\varepsilon{M_1}\,\varepsilon                      ~~&~~ \cdots  ~&~ -\varepsilon{M_{N-1}}\,\varepsilon                  \cr
\varepsilon{M_1}^T\varepsilon     ~~&~~ ({b}^2-{b}_1^2-1)\,\varepsilon                      ~~&~~ \cdots  ~&~ \varepsilon{M_1}^T\varepsilon{M}_{N-1}\,\varepsilon \cr
\vdots                             ~&~  \vdots                                               ~&~  \ddots  ~&~ \vdots                                              \cr
\varepsilon{M_{N-1}}^T\varepsilon ~~&~~ \varepsilon{M_{N-1}}^T\varepsilon{M}_1\,\varepsilon ~~&~~ \cdots  ~&~ ({b}^2-{b}_{N-1}^2-1)\,\varepsilon
\end{array}\right).     
\ee
The ellipsis symbol in the second row and third column represent $\;\varepsilon{M_1}^T\varepsilon{M}_i\,\varepsilon\;$ com $\,i=2, ..., N-2\,$  while, in the third row and fourth column, $\;\varepsilon{M_i}^T\varepsilon{M}_{N-1}\,\varepsilon\;$, and in the corresponding symmetrical position, they are obtained from the antisymmetry property of the matrix $\wt{f}$. In this non-commutative symplectic space the first-order Lagrangian is given by
\be
\label{0120}
\wt{L}=A_{\wt\xi_\alpha}\dot{\wt\xi_\alpha}-H(\wt\xi),
\ee
with 
\be
\label{0130}
A_{\wt\xi_\alpha}=\frac12\,\wt\xi_\beta\,\wt{f}_{\beta\alpha}=-\frac12\,\wt{f}_{\alpha\beta}\,\wt\xi_\beta. 
\ee
Since the transformation $f^{(0)}\rightarrow\wt{f}$ changes only the coefficients $A^{(0)}_\xi$: $A^{(0)}_\xi\rightarrow{A_{\wt\xi}}$, with $A_{\wt\xi}$ given by Eq.~(\ref{0130}), this transformation makes the new array elements $( \,\wt{f}\, )$ to be added in the Lagrangian, preserving the shape of the Hamiltonian , {\it{i.e}},

\be
\label{0140}
L=A^{(0)}_{\xi_\alpha}\,\dot\xi_\alpha-H(\xi) \longrightarrow
\wt{L}=A_{\wt\xi_\alpha}\dot{\wt\xi_\alpha}-H(\wt\xi)~,~~~  \wt{L} \neq  L.  
\ee
The next step is to make a transformation that allows the return to the originals symplectic variables $\xi$, \textit{i.e.}, $\Sigma(\tilde{\xi}_\alpha)\rightarrow\Sigma(\xi_\alpha)$. From Eq.~(\ref{0057}), we have
\be
\label{0150}
\{\wt\xi_\alpha,\wt\xi_\beta\}=\frac{\partial{\wt\xi_\alpha}}{\partial\xi_\kappa}\,\{\xi_\kappa,\xi_\lambda\}\,\frac{\partial\wt\xi_\beta}{\partial\xi_\lambda}. 
\ee
Defining the matrix that performs the transformation of symplectic coordinates, the old to the new one,
\be
\label{0160}
\wt\xi=R\,\xi,     
\ee
which it is the general presentation of the Bopp's shifts, Eq.~(\ref{0150}) renders to
\be
\label{0170}
\{\wt\xi_\alpha,\wt\xi_\beta\}=R_{\alpha\kappa}\,\{\xi_\kappa,\xi_\lambda\}\,R_{\beta\lambda}. 
\ee
Taking in account Eq.~(\ref{0030}) and Eq.~(\ref{0060}), Eq.~(\ref{0170}) is rewritten as
\be
\label{0180}
\wt{f}^{-1}=R\cdot{f^{(0)}}^{-1}\cdot R^T.                                
\ee
In accordance with the appendix A, this expression could be rewritten as 
\be
\label{0190}
R^{T}\cdot \wt{f}\cdot R=f^{(0)},                                          
\ee
and the matrix that satisfies Eq.~(\ref{0180}) or Eq.~(\ref{0190}) is
\be
\label{0200}
R=\sqrt{\wt{f}^{-1}\cdot f^{(0)}}.
\ee
Therefore, Eq.~(\ref{0160}) is rewritten as
\be
\label{0210}
\wt\xi_\alpha=\left(\sqrt{\wt{f}^{-1}\cdot f^{(0)}}\right)_{\alpha\beta}\xi_\beta.
\ee
From the Eq.~(\ref{0020}) and Eq.~(\ref{0100}) the term $\sqrt{\wt{f}^{-1}\cdot f^{(0)}}$ is obtained
\be
\arraycolsep=1.4pt
\def\arraystretch{1.5}
\wt{f}^{-1}f^{(0)}=\left(
\begin{array}{cccc}
    \I       ~~&~~     -M_1\,\varepsilon   ~~&~~  \cdots  ~&~ -M_{N-1}\,\varepsilon \cr
M_1^T\varepsilon  ~~&~~          \I        ~~&~~  \cdots  ~&~      0 \cr
\vdots~&~  \vdots~&~ \ddots  ~&~ \vdots\cr
M_{N-1}^T\varepsilon  ~~&~~          0         ~~&~~   \cdots  ~&~    \I
\end{array}\right) ~;
\ee
\be\arraycolsep=1.4pt\def\arraystretch{1.5}
\label{0215}
\sqrt{\wt{f}^{-1}f^{(0)}}=\left(
\begin{array}{cccc}
r\,\I                              ~~&~~ -\frac{1}{2r}{M_1}\,\varepsilon                             ~~&~~ \cdots  ~&~ -\frac{1}{2r}{M_{N-1}}\,\varepsilon                         \cr
\frac{1}{2r}{M_1}^T\varepsilon     ~~&~~ \left(r+\frac{(1-r)}{{b}^2}\,({b}^2-{b}_1^2)\right)\I       ~~&~~ \cdots  ~&~ \frac{(1-r)}{{b}^2}{M_1}^T\varepsilon{M}_{N-1}\,\varepsilon \cr
\vdots                              ~&~  \vdots                                                       ~&~  \ddots  ~&~ \vdots                                                      \cr
\frac{1}{2r}{M_{N-1}}^T\varepsilon ~~&~~ \frac{(1-r)}{{b}^2}{M_{N-1}}^T\varepsilon{M}_1\,\varepsilon ~~&~~ \cdots  ~&~ \left(r+\frac{(1-r)}{{b}^2}\,({b}^2-{b}_{N-1}^2)\right)\I
\end{array}\right), \nonumber       
\ee
with,
\be
\label{0217}
r=\frac12\left(\sqrt{1+{b}}+\sqrt{1-{b}}\right)~.
\ee
Ellipsis in the second row and third column represent $\;\frac{(1-r)}{{b}^2}\,{M_1}^T\varepsilon{M}_i\,\varepsilon\;$ com $\,i=2, ..., {N-2}$, and in the corresponding symmetrical position, they  represent $\;\frac{(1-r)}{{b}^2}\,{M_i}^T\varepsilon{M}_1\,\varepsilon$. In the third row and fourth column the ellipsis represent $\;\frac{(1-r)}{{b}^2}\,{M_i}^T\varepsilon{M}_{N-1}\,\varepsilon\;$, while in the corresponding symmetrical position, they are $\;\frac{(1-r)}{{b}^2}\,{M_{N-1}}^T\varepsilon{M}_i\,\varepsilon$.  
The transformation above induces a transformation in the one-form momenta $A_{\wt\xi_\alpha}$, Eq.(\ref{0130}),
\be
\label{0220}
A_{\wt\xi_\alpha}=-\frac12\,\wt{f}_{\alpha\beta}\wt\xi_\beta=-\frac12\,\left(\wt{f}\cdot \sqrt{\wt{f}^{-1}\cdot f^{(0)}}\right)_{\alpha\beta}\xi_\beta.  
\ee
From Eq.(\ref{0110}) and Eq.(\ref{0215}) the term $-\frac12\cdot\wt{f}\cdot \sqrt{\wt{f}^{-1}\cdot f^{(0)}}$ is straightforward computed and it is given by 
\be
\arraycolsep=1.0 pt\def\arraystretch{1.5}
-\frac12\cdot\wt{f}\cdot\sqrt{\wt{f}^{-1}\cdot f^{(0)}}=\frac{\Gamma}{2}\cdot \Lambda,
\ee
with $\,\Gamma=\frac{1}{\sqrt{1-{b}^2}}$ and $\Lambda$ as being
\be\arraycolsep=1.4pt\def\arraystretch{1.5}
\left(
\begin{array}{cccc}
r\,\varepsilon                                    ~~&~~ \frac{1}{2r}\,\varepsilon{M_1}\,\varepsilon                                                      ~~&~~ \cdots   ~&~ \frac{1}{2r}\,\varepsilon{M_{N-1}}\,\varepsilon \cr
-\frac{1}{2r}\,\varepsilon{M_1}^T\varepsilon      ~~&~~ \left(r-\frac{r-\sqrt{1-{b}^2}}{{b}^2}\,({b}^2-{b}_1^2)\right)\varepsilon                         ~~&~~ \cdots  ~&~ -\left(\frac{r-\sqrt{1-{b}^2}}{{b}^2}\right)\varepsilon{M}_1^T\varepsilon\,M_{N-1}\,\varepsilon \cr
\vdots                                             ~&~\vdots                                                                                               ~&~ \ddots   ~&~ \vdots \cr
-\frac{1}{2r}\,\varepsilon{M_{N-1}}^T\varepsilon  ~~&~~ -\left(\frac{r-\sqrt{1-{b}^2}}{{b}^2}\right)\varepsilon{M}_{N-1}^T\varepsilon\,M_1\,\varepsilon   ~~&~~ \cdots  ~&~ \left(r-\frac{r-\sqrt{1-{b}^2}}{{b}^2}\,({b}^2-{b}_{N-1}^2)\right)\varepsilon  
\end{array}\right).\nonumber
\ee
Ellipsis in the first row and third column mean $\frac{1}{2r}\,\varepsilon{M_i}\,\varepsilon$ with $\,i=2, ..., N-2\,$;  in the second row and third column the ellipsis represent  $\;-\left(\frac{r-\sqrt{1-{b}^2}}{{b}^2}\right)\varepsilon{M}_1^T\varepsilon\,M_i\,\varepsilon\;$ with $\,i=2, ..., N-2\,$; in third row e fourth column, ellipsis represent  $\;-\left(\frac{r-\sqrt{1-{b}^2}}{{b}^2}\right)\varepsilon{M}_i^T\varepsilon\,M_{N-1}\,\varepsilon\;$ with $\,i=2, ..., N-2\,$.  

From the Eq.(\ref{0220}) the product $A_{\wt\xi_i}\,\dot{\wt\xi_i}$ becomes
\be
\label{0230}
A_{\wt\xi_\alpha}\,\dot{\wt\xi_\alpha}=\frac12\,\wt\xi_\beta\,\wt{f}_{\beta\alpha}\,\dot{\wt\xi_\alpha}=\frac12\,{\wt\xi}^{\;T}\cdot \wt{f}\cdot \,\dot{\wt\xi}=\frac12\,\xi^T\cdot R^T\cdot \wt{f}\cdot \,R\cdot \,\dot{\xi}=\frac12\,\xi^T\cdot f^{(0)}\cdot \,\dot{\xi}=A^{(0)}_{\xi_\alpha}\,\dot\xi_\alpha,  
\ee
where the Eq.~(\ref{0190}) was used; it is an interesting result since the kinetic term in the original Lagrangian, Eq.~(\ref{0000}), is restored. The transformation given in Eq.~(\ref{0210}) permits that the newer system could be rewritten in terms of the original symplectic variables which, consequently, it affects both the shape of the Lagrangian as the Hamiltonian, without changing their values, as shown below,
\be
\label{0240}
\wt{L}=A_{\wt\xi_\alpha}\dot{\wt\xi_\alpha}-H(\tilde\xi) =A^{(0)}_{\xi_i}\,\dot\xi_i-\wt{H}(\xi) ~,
\ee
with
\be
\label{0250}
\wt{H}(\xi)=H(\tilde\xi)\,.
\ee
The equations of motion is given by 
\be
\label{0260}
\dot\xi_\alpha={f^{(0)}_{\alpha\beta}}^{-1}\,\frac{\partial{\wt{H}(\xi)}}{\partial{\xi_\beta}}, 
\ee
where NC parameters are embraced by the new Hamiltonian $(\wt{H}(\xi))$, given in Eq.(\ref{0250}).

\section{Applications}
\label{sec:3}

At this section, the NC symplectic induction formalism will be applied in order to introduce NC parameters into a general model. This will be done in order to illustrate and shed some light on the question about the role played by the NC parameter into the model. It is important to notice that for any potential, written in terms of inverse power of $|\vec {x}|$, the non-commutative induction approach could leads this potential to a very complex form, which it obstructs the straightforward calculation of the equation of motion and the implementation of any quantization process. However, the NC symplectic induction formalism gives an alternative way to introduce the NC algebra into the model, even that the potential presents inverse power of $|\vec {x}|$. As the first example, we will consider a two dimensional system where a particle, with mass \textit{m}, suffers the action of a general potential$(V(q_i))$. This particle has its dynamic governed by the following Lagrangian density,
\be
\label{0270A}
{\cal{L}}=\frac {m\cdot\dot q_i^2}{2} - V(q_i), \texttt{ with } i=1,2.
\ee
In agreement with the NC symplectic induction formalism, this Lagrangian should be written in a first-order form, then
\ba
\label{0270}
{\cal{L}}&=&p_i\cdot \dot q_i - \frac {p_i^2}{2m} - V(q_i),\nonumber\\
 &=&p_i\cdot \dot q_i - {\cal{H}},
\ea
where ${\cal{H}}=\frac {p_i^2}{2m} + V(q_i)$ is the Hamiltonian or, in the symplectic language, the symplectic potential. In this context, the variables are embraced by the symplectic variable, namely, $\xi_\alpha=(q_1, p_1,q_2,p_2)$. The respective symplectic matrix is 
\be
\label{0280}
f^{(0)}= \left(
\begin{array}{cccc}
0  & -1 & 0  & 0  \\
+1 & 0  & 0  & 0  \\
0  & 0  & 0  & -1 \\
0  & 0  & +1 & 0 
\end{array}
\right).
\ee
From the inverse of $f^{(0)}$, \textit{i.e.}, ${f^{(0)}}^{-1}$the Poisson brackets among the variables are obtained as
\ba
\label{0290}
\left\{q_i,q_j\right\}=\left\{p_i,p_j\right\}&=& 0,\nonumber\\
\left\{q_i,p_j\right\}=\delta_{ij}.
\ea

Now, we are ready to introduce the NC algebra into the model. To this end, we will propose the following NC brackets among the variables,
\ba
\label{0300}
\left\{\tilde{q}_i,\tilde{q}_j\right\}&=& \varepsilon_{ij}\theta,\nonumber\\
\left\{\tilde{p}_i,\tilde{p}_j\right\}&=& 0,\\
\left\{\tilde{q}_i,\tilde{p}_j\right\}&=&\delta_{ij},\nonumber
\ea
where the original phase-space $\Sigma(\xi_\alpha)$ changes to $\Sigma(\tilde{\xi}_\alpha)$,\textit{i.e.}, $\tilde{\xi}_\alpha=(\tilde{q}_1,\tilde{p}_1,\tilde{q}_2,\tilde{p}_2)$, with $\varepsilon_{12}= -\varepsilon_{21}~=~1$. Note that $\theta$ has canonical dimensions of $(\mathit{length})^2$ and, consequently, it introduces an ultra-violet scale to the problem if it is taken to be small. The inverse of symplectic matrix is,
\be
\label{0310}
{\wt{f}}^{-1}= \left(
\begin{array}{cccc}
0       & +1 & \theta & 0  \\
-1      & 0  & 0      & 0  \\
-\theta & 0  & 0      & +1 \\
0       & 0  & -1     & 0 
\end{array}
\right).
\ee
The respective symplectic matrix is
\be
\label{0320}
\wt{f}= \left(
\begin{array}{cccc}
0  & -1      & 0  & 0      \\
+1 & 0       & 0  & \theta \\
0  & 0       & 0  & -1     \\
0  & -\theta & +1 & 0 
\end{array}
\right).
\ee
Therefore, the NC first-order Lagrangian density, is written as
\ba
\label{0330}
\wt{\cal{L}}&=& A_{\tilde{q}_i}\cdot \dot{\tilde{q}}_i +A_{\tilde{p}_i}\cdot \dot {\tilde{p}}_i- \frac {{\tilde{p}_i}^2}{2m} - V(\tilde{q}_i),\nonumber\\
 &=& A_{\tilde{q}_i}\cdot \dot {\tilde{q}}_i + A_{\tilde{p}_i}\cdot \dot {\tilde{p}}_i - {\cal{H}}(\tilde{\xi}_\alpha).
\ea
From the definition of symplectic matrix, Eq.(\ref{0020}), the one-form canonical momenta, $A_{\tilde{\xi}_\alpha}$, are obtained as
\ba
\label{0340}
A_{\tilde{q}_1}&=& \tilde{p}_1,\nonumber\\
A_{\tilde{q}_2}&=& \tilde{p}_2,\nonumber\\
A_{\tilde{p}_1}&=& -\frac{\theta}{2}\tilde{p}_2,\\
A_{\tilde{p}_2}&=& \frac{\theta}{2}\tilde{p}_1,\nonumber
\ea
which allows us to rewritten the NC first-order Lagrangian, Eq.(\ref{0330}), as
\ba
\label{0350}
\wt{\cal{L}}&=& \tilde{p}_1\cdot \dot{\tilde{q}}_1 + \tilde{p}_2\cdot \dot{\tilde{q}}_2 - \frac{\theta}{2}\tilde{p}_2\cdot \dot {\tilde{p}}_1 + \frac{\theta}{2}\tilde{p}_1\cdot \dot {\tilde{p}}_2 - {\cal{H}}(\tilde{\xi}_\alpha),\nonumber\\
&=& -\left( {\tilde{q}}_1 + \frac{\theta}{2}\tilde{p}_2\right)\cdot\tilde{p}_1 - \left(+{\tilde{q}}_2 - \frac{\theta}{2}\tilde{p}_1\right)\cdot \dot {\tilde{p}}_2 - {\cal{H}}(\tilde{\xi}_\alpha).
\ea
Introducing the following transformation,
\ba
\label{0360}
+ {\tilde{q}}_1 + \frac{\theta}{2}\tilde{p}_2 &\rightarrow& q_1,\nonumber\\
+ {\tilde{q}}_2 - \frac{\theta}{2}\tilde{p}_1 &\rightarrow& q_2,\nonumber\\
{\tilde{p}}_1&\rightarrow& p_1,\\
{\tilde{p}}_2&\rightarrow& p_2,
\ea
the NC first-order Lagrangian renders to
\ba
\label{0371}
\wt{\cal{L}} &=& - q_1\cdot \dot p_1 - q_2\cdot \dot p_2 - \wt{\cal{H}}({\xi}_\alpha),\nonumber\\
&=&  p_1\cdot \dot q_1  +p_2\cdot \dot q_2 - \wt{\cal{H}}(q_1-\frac{\theta}{2}{p}_2\,,\,p_1\,,\,q_2+\frac{\theta}{2}{p}_1\,,\,p_2),\nonumber\\ 
&=&  p_1\cdot \dot q_1  +p_2\cdot \dot q_2 - \left[\frac{p_i^2}{2m}+V(q_1-\frac{\theta}{2}{p}_2\,,\,q_2+\frac{\theta}{2}{p}_1)\right].
\ea
This is the same result obtained when the $\ast$-product is implemented \cite{mezincescu}. 

The result obtained in Eq.(\ref{0371}) could be reproduced employing the general Bopp's shifts matrix method. In order to obtain directly the transformations that lead us to determine the NC Lagrangian and Hamiltonian that describe the system in two-dimensional configuration space, it is necessary to propose the following matrices:
\be
\label{0401}
\varepsilon = \left(
\begin{array}{cc}
~0  ~&~ 1 \cr
-1  ~&~ 0 
\end{array}
\right)~~~~,~~~~ 
M_1 = \left(
\begin{array}{cc}
\{\tilde{q}_1~,~\tilde{q}_2\} ~&~ \{\tilde{q}_1~,~\tilde{p}_2\} \cr
\{\tilde{p}_1~,~\tilde{q}_2\} ~&~ \{\tilde{p}_1~,~\tilde{p}_2\} 
\end{array}
\right).     
\ee
Considering the above $(2\times2)-$matrices, we can set up the $(4\times4)-$matrices as
\be\arraycolsep=1.4pt\def\arraystretch{1.3}
\label{0402}
\wt{f}^{-1}=\left(
\begin{array}{cc}
\varepsilon ~&~  M_1  \cr
-M_1^T ~&~ \varepsilon  
\end{array}\right)~~~~,~~~~  
\wt{f}=\frac{1}{1-{b}^2}\left(
\begin{array}{cc}
-\varepsilon                  ~~&~~  -\varepsilon{M_1}\varepsilon\cr
\varepsilon{M_1}^T\varepsilon ~~&~~  -\varepsilon
\end{array}\right).     
\ee
\be\arraycolsep=1.4pt\def\arraystretch{1.3}
\label{0403}
R=\left(
\begin{array}{cc}
r\,\I                          ~~&~~ -\frac{1}{2r}{M_1}\,\varepsilon  \cr
\frac{1}{2r}{M_1}^T\varepsilon ~~&~~ r\,\I
\end{array}\right)~~~~,~~~~  
S=\frac{1}{2\sqrt{1-{b}^2}}\left(
\begin{array}{cc}
 r\,\varepsilon                               ~~&~~ \frac{1}{2r}\,\varepsilon{M_1}\,\varepsilon\cr
-\frac{1}{2r}\,\varepsilon{M_1}^T\varepsilon  ~~&~~ r\,\varepsilon
\end{array}\right)~.     
\ee
The transformation of symplectic coordinates, the older to the newer one,  
and the transformation of the older symplectic coordinates to the newer kinetic coefficients are performed by
\be
\label{0404}
\wt\xi_\alpha=R_{\alpha\beta}\,\xi_\beta~~~~,~~~~
A_{\wt\xi_\alpha}=S_{\alpha\beta}\,\xi_\beta.  
\ee
We can treat this first case taking $\;{b}^2=\det(M_1)=0\,$, $\;r=1\;$ and filling the 
block matrices given in Eq.(\ref{0402}) and Eq.(\ref{0403}) with the following matrix,
\be
M_1=\left(
\begin{array}{cc}
\theta ~&~ 0 \cr
0      ~&~ 0
\end{array}
\right).     
\ee
Then, using the Eq.~(\ref{0404}), we get the necessary transformations that lead us 
to the Lagrangian given in Eq.(\ref{0371}).


After that, we will explore another possibility. To this end, the following NC brackets among the variables will be proposed,
\ba
\label{0410}
\left\{\tilde{q}_i,\tilde{q}_j\right\}&=& 0 ,\nonumber\\
\left\{\tilde{p}_i,\tilde{p}_j\right\}&=& \varepsilon_{ij}\theta,\\
\left\{\tilde{q}_i,\tilde{p}_j\right\}&=&\delta_{ij}.\nonumber
\ea
Now, the NC parameter $(\theta)$ has canonical dimensions of $(\mathit{mass})^2$ and, consequently, it introduces an infra-red scale if it is taken to be small. Repeating the receipt used before, the NC first-order Lagrangian is obtained as
\be
\label{0420}
\wt{\cal{L}} = (\tilde{p}_1-\frac{\theta}{2}\tilde{q}_2)\cdot \dot{\tilde{q}}_1 + (\tilde{p}_2+\frac{\theta}{2}\tilde{q}_1)\cdot \dot{\tilde{q}}_2 - {\cal{H}}(\tilde{q}_i,\tilde{p}_i).
\ee
Introducing the following transformation,
\ba
\label{0430}
\tilde{q}_i  &\rightarrow& q_i,\nonumber\\
\tilde{p}_1-\frac{\theta}{2}\tilde{q}_2 &\rightarrow&p_1 ,\nonumber\\
\tilde{p}_2+\frac{\theta}{2}\tilde{q}_1 &\rightarrow&p_2 ,\nonumber
\ea
into the NC first-order Lagrangian given in Eq.(\ref{0420}),
\be
\label{0440}
\wt{\cal{L}} = p_i\cdot \dot q_i -\wt{\cal{H}}({q}_1\,,\,p_1+\frac{\theta}{2}q_2\,,\,q_2\,,\,p_2-\frac{\theta}{2}q_1),
\ee
where the NC Hamiltonian is
\ba
\label{0445}
\wt{\cal{H}}&=& \frac{p_i^2}{2m}+ \frac{\theta}{2m}(p_1q_2-p_2q_1)+\frac{\theta^2}{8m}q_i^2+V(q_i),\nonumber\\
\mbox{}&=& \frac{p_i^2}{2m}+ \frac{\theta}{2m}(p_i\varepsilon_{ij}q_j)+\frac{\theta^2}{8m}q_i^2+V(q_i)
\ea
From the Euler-Lagrange equation of motion for the canonical momenta, we get,
\ba
p_1 &=& m\cdot \dot q_1 - \frac{\theta}{2} q_2,\nonumber\\
p_2 &=& m\cdot \dot q_2 + \frac{\theta}{2} q_1.
\ea
Introducing these canonical momenta into the NC first-order Lagrangian, Eq. (\ref{0440}), we get
\be
\label{0450}
\wt{\cal{L}} = \frac{m\dot q_i^2}{2} + \frac{\theta}{2}q_i\varepsilon_{ij}\dot q_j - V(q_i).
\ee
We can access directly the transformations that lead us to the Hamiltonian given in Eq.(\ref{0445}) or the Lagrangian given in Eq.(\ref{0450}) by taking $\;{b}^2=\det(M_1)=0\,$,  $\;r=1$, and filling the block matrices given in Eq.(\ref{0402}) and Eq.(\ref{0403}) with the following matrix,
\be
M_1=\left(
\begin{array}{cc}
0 ~&~ 0 \cr
0 ~&~ \theta
\end{array}
\right)     
\ee
and, in the sequel, using the Eq.(\ref{0404}).
Assuming that this system is a two-dimensional harmonic oscillator, the potential is
\be
\label{0460}
V(q_i) = \frac{m\omega^2}{2}(q_i^2),
\ee
where $\omega$ is the frequency, the NC first-order Lagrangian, Eq. (\ref{0450}), renders to
\ba
\label{0465}
\wt{\cal{L}} &=& \frac{m\dot q_i^2}{2} + \frac{\theta}{2}q_i\varepsilon_{ij}\dot q_j - \frac{m\omega^2}{2}(q_i^2),\nonumber\\
&=& \frac{m\dot q_i^2}{2} + \frac{\theta}{2}q_i\varepsilon_{ij}\dot q_j - \frac{\theta^2}{8m}q_i^2 - \frac{m\omega^+\cdot\omega^-}{2}(q_i^2),
\ea
with
\be
\label{0466}
\omega^\pm=\omega\pm\frac{\theta}{2m}.
\ee
Writing the Lagrangian above in a first-order form, namely,
\be
\wt{\cal{L}} = p_i\cdot \dot{q}_i - \left(\frac{p_i^2}{2m} + \frac{\theta (p_i\varepsilon_{ij}q_j)}{2m} + \frac{\theta^2q_i^2}{4m} + \frac{m\omega^+\cdot\omega^-}{2}(q_i^2)\right),
\ee
the respective Hamiltonian is identified as being
\be
\wt{\cal{H}} = \frac{p_i^2}{2m} + \frac{\theta (p_i\varepsilon_{ij}q_j)}{2m} + \frac{\theta^2q_i^2}{4m} + \frac{m\omega^+\cdot\omega^-}{2}(q_i^2),
\ee
or
\be
\label{0466a}
\wt{\cal{H}} = \frac{p_i^2}{2m} + \frac{\theta (p_i\varepsilon_{ij}q_j)}{2m} + \frac{\theta^2q_i^2}{8m} + \frac{m\omega^2}{2}(q_i^2),
\ee
since the relation given in Eq.(\ref{0466}) was used. 

From the investigation done by Banerjee and Ghosh\cite{BG}, the Hamiltonian of a charged harmonic oscillator in an axially symmetric magnetic field is
\be
\label{0470}
{\cal{H}}=\frac{p_i^2}{2m}+\frac{e\cdot B(t)}{2m\cdot c}(p_i\varepsilon_{ij}q_j) + \frac{e^2\cdot B^2(t)}{8m\cdot c^2}q_i^2 +\frac{m\omega^2}{2}(q_i^2).
\ee
Comparing the Hamiltonian above with the one given in Eq.(\ref{0466a}), we obtain the following identity,
\be
\label{0480}
\theta = \frac{e\cdot B(t)}{c}.
\ee
Therefore, the NC version of two-dimensional harmonic oscillator represents an electron in a very simple Bohr model of the hydrogen atom, where there is a background magnetic field. Due to this, the motion of the electron can be split into components parallel and perpendicular to the magnetic field. In this setup, the NC parameter or the background magnetic field splits the original quantum level into three-levels, where one of the frequency remains unchanged and the others two frequencies changed to
\ba
\label{0490}
\omega^\pm&=& \omega \pm \frac{\theta}{2m},\nonumber\\
&=& \omega \pm \frac{e\cdot B(t)}{2m\cdot c}.
\ea
The Hamiltonian above has a structure very similar to the model of a charged particle in a specified electromagnetic field discussed in \cite{Lewis1,Lewis2}, where the eigenstates are constructed of the invariant operator. In accordance with Banerjee and Ghosh\cite{BG}, the standard Zeeman level is given by
\ba
\label{0500}
E^{\pm}_n &=& (n+\frac 12)\hbar\omega \pm \left[ n+ (j+1/2)\right] \frac{\theta}{m},\nonumber\\
&=& (n+\frac 12)\hbar\omega \pm \left[ n+ (j+1/2)\right] \frac{e\cdot B}{m\cdot c},
\ea
where $B$ is a constant. With this approach, the presence of NC parameter into the system might be interpreted as being the origin of charge property of the particle and the interaction of this one with a background magnetic field. As expected, this changes the energy spectrum of the model, which it was shown above.

After that, we will explore new features from the NC approach discussed here. In order to do this, we will consider a system with two particles where each one presents mass \textit{m}. This system has its dynamics governed by the following Lagrangian density,
\be
\label{0510}
{\cal{L}}=\frac {m\cdot(\dot q_{1,i}^2+ \dot q_{2,i}^2)}{2} - V(q_{1,i},q_{2,i}), \texttt{ with } i=1,2.
\ee
The potential above is, for a while, a general interactive potential. The respective Hamiltonian is
\be
\label{0520}
{\cal{H}}=\frac {(p_{1,i}^2+ p_{2,i}^2)}{2m} + V(q_{1,i},q_{2,i}), \texttt{ with } i=1,2.
\ee
Introducing the NC brackets among the variables,
\ba
\label{0530}
\left\{\tilde{q}_{a,i},\tilde{q}_{b,j}\right\}&=& 0 ,\nonumber\\
\left\{\tilde{p}_{a,i},\tilde{p}_{b,j}\right\}&=& \varepsilon_{ij}\cdot \varepsilon_{ab}\cdot\theta,\\
\left\{\tilde{q}_{a,i},\tilde{p}_{b,j}\right\}&=& \delta_{ij}\cdot \delta_{ab}.\nonumber
\ea
The inverse of symplectic matrix is
\be
\label{0540}
\wt{f}^{-1} = \left(
\begin{array}{cccc}
0             & \delta_{ij}             &  0           & 0                      \\
-\delta_{ij}  & 0                       & 0            & \varepsilon_{ij}\theta \\
0             & 0                       & 0            & \delta_{ij}            \\
0             & \varepsilon_{ij}\theta  & -\delta_{ij} & 0                    
\end{array}
\right).
\ee
The respective symplectic matrix is
\be
\label{0550}
\wt{f} = \left(
\begin{array}{cccc}
0                       & -\delta_{ij} &  \varepsilon_{ij}\theta  & 0            \\
+\delta_{ij}            & 0            & 0                        & 0            \\
\varepsilon_{ij}\theta  & 0            & 0                        & -\delta_{ij} \\
0                       & 0            & +\delta_{ij}             & 0 
\end{array}
\right).
\ee
Repeating the procedure used before, we get
\be
\label{0560}
\tilde{\cal{L}}= (\tilde {p}_{1,i}-\frac{\theta}{2}\varepsilon_{ij}\tilde{q}_{2,j})\cdot \dot{\tilde{q}}_{1,i} + (\tilde {p}_{2,i}-\frac{\theta}{2}\varepsilon_{ij}\tilde{q}_{1,j})\cdot \dot{\tilde{q}}_{2,i} - {\cal{H}}(\tilde{q}_{a,i},\tilde{p}_{a,i}).
\ee
Introducing the transformation given by
\ba
\label{0570}
\tilde{q}_{a,i} &\rightarrow&{q}_{a,i},\nonumber\\
\tilde {p}_{1,i}-\frac{\theta}{2}\varepsilon_{ij}\tilde{q}_{2,j}&\rightarrow&{p}_{1,i},\\
\tilde {p}_{2,i}-\frac{\theta}{2}\varepsilon_{ij}\tilde{q}_{1,j}&\rightarrow&{p}_{2,i},\nonumber
\ea
the NC first-order Lagrangian, Eq.(\ref{0560}), becomes
\be
\label{0580}
\tilde{\cal{L}}= {p}_{a,i}\cdot \dot{{q}}_{a,i} - \tilde{{\cal{H}}}({q}_{1,i},{p}_{1,i}+\varepsilon_{ij}\frac{\theta}{2}q_{2,j},{q}_{2,i},{p}_{2,i}+\varepsilon_{ij}\frac{\theta}{2}q_{1,j}),
\ee
where the NC Hamiltonian is
\ba
\label{0590}
\tilde{{\cal{H}}}({q}_{1,i},{p}_{1,i}+\varepsilon_{ij}\frac{\theta}{2}q_{2,j},{q}_{2,i},{p}_{2,i}+\varepsilon_{ij}\frac{\theta}{2}q_{1,j})&=& \frac{p_{a,i}^2}{2m}+ \frac{\theta}{2m}(p_{1,i}\varepsilon_{ij}q_{2,j}+p_{2,i}\varepsilon_{ij}q_{1,j})\nonumber\\
&+&\frac{\theta^2}{8m}q_{a,i}^2+V(q_{a,i}).
\ea
From Eq.(\ref{0560}) the canonical momenta are computed as
\ba
\label{0600}
{p}_{1,i} &=& m\dot q_{1,i}-\frac{\theta}{2}\varepsilon_{ij}q_{2,j},\nonumber\\
{p}_{2,i} &=& m\dot q_{2,i}-\frac{\theta}{2}\varepsilon_{ij}q_{1,j}.
\ea
Introducing these canonical momenta into Eq.(\ref{0560}), we get
\be
\label{0610}
\tilde{\cal{L}}=\frac{m\dot q_{a,i}^2}{2}- \theta \dot q_{1,i}\varepsilon_{ij}q_{2,j} - V(q_{a,i}).
\ee
The results above could also be obtained using the general Bopp's shifts method matrix, since the transformations that lead us to Hamiltonian given in Eq.(\ref{0590}) and Lagrangian of Eq.(\ref{0610}) are assumed to be
\be
R = \left(
\begin{array}{cccc}
\delta_{ij}                        ~&~ 0           ~&~  0                                 ~&~ 0  \\
0                                  ~&~ \delta_{ij} ~&~ \frac{\theta}{2}\,\varepsilon_{ij} ~&~ 0 \\
0                                  ~&~ 0           ~&~ \delta_{ij}                        ~&~ 0 \\
\frac{\theta}{2}\,\varepsilon_{ij} ~&~ 0           ~&~ 0                                  ~&~ \delta_{ij} 
\end{array}
\right)
~~~~,~~~~
S= \frac12\left(
\begin{array}{cccc}
0                                   ~&~ \delta_{ij}           ~&~ -\frac{\theta}{2}\,\varepsilon_{ij} ~&~ 0  \\
-\delta_{ij}                        ~&~ 0                     ~&~ 0                                   ~&~ 0 \\
-\frac{\theta}{2}\,\varepsilon_{ij} ~&~ 0                     ~&~ 0                                   ~&~ \delta_{ij} \\
 0                                  ~&~ 0                     ~&~ -\delta_{ij}                        ~&~ 0 
\end{array}
\right)
\ee
and, in the sequel, the Eq.(\ref{0404}) should be used.
Since we consider that the two particle behave as two harmonic oscillator, initially independent, we have
\be
\label{0620}
V(q_{a,i})= \frac{m\omega_1^2q_{1,i}^2}{2}+\frac{m\omega_2^2q_{2,i}^2}{2},
\ee
where $\omega_1$ and $\omega_2$ are the frequencies of those harmonic oscillators. Due to this, the NC Lagrangian density given in Eq.(\ref{0610}) becomes
\ba
\label{0630}
\tilde{\cal{L}}&=&\frac{m\dot q_{a,i}^2}{2}- \theta \dot q_{1,i}\varepsilon_{ij}q_{2,j}  - \frac{m\omega_1^2q_{1,i}^2}{2}-\frac{m\omega_2^2q_{2,i}^2}{2}\nonumber\\
&=&\frac{m\dot q_{a,i}^2}{2}- \theta \dot q_{1,i}\varepsilon_{ij}q_{2,j} - \frac{m\omega_1^+\cdot \omega_1^-q_{1,i}^2}{2}-\frac{m\omega_2^+\cdot \omega_2^-q_{2,i}^2}{2} -\frac {\theta^2q_{a,i}^2}{8m},
\ea
where 
\ba
\label{0640}
\omega_1^\pm&=&\omega_1\pm \frac{\theta}{{2}m},\nonumber\\
\omega_2^\pm&=&\omega_2\pm \frac{\theta}{{2}m}.
\ea
Assuming that $\theta = \frac{e\cdot B(t)}{c}$, the NC Lagrangian density given in Eq.(\ref{0630}) is rewritten as
\be
\label{0650}
\tilde{\cal{L}}=\frac{m\dot q_{a,i}^2}{2}- \frac{e\cdot B(t)}{c} \dot q_{1,i}\varepsilon_{ij}q_{2,j} - \frac{m(\omega_1^+\cdot \omega_1^-)q_{1,i}^2}{2}-\frac{m(\omega_2^+\cdot \omega_2^-)q_{2,i}^2}{2}-\frac {e^2 \cdot B(t)^2 \cdot {q}_{a,i}^2}{8mc^2} ,
\ee
while the Hamiltonian given in Eq.(\ref{0590}) renders to
\ba
\wt{\cal{H}}({\xi}_\alpha) &=& \frac{p_{a,i}^2}{2m}+ \frac{\theta}{2m}(p_{1,i}\varepsilon_{ij}q_{2,j}+p_{2,i}\varepsilon_{ij}q_{1,j})+\frac{\theta^2}{8m}q_{a,i}^2+ \frac{m\omega_1^2q_{1,i}^2}{2}+\frac{m\omega_2^2q_{2,i}^2}{2},\nonumber\\
&=& \frac{p_{a,i}^2}{2m}+ \frac{e\cdot B(t)}{2mc}(p_{1,i}\varepsilon_{ij}q_{2,j}+p_{2,i}\varepsilon_{ij}q_{1,j})+\frac {e^2 \cdot B(t)^2 \cdot {q}_{a,i}^2}{4mc^2} \nonumber\\ 
&+&\frac{m(\omega_1^+\cdot \omega_1^-)q_{1,i}^2}{2}+\frac{m(\omega_2^+\cdot \omega_2^-)q_{2,i}^2}{2}.
\ea
Following the procedure developed before the energy spectrum is obtained as being
\ba
\label{0650A}
E^{+-}_{n_1,n_2} &=& (n_1+\frac 12)\hbar\omega_1 + (n_2+ \frac 12)\hbar\omega_2 + \left[ (n_1+(j_1+1/2))-(n_2+(j_2+1/2))\right] 
\frac{e\cdot B}{m\cdot c},\nonumber\\
E^{\pm\pm}_{n_1,n_2} &=& (n_1+\frac 12)\hbar\omega_1 +(n_2+ \frac 12)\hbar\omega_2 \pm \left[ n_1+n_2 +(j_1+j_2+1)\right] 
\frac{e\cdot B}{m\cdot c},\nonumber\\
\ea
This could be interpreted as the two interacting electrons presented into the helium atom. Since this atom has its fundamental state $1s^2$, \textit{i.e.}, $\omega_1=\omega_2$, which it splits the quantum level to four-levels: singlet state: $(\omega^o\cdot\omega^o)$ and three triplet states: $(\omega^-\cdot\omega^-),\,(\omega^-\cdot\omega^+),\,(\omega^+\cdot\omega^+)$. In accordance with the Pauli's exclusion principle, which states that two electrons in an atom can not have the same set of four quantum numbers, the two 1s electrons must have different spins, one spin-up $(\alpha)$ and, the other, spin-down$(\beta)$. A product of an orbital with a spin function is called a spinorbital. For example, electron 1 might occupy a spinorbital which we designate
\be
\label{0660}
\phi(1)=\psi_{1s}(1)\alpha(1) \texttt{ or } \phi(1)=\psi_{1s}(1)\beta(1).
\ee
Spinorbitals can be designated by a single subscript, for example, $\phi_a$ or $\phi_b$, where the subscript stands for a set of four quantum numbers. In a two electron system the occupied spinorbitals $\phi_a$ or $\phi_b$ must be different, meaning that at least one of their four quantum numbers must be unequal. A two-electron spinorbital function of the form
\be
\label{0670}
\Psi(1,2)=\frac{1}{\sqrt{2}}(\phi_a(1)\phi_b(2)-\phi_b(1)\phi_a(2)),
\ee
automatically fulfills the Pauli principle since it vanishes if $a=b$. Moreover, this function associates each electron equally with each orbital, which is consistent with the indistinguishability of identical particles in quantum mechanics. The factor $\frac{1}{\sqrt{2}}$ normalizes the two-particle wave function, assuming that $\phi_a$ and $\phi_b$ are normalized and mutually orthogonal. For the $1s^2$ configuration of helium, the two orbital functions are the same and Eq.(\ref{0670}), for two-electron only, it can be factored into an orbital function times a spin function, consequently it  can be written as
\be
\label{0680}
\Psi(1,2)=\frac{\psi_{1s}(1)\,\psi_{1s}(2)}{\sqrt{2}}\sigma_{nl}(1,2),
\ee
with $\sigma_{nl}(1,2)$ is the two-electron spin function, where the two subscripts are the quantum numbers \textit{n} and \textit{l} for the total electron spin. For the two electron spins in opposing directions (antiparallel) with a total spin angular momentum of zero, we have

\be
\label{0690}
\sigma_{00}(1,2)=\alpha(1)\beta(2)-\beta(1)\alpha(2),
\ee
which it represents a singlet, since the total spin angular momentum is null. However, it is also possible to have both spins in the same state, provided the orbitals are different. There are three possible states for two parallel spins:
\ba
\label{0692}
\sigma_{11}(1,2)&=&\alpha(1)\alpha(2),\nonumber\\
\sigma_{10}(1,2)&=&(\alpha(1)\beta(2)+\beta(1)\alpha(2)),\\
\sigma_{1-1}(1,2)&=&\beta(1)\beta(2).\nonumber
\ea
These make up the triplet spin states, which have the three possible orientations
of a total angular momentum of 1.

The lowest excited state of helium is represented by the electron configuration $1s 2s$ and, as a consequence, $\omega_1\neq\omega_2$. Therefore, it splits the quantum excited level to four-levels: singlet state $(\omega_1^o\cdot\omega_2^o)$ and three triplet states $(\omega_1^-\cdot\omega_2^-),\,(\omega_1^-\cdot\omega_2^+\equiv\omega_1^+\cdot\omega_2^-),\,(\omega_1^+\cdot\omega_2^+)$. When electrons are in different orbitals, their spins can be either parallel or antiparallel. In order that the wave-function satisfies the antisymmetry requirement, the two-electron orbital and spin functions must have opposite behavior under exchange of electron labels. There are four possible states from the 1s 2s configuration: a singlet state
\be
\label{0700}
\Psi^+(1,2)=\frac{1}{\sqrt{2}}\left[\psi_{1s}(1)\psi_{2s}(2)+\psi_{1s}(2)\psi_{2s}(2)\right]\sigma_{00}(1,2),
\ee
and three triplet states
\be
\label{0710}
\Psi^-(1,2)=\frac{1}{\sqrt{2}}\left[\psi_{1s}(1)\psi_{2s}(2)-\psi_{1s}(2)\psi_{2s}(2)\right]\sigma_{nl}(1,2), \texttt{ with } i=1, j=-1,0,+1
\ee
From this investigation, we can argue that the NC parameter plays a role into the two electron system (helium atom) since it can be interpreted as being the electromagnetic interaction between the electrons, which it splits the original quantum level.

At this point, we departure from the Lagrangian density, Eq.(\ref{0510}), and propose a different NC brackets among the variables, namely,
\ba
\label{0720}
\left\{\tilde{q}_{a,i},\tilde{q}_{b,j}\right\}&=& 0 ,\nonumber\\
\left\{\tilde{p}_{a,i},\tilde{p}_{b,j}\right\}&=& \varepsilon_{ij}\cdot \delta_{ab}\cdot(-1)^{\delta_{b1}}\cdot\theta,\\
\left\{\tilde{q}_{a,i},\tilde{p}_{b,j}\right\}&=& \delta_{ij}\cdot \delta_{ab}.\nonumber
\ea
The inverse of symplectic matrix is
\be
\label{0730}
\wt{f}^{-1} = \left(
\begin{array}{cccc}
0            & \delta_{ij}              & 0            & 0                       \\
-\delta_{ij} & -\varepsilon_{ij}\theta  & 0            & 0                       \\
0            & 0                        & 0            & \delta_{ij}             \\
0            & 0                        & -\delta_{ij} & +\varepsilon_{ij}\theta 
\end{array}
\right).
\ee
The respective symplectic matrix is
\be
\label{0740}
\wt{f} = \left(
\begin{array}{cccc}
-\varepsilon_{ij}\theta & -\delta_{ij} & 0                       & 0            \\
+\delta_{ij}            & 0            & 0                       & 0            \\
0                       & 0            & +\varepsilon_{ij}\theta & -\delta_{ij} \\
0                       & 0            & +\delta_{ij}            & 0 
\end{array}
\right).
\ee
Repeating the procedure used before, we get the set of one-form canonical momenta
\ba
\label{0750}
A_{\tilde{q}_{1,i}}&=& \tilde{p}_{1,i}+\frac{\theta}{2}\varepsilon_{ij}\tilde{q}_{1,j},\nonumber\\
A_{\tilde{q}_{2,i}}&=& \tilde{p}_{2,i}-\frac{\theta}{2}\varepsilon_{ij}\tilde{q}_{2,j},\\
A_{\tilde{p}_{1,i}}&=& A_{\tilde{p}_{2,i}}=0.\nonumber
\ea
The corresponding NC first-order Lagrangian is 
\be
\label{0760}
\tilde{\cal{L}}=(\tilde{p}_{1,i}+\frac{\theta}{2}\varepsilon_{ij}\tilde{q}_{1,j})\cdot \dot{\tilde{q}}_{1,i}+ (\tilde{p}_{2,i}-\frac{\theta}{2}\varepsilon_{ij}\tilde{q}_{2,j})\cdot \dot{\tilde{q}}_{2,i} -\left[\frac{\tilde{p}_{a,i}^2}{2m}+V(\tilde{q}_{a,j})\right].
\ee
Introducing the following transformation
\ba
\label{0770}
\tilde{q}_{a,j} &\rightarrow& {q}_{a,j},\nonumber\\
\tilde{p}_{1,i}+\frac{\theta}{2}\varepsilon_{ij}\tilde{q}_{1,j}&\rightarrow& {p}_{1,j},\\
\tilde{p}_{2,i}-\frac{\theta}{2}\varepsilon_{ij}\tilde{q}_{2,j}&\rightarrow& {p}_{2,j},
\ea
into the NC first-order Lagrangian, Eq.(\ref{0760}), we get
\be
\label{0780}
\tilde{\cal{L}}= p_{a,i}\cdot q_{a,i}- \left[\frac{p_{a,i}^2}{2m} +\frac{\theta}{2m}(-p_{1,i}\varepsilon_{ij}q_{1,j}+p_{2,i}\varepsilon_{ij}q_{2,j}) - \frac{\theta^2}{8m}q_{a,i}^2+V({q}_{a,j})\right],
\ee
where the term embraced by the brackets is the Hamiltonian, namely:
\be
\label{0790}
\tilde{\cal{H}}=\frac{p_{a,i}^2}{2m} +\frac{\theta}{2m}(-p_{1,i}\varepsilon_{ij}q_{1,j}+p_{2,i}\varepsilon_{ij}q_{2,j}) - \frac{\theta^2}{8m}q_{a,i}^2+V({q}_{a,j}).
\ee
The canonical momenta can be computed from Eq.(\ref{0780}),
\ba
\label{0800}
p_{1,i}&=& m\cdot \dot{q}_{1,i}+\frac{\theta}{2}\varepsilon_{ij}q_{1,j},\\
p_{2,i}&=& m\cdot \dot{q}_{2,i}-\frac{\theta}{2}\varepsilon_{ij}q_{2,j}.\nonumber
\ea
Introducing the canonical momenta given above into the first-order Lagrangian, Eq.(\ref{0780}), we get
\be
\label{0810}
\tilde{\cal{L}}= \frac{m}{2}\cdot (\dot{q}^2_{1,i}+\dot{q}^2_{2,i}) + \frac{\theta}{2}(\dot{q}_{1,i}\varepsilon_{ij}q_{1,j}-\dot{q}_{2,i}\varepsilon_{ij}q_{2,j})- V({q}_{a,j}).
\ee
Assuming that this model is a pair of an unitary charges of opposite sign in a magnetic field, where the NC parameter is interpreted as being the magnetic field, \textit{i.e.}, $\theta=B$, and doing an \textit{educated guess} for the potential,
\be
\label{0820}
V({q}_{a,j})=\frac{K}{2}(q_{1,j}-q_{2,j})^2 ,
\ee
the Lagrangian density, Eq.(\ref{0810})  renders to
\be
\label{0830}
\tilde{\cal{L}}= \frac{m}{2}\cdot (\dot{q}^2_{1,i}+\dot{q}^2_{2,i}) + \frac{B}{2}(\dot{q}_{1,i}\varepsilon_{ij}q_{1,j}-\dot{q}_{2,i}\varepsilon_{ij}q_{2,j})-\frac{K}{2}(q_{1,j}-q_{2,j})^2 ,
\ee
which it is the result discussed in Ref.\cite{BS}, section I.A.

This last case demands a slightly different treatment for set up  $R$ and $S$ matrices. We can obtain the transformations that lead us to the Hamiltonian of Eq.(\ref{0790}) and Lagrangian of Eq.(\ref{0810}) by means of
\be
R = \left(
\begin{array}{cccc}
\delta_{ij}                         ~&~ 0           ~&~  0                                  ~&~ 0  \\
-\frac{\theta}{2}\,\varepsilon_{ij} ~&~ \delta_{ij} ~&~  0                                  ~&~ 0 \\
0                                   ~&~ 0           ~&~ \delta_{ij}                         ~&~ 0 \\
0                                   ~&~ 0           ~&~ \frac{\theta}{2}\,\varepsilon_{ij}  ~&~ \delta_{ij} 
\end{array}
\right)
~~~~,~~~~
S= \frac12\left(
\begin{array}{cccc}
\frac{\theta}{2}\,\varepsilon_{ij}  ~&~ \delta_{ij}           ~&~ 0                                     ~&~ 0  \\
-\delta_{ij}                        ~&~ 0                     ~&~ 0                                     ~&~ 0 \\
0                                   ~&~ 0                     ~&~ -\frac{\theta}{2}\,\varepsilon_{ij}   ~&~ \delta_{ij} \\
0                                   ~&~ 0                     ~&~ -\delta_{ij}                          ~&~ 0 
\end{array}
\right)~,
\ee
where Eq.(\ref{0404}) also should be used.


\section{Conclusion}
\label{sec:4}
In this article, the NC symplectic induction formalism\cite{ANO2} was used to generalize the Bopp's shifts, now it is called general Bopp's shifts matrix method, which it introduces $4\cdot(N-1)$ NC parameters into a $N$-dimensional system. This result driven us to conclude that $\ast$-product was also generalized, since the usual $\ast$-product induces simple Bopp's shifts\cite{bopp1,bopp2}. Furthermore, we explore the application of both NC symplectic induction formalism and general Bopp's shifts matrix method in some very simple models: the one and two-dimensional harmonic oscillator. From this application it was possible to show that distinct choices for the NC algebra among the symplectic variables generates distinct, in a dynamically point of view, NC systems that could be interpreted as being \textit{quantum} versions for the harmonic oscillator, since the role played by the NC parameters in these systems were responsible to split the quantum states of the hydrogen and helium atom. At this point, the canonical quantized process can be used, since the Poisson brackets, now, might correspond to the quantum commutators, because the NC or \textit{quantum} nature is already incorporated into the system.

