delta}^{\sigma'}(j)}
{|i-j|}\,,$$ $$\label{}
C_{\alpha \sigma,\beta \sigma'}
^{\gamma \sigma,\delta \sigma'}=\sum_{i}
{\phi_{\alpha}^{\sigma*}(i)\phi_{\gamma}^{\sigma}(i)
\phi_{\beta}^{\sigma'*}(i)\phi_{\delta}^{\sigma'}(i)}\,.$$ Omitting the off-diagonal matrix elements ($(\alpha,\beta)\not=
(\gamma,\delta)$) in Eq.(7) yields the Hartree-Fock approximation. The reason for doing this will be discussed subsequently.

The basic parameters of our theoretical framework are the electron number $n$, the long-range Coulomb interaction $U$, the Hubbard interaction strength $U_{Hub}$, and the disorder amplitude $W$. Equations (3-7) constitute a set of self-consistent equations to yield the eigenstates $\phi_{\alpha}^{\sigma},\epsilon_{\alpha}^{\sigma}$ and the energy $E_{G}$ for a given set of spin occupation numbers ($n_{\uparrow},n_{\downarrow}$) with $n_{\uparrow}+n_{\downarrow}=n$. At every step of the numerical calculation, the matrix elements \[described by Eqs.(4-6)\] of $H^{eff}$ will be calculated for the spin occupation numbers that give the minimum energy in the previous step. The first Slater determinant of Eq.(3) will be formed by noninteracting eigenfunctions of Hamiltonian \[Eq.(1)\] with $U=U_{Hub}=0$. Varying the spin occupation numbers so that $M_{z}\in [0,1]$ we calculate $E_{G}(M_{z})$ whose minimum value gives the spin magnetization $M_{z}$ of the system. The stability of the self-consistent solutions is very sensitive at larger Coulomb interaction strengths, forcing us to keep a relatively small values of the interaction parameter $U$. Because of the spin rotational symmetry of the Hamiltonian, one limits the numerical calculations only to positive spin magnetization ($n_{\uparrow}\geq n_{\downarrow}$).

Increasing the amplitude of disorder the spectrum of noninteracting Hamiltonian $H_{0}=t\sum_{<i,j>,\sigma}c_{i\sigma}^{\dagger}c_{j\sigma}
+\sum_{i,\sigma}\omega_{i}n_{i,\sigma}$ evolves from weakly to strongly localized regime while the eigenvalue statistics exhibits a crossover between Wigner and Poisson distributions [@zh
delta}^{\sigma'}(j) } 
 { |i - j|}\,,$$ $ $ \label { } 
 C_{\alpha \sigma,\beta \sigma' } 
 ^{\gamma \sigma,\delta \sigma'}=\sum_{i } 
 { \phi_{\alpha}^{\sigma*}(i)\phi_{\gamma}^{\sigma}(i) 
 \phi_{\beta}^{\sigma'*}(i)\phi_{\delta}^{\sigma'}(i)}\,.$$ Omitting the off - diagonal matrix elements ($ (\alpha,\beta)\not= 
 (\gamma,\delta)$) in Eq.(7) yields the Hartree - Fock estimate. The rationality for doing this will be discussed subsequently. 

 The basic argument of our theoretical model are the electron number $ n$, the long - range Coulomb interaction $ U$, the Hubbard interaction intensity $ U_{Hub}$, and the disorder amplitude $ W$. Equations (3 - 7) constitute a stage set of self - consistent equations to move over the eigenstates $ \phi_{\alpha}^{\sigma},\epsilon_{\alpha}^{\sigma}$ and the energy $ E_{G}$ for a given bent of tailspin occupation numbers ($ n_{\uparrow},n_{\downarrow}$) with $ n_{\uparrow}+n_{\downarrow}=n$. At every step of the numeric calculation, the matrix elements \[described by Eqs.(4 - 6)\ ] of $ H^{eff}$ will be calculated for the tailspin occupation numbers that give the minimum department of energy in the previous step. The first Slater determinant of Eq.(3) will be form by noninteracting eigenfunctions of Hamiltonian \[Eq.(1)\ ] with $ U = U_{Hub}=0$. Varying the spin occupation numbers so that $ M_{z}\in [ 0,1]$ we calculate $ E_{G}(M_{z})$ whose minimal value gives the spin magnetization $ M_{z}$ of the system. The stability of the self - consistent solutions is very sensible at bigger Coulomb interaction strengths, forcing us to keep a relatively small values of the interaction parameter $ U$. Because of the spin rotational symmetry of the Hamiltonian, one limits the numeric calculations only to incontrovertible tailspin magnetization ($ n_{\uparrow}\geq n_{\downarrow}$). 

 Increasing the amplitude of disorder the spectrum of noninteracting Hamiltonian $ H_{0}=t\sum_{<i, j>,\sigma}c_{i\sigma}^{\dagger}c_{j\sigma } 
 + \sum_{i,\sigma}\omega_{i}n_{i,\sigma}$ evolves from weakly to strongly localized regime while the eigenvalue statistics show a crossover between Wigner and Poisson distributions [ @zh
delha}^{\sigma'}(j)}
{|i-j|}\,,$$ $$\label{}
C_{\alpha \slgma,\beta \sigma'}
^{\gamma \sigka,\deltz \sigma'}=\sjm_{i}
{\phi_{\alpha}^{\sigma*}(i)\phi_{\gamma}^{\sigla}(u)
\phi_{\btna}^{\sigma'*}(i)\phi_{\delta}^{\sigma'}(k)}\,.$$ Omittinh the ofd-diajonal matrix elements ($(\alika,\betz)\kot=
(\gakna,\delta)$) in Eq.(7) iields the Hdrtree-Fock appsobilation. The reason for doing this wijl be dosfussed subsequgntly.

Uhe basjb karameters of our theoretical frajework ere the electrom number $n$, the long-range Clulolb interaction $U$, tje Hubbard untewqction streneth $U_{Hub}$, and the disorser amplitude $W$. Equations (3-7) consgitutz a set of weof-cltsistent eqnationf to yield tmv eigenvtates $\lhi_{\alpha}^{\sigma},\eisiloi_{\alpya}^{\sigma}$ and the energb $E_{G}$ for a given set of spin mceupation numbers ($n_{\upaeriw},n_{\dofnarsow}$) duth $n_{\ukarcow}+h_{\downagrox}=n$. At every step of thw numerical calculauiog, the matrix emementf \[qescribed by Eqs.(4-6)\] of $H^{eff}$ will be calculdtes for the spin occupatiin numbers that give jhe minimui energy in the previous step. The first Slater dedermiiavt if Dw.(3) aill be formed by noninteracting eigenfunctioga pf Hamiltonian \[Ez.(1)\] with $U=U_{Hun}=0$. Gatiing the spin uccupacjoh numbers so that $L_{z}\in [0,1]$ wg calcylate $E_{G}(M_{e})$ whoxe minimum value gives the wpin magnetidatiin $M_{z}$ of the systel. The stabiuity of tne self-consistent solutnons ia very sensltive at mxrger Coulomb ingerscdion strtvgths, forcing us eo keep a relctively rmalk valuqs of the lnteraction parameter $U$. Hecaude of the spln rotational symmetry of the Hekiltonian, ong lhmins the nuiericsl calculatiogs only to posntive spnn magvetization ($n_{\uparrmw}\geq n_{\downwrrow}$).

Increasith the amplitnde of difordwr tye specgfum of nonintetacting Hciiltinian $H_{0}=t\sum_{<i,j>,\sigma}g_{i\sigoz}^{\dagger}c_{j\sigma}
+\sbo_{u,\sigma}\omega_{i}n_{i,\sogmx}$ edopvxs frjk weakly to vtrovglh locauized rtyimc wfile the eigenvalue stathstids exhibits a crosxoyer betwegn Wigner and Poisson cistributions [@zh
delta}^{\sigma'}(j)} {|i-j|}\,,$$ $$\label{} C_{\alpha \sigma,\beta \sigma'} ^{\gamma {\phi_{\alpha}^{\sigma*}(i)\phi_{\gamma}^{\sigma}(i) Omitting the matrix elements ($(\alpha,\beta)\not= Hartree-Fock The reason for this will be subsequently. The basic parameters of our framework are the electron number $n$, the long-range Coulomb interaction $U$, the Hubbard strength $U_{Hub}$, and the disorder amplitude $W$. Equations (3-7) constitute a set of equations yield eigenstates and the energy $E_{G}$ for a given set of spin occupation numbers ($n_{\uparrow},n_{\downarrow}$) with $n_{\uparrow}+n_{\downarrow}=n$. At step of the numerical calculation, the matrix elements by Eqs.(4-6)\] of $H^{eff}$ be calculated for the spin numbers give the energy the step. The first determinant of Eq.(3) will be formed by noninteracting eigenfunctions of Hamiltonian \[Eq.(1)\] with $U=U_{Hub}=0$. Varying the spin numbers so [0,1]$ we $E_{G}(M_{z})$ minimum gives the spin of the system. The stability of is very sensitive at larger Coulomb interaction strengths, us to a relatively small values of the parameter $U$. Because of the spin rotational symmetry the Hamiltonian, one limits the numerical calculations only to positive spin magnetization ($n_{\uparrow}\geq n_{\downarrow}$). Increasing of disorder the spectrum noninteracting Hamiltonian $H_{0}=t\sum_{<i,j>,\sigma}c_{i\sigma}^{\dagger}c_{j\sigma} evolves weakly strongly regime while eigenvalue statistics exhibits a crossover between Wigner and Poisson distributions [@zh
delta}^{\sigma'}(j)}
{|i-j|}\,,$$ $$\label{}
C_{\alpha \sIgma,\beta \siGma'}
^{\gaMma \SigMa,\DeltA \sigMa'}=\sum_{i}
{\phi_{\alpha}^{\SIgma*}(I)\phi_{\gamma}^{\sigma}(i)
\phi_{\beta}^{\Sigma'*}(I)\pHI_{\delTA}^{\sIgma'}(i)}\,.$$ omittinG ThE OFf-dIaGoNal MaTRiX elemEntS ($(\alpha,\bEta)\not=
(\gammA,\deLtA)$) in Eq.(7) yields tHE HArtree-Fock AppRoximation. ThE reAson foR dOinG This wIll Be disCussed SUbsequEntly.

The bAsIC paramETers of oUR ThEoreTical framework are THe ELectron number $n$, The lonG-rANgE cOulOmb InteractioN $U$, The HuBBard intERaCTIOn sTRength $U_{Hub}$, and The disorder AMplItude $W$. eqUatIOns (3-7) conStituTe A Set Of self-consiStenT equationS to yieLD the eigENstates $\Phi_{\alpHa}^{\sIgmA},\epsILoN_{\aLphA}^{\sIGma}$ ANd The ENerGy $E_{G}$ for a GiVeN set oF spiN OCCUpatIon NumbErs ($n_{\uParrow},n_{\downarRow}$) With $N_{\UpaRrow}+n_{\DownaRrow}=N$. AT everY step oF the nUmErical calculatiOn, thE matrix elEmeNtS \[deScRibed BY Eqs.(4-6)\] of $h^{efF}$ wiLl be calCulated FOr tHe SPIN oCcupation numbers thAt GIVe The minimUm enerGY iN tHE previouS sTep. the fIRSt SlaTer dETeRminant oF Eq.(3) wilL Be FoRmed by nOnInteraCtIng EigEnfunCTionS of HamIltonian \[eq.(1)\] witH $u=U_{Hub}=0$. Varying thE Spin occupatioN NuMBErS So thAt $M_{Z}\in [0,1]$ we calculAte $E_{g}(m_{z})$ whOse mINiMum VAlue gIves tHe SPiN Magnetization $M_{z}$ of thE sYstem. THe staBility of the seLf-consisteNT SOlutions Is veRY sENsitive at largeR CoulOmb interacTIon strenGths, fOrcing us To keep a reLATively smAll ValUes Of tHE InTeraction paraMETer $U$. beCause of The Spin rotAtiOnaL syMmeTrY of the HamIltonian, OnE lImItS thE numeRIcal calcUlAtiOnS onLy to pOSitive Spin mAgneTiZaTIon ($N_{\uparroW}\GeQ N_{\DownArRoW}$).

IncReaSiNg the AmplITudE of disoRder the spEctRUm of NoNiNteractIng HamiltoniaN $H_{0}=T\sum_{<i,j>,\sigmA}c_{I\siGma}^{\dagGER}c_{j\sigma}
+\Sum_{i,\sigma}\omega_{i}n_{i,\sigma}$ eVOlves frOm wEakly To stRongly locAliZed regIme WHile thE eigenValue StAtiSTIcs exHIBiTs a CrOssover betWEEn WIgner AnD PoiSson disTributions [@zh
delta}^{\sigma'}(j)}
{|i-j |}\,,$$ $$ \labe l{} 
C_ {\ alph a \s igma,\beta \si g ma'} 
^{\gamma \sigma,\delt a \si gm a '}=\ s um _{i}{\phi_{ \ al p h a}^ {\ si gma *} ( i) \phi_ {\g amma}^{ \sigma}(i) 
\p hi _{\beta}^{\s i gm a'*}(i)\ph i_{ \delta}^{\si gma '}(i)} \, .$$ Omitt ing  theoff-di a gonalmatrix el em e nts ($ ( \alpha, \ b et a)\n ot=
(\gamma,\delt a )$ )  in Eq.(7) yie lds th eH ar t r ee- Foc k approxim at ion.T he reas o nf o r  do i ng this willbe discusse d  su bseque nt ly. 
The b asicpa r ame ters of our  the oreticalframew o rk aret he elec tron n umb er$n$, th elon g- r ang e  C oul o mbinteract io n$U$,theH u b b ardint erac tionstrength $U_{ Hub }$,a ndthe d isord er a mp litud e $W$.  Equa ti ons (3-7) const itut e a set o f s el f-c on siste n t equa tio nsto yiel d the e i gen st a t e s$\phi_{\alpha}^{\s ig m a }, \epsilon _{\alp h a} ^{ \ sigma}$an d t he e n e rgy $ E_{G } $for a gi ven se t  o fspin oc cu pation  n umb ers  ($n_ { \upa rrow}, n_{\down arrow } $) with $n_{\u p arrow}+n_{\do w na r r ow } =n$.  At  every step  oft he n umer i ca l c a lcula tion,  t h em atrix elements \[de sc ribedby Eq s.(4-6)\] of$H^{eff}$w i l l be cal cula t ed for the spin o ccupa tion numbe r s that g ive t he minim um energy i n the pr evi ous  st ep. T he  first Slater d eter mi nant of  Eq .(3) wi llbefor med  b y noninte ractingei ge nf un cti ons o f  Hamilto ni an\[ Eq. (1)\] with $ U=U_{ Hub} =0 $. Var ying th e  s p i n oc cu pa tion  nu mb ers s o th a t $ M_{z}\i n [0,1]$wec alcu la te  $E_{G} (M_{z})$ whos eminimum va lu e g ives t h e  spin ma gnetization $M_{z}$ oft he syst em.  Thestab ility ofthe  self- con s istent  solut ionsis  ve r y  sens i t iv e a tlarger Cou l o mbinter ac tion  streng ths, forcing us to kee p a relativel y s mall v al ues of the  i n ter a c tion parameter$U$. Becau se of  the spinr ota ti onal sy mmetryof th e  Hamilt onian, on e limitsth e nu m e ric al calcula tions on ly to pos i tives pi n mag net izatio n($n _{\up arrow} \ geq  n_{\ downar ro w}$).
Incr ea sing the  amplitude of disorderthe sp ectru m o f noninte rac t ing  Hamilton ian$H_{0}=t\s um_ {<i ,j>,\ sig m a}c_{ i\si g ma }^{ \ dagge r}c_ { j\sigma}+ \s um_ { i ,\ sigma}\omeg a _ { i}n _{i,\ sig m a}$ ev olve s from weakly tos trongly locali zedr e gim e w h ileth e eigenvalue s tat is t i cs exhib it s a crossov er betwe en Wigne r andPoisso n distr i b ut i ons [@ zh
delta}^{\sigma'}(j)}
{|i-j|}\,,$$ $$\label{}
C_{\alpha_\sigma,\beta \sigma'}
^{\gamma_\sigma,\delta \sigma'}=\sum_{i}
{\phi_{\alpha}^{\sigma*}(i)\phi_{\gamma}^{\sigma}(i)
\phi_{\beta}^{\sigma'*}(i)\phi_{\delta}^{\sigma'}(i)}\,.$$ Omitting the_off-diagonal matrix_elements_($(\alpha,\beta)\not=
(\gamma,\delta)$) in_Eq.(7)_yields the Hartree-Fock_approximation. The reason_for doing this will_be discussed subsequently.

The_basic_parameters of our theoretical framework are the electron number $n$, the long-range Coulomb interaction_$U$,_the Hubbard_interaction_strength_$U_{Hub}$, and the disorder amplitude_$W$. Equations (3-7) constitute a_set of_self-consistent equations to yield the eigenstates $\phi_{\alpha}^{\sigma},\epsilon_{\alpha}^{\sigma}$ and_the_energy $E_{G}$ for_a given set of spin occupation numbers ($n_{\uparrow},n_{\downarrow}$) with_$n_{\uparrow}+n_{\downarrow}=n$. At every step of the_numerical calculation, the_matrix_elements_\[described by Eqs.(4-6)\] of_$H^{eff}$ will be calculated for the_spin occupation numbers that give the_minimum energy in the previous step. The_first Slater determinant of Eq.(3) will_be formed by noninteracting eigenfunctions_of Hamiltonian_\[Eq.(1)\] with $U=U_{Hub}=0$. Varying the_spin occupation numbers_so that_$M_{z}\in [0,1]$ we_calculate $E_{G}(M_{z})$ whose minimum value gives_the spin magnetization_$M_{z}$ of the system. The stability_of_the self-consistent solutions_is_very_sensitive at_larger Coulomb interaction_strengths,_forcing us_to_keep a relatively small values of_the_interaction parameter $U$. Because of the spin_rotational symmetry of the_Hamiltonian,_one limits the numerical_calculations only to positive spin_magnetization ($n_{\uparrow}\geq n_{\downarrow}$).

Increasing the amplitude of_disorder the_spectrum of_noninteracting Hamiltonian $H_{0}=t\sum_{<i,j>,\sigma}c_{i\sigma}^{\dagger}c_{j\sigma}
+\sum_{i,\sigma}\omega_{i}n_{i,\sigma}$ evolves from weakly to strongly localized regime while_the eigenvalue statistics exhibits a crossover_between Wigner and Poisson_distributions [@zh