D}$ Hermitian operators, $\Lambda^i_{\mu}$, with $\mu=1,... D$, using the generalized Gell-Mann matrices (GGM) and the identity matrix [@gmm2008]: $$\begin{aligned}
    \Lambda^i_{\mu=1,...\mathcal{N}(\mathcal{N}-1)/2}&=&\frac{1}{\sqrt{2}}(\ket{\beta}\bra{\alpha}+h.c.), \end{aligned}$$ for $\alpha>\beta$, $1\le \alpha$,$\beta\le \mathcal{N}$, $$\begin{aligned}
    \Lambda^i_{\mu=\mathcal{N}(\mathcal{N}-1)/2+1,...\mathcal{N}(\mathcal{N}-1)}&=&\frac{1}{\sqrt{2}i}(\ket{\beta}\bra{\alpha}-h.c.),\end{aligned}$$ for $\alpha>\beta$, $1\le \alpha$,$\beta\le \mathcal{N}$, $$\begin{aligned}
    \Lambda^i_{\mu=\mathcal{N}(\mathcal{N}-1)+1,...\mathcal{N}^2-1}&=& \frac{1}{\sqrt{\alpha(\alpha+1)}} \nonumber \\ &&\times (\sum_{\beta=1}^\alpha\ket{\beta}\bra{\beta}\nonumber \\&&-\alpha\ket{\alpha+1}\bra{\alpha+1}),\end{aligned}$$ for $1\leq\alpha<\mathcal{N}$ $$\begin{aligned}
    \Lambda^i_{\mathcal{D}}&=&\sqrt{\frac{1}{\mathcal{D}}} \mathbb{I}.\end{aligned}$$ With these operators, the local density matrix $\hat\rho_i$, as well as any operator $\hat O^i$ of local observables can be represented as $$\begin{aligned}
      \hat O^i&=&\sum_\mu c^i_\mu \Lambda^i_\mu, ~~~\text{with}\\
      c^i_\mu&=&{\rm Tr}[\Lambda^i_\mu\hat O^i],
    \end{aligned}$$ and $\mu=1,2,...\mathcal{D}$. This allows expressing both one-body and two-body Hamiltonians in the form $ \hat H_i=\sum_\mu c^i_\mu\Lambda^i_\mu$, and $\hat H_{ij}=\sum_{\mu,\nu
D}$ Hermitian operators, $ \Lambda^i_{\mu}$, with $ \mu=1, ... D$, using the generalized Gell - Mann matrices (GGM) and the identity matrix   [ @gmm2008 ]: $ $ \begin{aligned } 
     \Lambda^i_{\mu=1,... \mathcal{N}(\mathcal{N}-1)/2}&=&\frac{1}{\sqrt{2}}(\ket{\beta}\bra{\alpha}+h.c .), \end{aligned}$$ for $ \alpha>\beta$, $ 1\le \alpha$,$\beta\le \mathcal{N}$, $ $ \begin{aligned } 
     \Lambda^i_{\mu=\mathcal{N}(\mathcal{N}-1)/2 + 1,... \mathcal{N}(\mathcal{N}-1)}&=&\frac{1}{\sqrt{2}i}(\ket{\beta}\bra{\alpha}-h.c.),\end{aligned}$$ for $ \alpha>\beta$, $ 1\le \alpha$,$\beta\le \mathcal{N}$, $ $ \begin{aligned } 
     \Lambda^i_{\mu=\mathcal{N}(\mathcal{N}-1)+1,... \mathcal{N}^2 - 1}&= & \frac{1}{\sqrt{\alpha(\alpha+1) } } \nonumber \\ & & \times (\sum_{\beta=1}^\alpha\ket{\beta}\bra{\beta}\nonumber \\&&-\alpha\ket{\alpha+1}\bra{\alpha+1}),\end{aligned}$$ for $ 1\leq\alpha<\mathcal{N}$ $ $ \begin{aligned } 
     \Lambda^i_{\mathcal{D}}&=&\sqrt{\frac{1}{\mathcal{D } } } \mathbb{I}.\end{aligned}$$ With these operators, the local concentration matrix $ \hat\rho_i$, equally well as any hustler $ \hat O^i$ of local observables can be represented as $ $ \begin{aligned } 
       \hat O^i&=&\sum_\mu c^i_\mu \Lambda^i_\mu, ~~~\text{with}\\ 
       c^i_\mu&=&{\rm Tr}[\Lambda^i_\mu\hat O^i ], 
     \end{aligned}$$ and $ \mu=1,2,... \mathcal{D}$. This allows expressing both one - consistency and two - body Hamiltonians in the form $ \hat H_i=\sum_\mu c^i_\mu\Lambda^i_\mu$, and $ \hat H_{ij}=\sum_{\mu,\nu
D}$ Hfrmitian operators, $\Lambdx^i_{\mu}$, with $\mu=1,... D$, osung thx generzlized Gdll-Mann matrices (GGM) and the ieentiuj matrix [@gmm2008]: $$\begin{aligved}
    \Lalbda^i_{\mu=1,...\mqthcel{N}(\mathcal{N}-1)/2}&=&\frac{1}{\sqrt{2}}(\ket{\bebc}\bra{\amiha}+h.c.), \xnd{aligned}$$ for $\slpha>\beta$, $1\ne \alpha$,$\beta\le \mxtkcal{N}$, $$\begin{aligned}
    \Lambda^i_{\mu=\mathcwl{N}(\mathvap{N}-1)/2+1,...\mathcal{N}(\mathsal{N}-1)}&=&\gwac{1}{\ssgt{2}l}(\ket{\beta}\bra{\alpha}-h.c.),\end{aligned}$$ for $\alpha>\btta$, $1\le \alpha$,$\beta\le \mathcal{N}$, $$\begin{aligned}
    \Pambfa^i_{\mu=\mathcal{N}(\mathcwl{N}-1)+1,...\mathcal{N}^2-1}&=& \fras{1}{\wqrt{\alpha(\alpfa+1)}} \nonumbeg \\ &&\times (\suj_{\beta=1}^\alpha\ket{\beta}\bra{\beta}\nonumbef \\&&-\alpka\ket{\alpha+1}\btc{\qlpjd+1}),\end{aligned}$$ for $1\jeq\alpha<\mathgsl{N}$ $$\becin{aligmed}
    \Lambda^i_{\katical{E}}&=&\sqrt{\frac{1}{\mathcal{D}}} \matibb{I}.\end{aligned}$$ With jhese operdtkrs, the local denwiry majrix $\vat\rfi_i$, xs sekl as anj o'erator $\hat O^i$ of locao observables can bt rqitesented as $$\gegin{ajidned}
      \hat O^i&=&\sum_\mu c^i_\mu \Lambda^i_\mu, ~~~\teqt{wifh}\\
      c^i_\mu&=&{\rm Tr}[\Lambda^u_\mu\hat O^i],
    \end{aligngd}$$ and $\mu=1,2,...\mwthcal{D}$. This allows expressing both one-body and tfo-bodb Famnotonixbs in the form $ \hat H_i=\sum_\mu c^i_\mu\Lambda^i_\mu$, and $\hwf N_{im}=\sum_{\mu,\nu
D}$ Hermitian operators, $\Lambda^i_{\mu}$, with $\mu=1,... D$, generalized matrices (GGM) the identity matrix $\alpha>\beta$, \alpha$,$\beta\le \mathcal{N}$, $$\begin{aligned} for $\alpha>\beta$, $1\le \mathcal{N}$, $$\begin{aligned} \Lambda^i_{\mu=\mathcal{N}(\mathcal{N}-1)+1,...\mathcal{N}^2-1}&=& \frac{1}{\sqrt{\alpha(\alpha+1)}} \nonumber \\ (\sum_{\beta=1}^\alpha\ket{\beta}\bra{\beta}\nonumber \\&&-\alpha\ket{\alpha+1}\bra{\alpha+1}),\end{aligned}$$ for $1\leq\alpha<\mathcal{N}$ $$\begin{aligned} \Lambda^i_{\mathcal{D}}&=&\sqrt{\frac{1}{\mathcal{D}}} \mathbb{I}.\end{aligned}$$ With these operators, the local density $\hat\rho_i$, as well as any operator $\hat O^i$ of local observables can be as \hat c^i_\mu ~~~\text{with}\\ c^i_\mu&=&{\rm Tr}[\Lambda^i_\mu\hat O^i], \end{aligned}$$ and $\mu=1,2,...\mathcal{D}$. This allows expressing both one-body and two-body Hamiltonians in form $ \hat H_i=\sum_\mu c^i_\mu\Lambda^i_\mu$, and $\hat H_{ij}=\sum_{\mu,\nu
D}$ Hermitian operators, $\Lambda^I_{\mu}$, with $\mu=1,... D$, Using The GenErAlizEd GeLl-Mann matrices (ggM) anD the identity matrix [@gmm2008]: $$\bEgin{aLiGNed}
    \LAMbDa^i_{\mu=1,...\Mathcal{n}(\MaTHCal{n}-1)/2}&=&\fRaC{1}{\sqRt{2}}(\KEt{\Beta}\bRa{\aLpha}+h.c.), \eNd{aligned}$$ fOr $\aLpHa>\beta$, $1\le \alphA$,$\BeTa\le \mathcaL{N}$, $$\bEgin{aligned}
    \LAmbDa^i_{\mu=\mAtHcaL{n}(\mathCal{n}-1)/2+1,...\mathCal{N}(\maTHcal{N}-1)}&=&\fRac{1}{\sqrt{2}i}(\kEt{\BEta}\bra{\ALpha}-h.c.),\eND{AlIgneD}$$ for $\alpha>\beta$, $1\le \alPHa$,$\BEta\le \mathcal{N}$, $$\bEgin{alIgNEd}
    \lAMbdA^i_{\mU=\mathcal{N}(\mAtHcal{N}-1)+1,...\MAthcal{N}^2-1}&=& \FRaC{1}{\SQRt{\aLPha(\alpha+1)}} \nonumBer \\ &&\times (\sum_{\BEta=1}^\Alpha\kEt{\BetA}\Bra{\betA}\nonuMbER \\&&-\alPha\ket{\alpha+1}\Bra{\aLpha+1}),\end{alIgned}$$ fOR $1\leq\alpHA<\mathcaL{N}$ $$\begiN{alIgnEd}
    \LaMBdA^i_{\MatHcAL{D}}&=&\sQRt{\FraC{1}{\MatHcal{D}}} \matHbB{I}.\End{alIgneD}$$ wITH theSe oPeraTors, tHe local densitY maTrix $\HAt\rHo_i$, as Well aS any OpEratoR $\hat O^i$ Of locAl Observables can bE repResented aS $$\beGiN{alIgNed}
      \haT o^i&=&\sum_\mU c^i_\Mu \LAmbda^i_\mU, ~~~\text{wiTH}\\
      c^i_\Mu&=&{\RM tR}[\LAmbda^i_\mu\hat O^i],
    \end{alIgNED}$$ aNd $\mu=1,2,...\mathCal{D}$. ThIS aLlOWs expresSiNg bOth oNE-Body aNd twO-BoDy HamiltOnians IN tHe Form $ \hat h_i=\Sum_\mu c^I_\mU\LaMbdA^i_\mu$, aND $\hat h_{ij}=\sum_{\Mu,\nu
D}$ Hermitian operators, $ \Lambda^i_ {\mu} $,wit h$\mu =1,. .. D$, using t h e ge neralized Gell-Mann ma trice s( GGM) an d the  identi t ym a tri x[@ gmm 20 0 8] : $$\ beg in{alig ned}
    \ Lam bd a^i_{\mu=1,. . .\ mathcal{N} (\m athcal{N}-1) /2} &=&\fr ac {1} { \sqrt {2} }(\ke t{\bet a }\bra{ \alpha}+h .c . ), \en d {aligne d } $$  for  $\alpha>\beta$,$ 1\ l e \alpha$,$\be ta\le\m a th c a l{N }$,  $$\begin{ al igned } 
    \L a mb d a ^ i_{ \ mu=\mathcal{N }(\mathcal{ N }-1 )/2+1, .. .\m a thcal{ N}(\m at h cal {N}-1)}&=&\ frac {1}{\sqrt {2}i}( \ ket{\be t a}\bra{ \alpha }-h .c. ),\e n d{ al ign ed } $$f or  $\ a lph a>\beta$ ,$1 \le \ alph a $ , $ \bet a\l e \m athca l{N}$, $$\beg in{ alig n ed} 
   \Lamb da^i _{ \mu=\ mathca l{N}( \m athcal{N}-1)+1, ...\ mathcal{N }^2 -1 }&= &\frac { 1}{\sq rt{ \al pha(\al pha+1)} }  \n on u m b er  \\ &&\times (\sum _{ \ b et a=1}^\al pha\ke t {\ be t a}\bra{\ be ta} \non u m ber \ \&&- \ al pha\ket{ \alpha + 1} \b ra{\alp ha +1}),\ en d{a lig ned}$ $  for  $1\le q\alpha< \math c al{N}$ $$\begi n {aligned}
  \L a m bd a ^i_{ \ma thcal{D}}&= &\sq r t{\f rac{ 1 }{ \ma t hcal{ D}}}\m a th b b{I}.\end{aligned}$ $With t heseoperators, th e local de n s i ty matri x $\ h at \ rho_i$, as wel l asany operat o r $\hatO^i$of local  observab l e s can be  re pre sen ted a s$$\begin{alig n e d}
    \hat  O^ i&=&\su m_\ muc^i _\m u\Lambda^i _\mu, ~~ ~\ te xt {w ith }\\
    c^i_ \m u&= &{ \rm  Tr}[ \ Lambda ^i_\m u\ha tO^ i ],    \en d {a l i gned }$ $and$\m u= 1,2,. ..\m a thc al{D}$.  This all ows expr es si ng both  one-body and  t wo-body Ha mi lto niansi n  the for m $ \hat H_i=\sum_\mu c ^ i_\mu\L amb da^i_ \mu$ , and $\h atH_{ij} =\s u m_{\mu ,\nu
D}$ Hermitian_operators, $\Lambda^i_{\mu}$,_with $\mu=1,... D$, using_the generalized_Gell-Mann_matrices (GGM)_and_the identity matrix [@gmm2008]:_$$\begin{aligned}
  _ \Lambda^i_{\mu=1,...\mathcal{N}(\mathcal{N}-1)/2}&=&\frac{1}{\sqrt{2}}(\ket{\beta}\bra{\alpha}+h.c.), \end{aligned}$$ for_$\alpha>\beta$, $1\le \alpha$,$\beta\le_\mathcal{N}$,_$$\begin{aligned}
    \Lambda^i_{\mu=\mathcal{N}(\mathcal{N}-1)/2+1,...\mathcal{N}(\mathcal{N}-1)}&=&\frac{1}{\sqrt{2}i}(\ket{\beta}\bra{\alpha}-h.c.),\end{aligned}$$ for $\alpha>\beta$, $1\le \alpha$,$\beta\le \mathcal{N}$, $$\begin{aligned}
   _\Lambda^i_{\mu=\mathcal{N}(\mathcal{N}-1)+1,...\mathcal{N}^2-1}&=&_\frac{1}{\sqrt{\alpha(\alpha+1)}} \nonumber_\\_&&\times_(\sum_{\beta=1}^\alpha\ket{\beta}\bra{\beta}\nonumber \\&&-\alpha\ket{\alpha+1}\bra{\alpha+1}),\end{aligned}$$ for $1\leq\alpha<\mathcal{N}$ $$\begin{aligned}
_   \Lambda^i_{\mathcal{D}}&=&\sqrt{\frac{1}{\mathcal{D}}} \mathbb{I}.\end{aligned}$$_With these_operators, the local density matrix $\hat\rho_i$, as well_as_any operator $\hat_O^i$ of local observables can be represented as $$\begin{aligned}
_     \hat_O^i&=&\sum_\mu c^i_\mu \Lambda^i_\mu,_~~~\text{with}\\
__   _c^i_\mu&=&{\rm Tr}[\Lambda^i_\mu\hat O^i],
   _\end{aligned}$$ and $\mu=1,2,...\mathcal{D}$. This allows expressing_both one-body and two-body Hamiltonians in the_form $ \hat H_i=\sum_\mu c^i_\mu\Lambda^i_\mu$, and_$\hat H_{ij}=\sum_{\mu,\nu