} - \frac{1} {4} n^f_j {\hat t}_{jl} ),
\label{heff2}\end{aligned}$$ where ${\bf S}_j=\frac{1}{2}\sum_{ss^\prime} {\tilde f}^\dagger_{js}  {\boldsymbol \sigma}_{ss^\prime}  {\tilde f}_{js^\prime}$ is the $f$-electron spin on site $j$, $
{\bf s}_{jl} = \frac{1}{2}\sum_{ss^\prime} (c^\dagger_{js} {\boldsymbol \sigma}_{ss^\prime} c_{l s^\prime}
+c^\dagger_{ls} {\boldsymbol \sigma}_{ss^\prime} c_{j s^\prime})
$ is the conduction electron “bond" spin (${\boldsymbol \sigma}$ are the Pauli Matrices), and $
{\hat t}_{jl} =  \sum_{s} (c^{\dagger}_{js} c^{\;}_{ls} + c^{\dagger}_{ls} c^{\;}_{js}).
$ The constrained operators ${\tilde f}^{\dagger}_{j \sigma}= f^{\dagger}_{j \sigma} (1-n^f_{j {\bar \sigma}})$ do not allow for double occupancy of the $f$-orbitals. ${\tilde f}^{\dagger}_{{\bf k} \sigma}$ creates an electron in the $f$-band with well defined momentum ${\bf k}$: ${\tilde f}^{\dagger}_{{\bf k} \sigma} = \frac{1}{\sqrt{N}} \sum_{j} e^{i {\bf k} \cdot {\bf r}_j} {\tilde f}^{\dagger}_{j \sigma}$.

The original broad band dispersion, $\epsilon_{\bf k}$, is renormalized to $
{\tilde \epsilon}_{\bf k}=\epsilon_{\bf k} + \frac{V^2_{\bf k}}{\epsilon_{\bf k}} = \epsilon_{\bf k} - 2 \frac{V^2_0}{t} \gamma_{\bf k}$. The $f$-electrons acquire an effective dispersion $
{\tilde\epsilon}_{f \bf k}=\epsilon_f - 2 {\tilde t}_f \gamma_{\bf k},
$ where $
{\tilde{t}_f} = - V_0^2/t
$ is the effective hopping between nearest-neighbor $f$-orbitals (here we neglected terms ${\
} - \frac{1 } { 4 } n^f_j { \hat t}_{jl }), 
 \label{heff2}\end{aligned}$$ where $ { \bf S}_j=\frac{1}{2}\sum_{ss^\prime } { \tilde f}^\dagger_{js }   { \boldsymbol \sigma}_{ss^\prime }   { \tilde f}_{js^\prime}$ is the $ f$-electron spin on site $ j$, $ 
 { \bf s}_{jl } = \frac{1}{2}\sum_{ss^\prime } (c^\dagger_{js } { \boldsymbol \sigma}_{ss^\prime } c_{l s^\prime } 
 + c^\dagger_{ls } { \boldsymbol \sigma}_{ss^\prime } c_{j s^\prime }) 
 $ is the conduction electron “ bond " spin ($ { \boldsymbol \sigma}$ are the Pauli Matrices), and $ 
 { \hat t}_{jl } =   \sum_{s } (c^{\dagger}_{js } c^{\;}_{ls } + c^{\dagger}_{ls } c^{\;}_{js }). 
 $ The constrained operator $ { \tilde f}^{\dagger}_{j \sigma}= f^{\dagger}_{j \sigma } (1 - n^f_{j { \bar \sigma}})$ do not leave for double occupancy of the $ f$-orbitals. $ { \tilde f}^{\dagger}_{{\bf k } \sigma}$ creates an electron in the $ f$-band with well define momentum $ { \bf k}$: $ { \tilde f}^{\dagger}_{{\bf k } \sigma } = \frac{1}{\sqrt{N } } \sum_{j } e^{i { \bf k } \cdot { \bf r}_j } { \tilde f}^{\dagger}_{j \sigma}$. 

 The original broad band distribution, $ \epsilon_{\bf k}$, is renormalized to $ 
 { \tilde \epsilon}_{\bf k}=\epsilon_{\bf k } + \frac{V^2_{\bf k}}{\epsilon_{\bf k } } = \epsilon_{\bf k } - 2 \frac{V^2_0}{t } \gamma_{\bf k}$. The $ f$-electrons acquire an effective dispersion $ 
 { \tilde\epsilon}_{f \bf k}=\epsilon_f - 2 { \tilde t}_f \gamma_{\bf k }, 
 $ where $ 
 { \tilde{t}_f } = - V_0 ^ 2 / thymine 
 $ is the effective hopping between nearest - neighbor $ f$-orbitals (here we neglected term $ { \
} - \fgac{1} {4} n^f_j {\hat t}_{jl} ),
\label{henf2}\end{aligned}$$ whete ${\bf S}_j=\hrac{1}{2}\sum_{as^\prime} {\gilde f}^\dagger_{js}  {\boldsymbol \ditma}_{ss^\kgime}  {\tilde f}_{js^\prime}$ ks the $f$-epectron wpin in site $j$, $
{\ur s}_{jl} = \nxac{1}{2}\suj_{ds^\prnmx} (c^\dagger_{js} {\bolcsymbol \sicma}_{ss^\prime} c_{l v^\pfile}
+c^\dagger_{ls} {\boldsymbol \sigma}_{ss^\prime} c_{j s^\promf})
$ is the condustiom elednrin “bond" spin (${\boldsymbol \sigma}$ are tht Pauli Matrices), amd $
{\hat t}_{jl} =  \sum_{s} (c^{\dagger}_{jd} c^{\;}_{ld} + c^{\dagger}_{ls} c^{\;}_{js}).
$ Tje constraibed jperators ${\tilae f}^{\dagger}_{j \sigma}= f^{\dagfer}_{j \sigma} (1-n^f_{j {\bar \sigma}})$ do not alloc for doublg kcfopancy of thx $f$-orbptals. ${\tilde f}^{\dagger}_{{\bf k} \sigms}$ creates an ekecvron in the $f$-band with wenl defined momentui ${\bf k}$: ${\tindz f}^{\dagger}_{{\bf k} \sigma} = \drqc{1}{\sqrj{N}} \suk_{j} e^{k {\bf k} \ddpt {\bf r}_j} {\timde f}^{\dagged}_{j \sigma}$.

The original broad banc qpxpersion, $\epsjlon_{\bf k}$, is renormalized to $
{\tilde \epsilon}_{\bf k}=\epvilkn_{\bf k} + \frac{V^2_{\bf k}}{\epsilob_{\bf k}} = \epsilon_{\bf k} - 2 \vrac{V^2_0}{t} \gaima_{\bf k}$. The $f$-electrons acquire an effective dispession $
{\gilbc\cpsiuin}_{v \bf k}=\epsilon_f - 2 {\tilde t}_f \gamma_{\bf k},
$ where $
{\tilqs{t}_g} = - V_0^2/t
$ is the efnective hopping beywfem nearest-neighcor $f$-oxgifals (here we neglefted tetms ${\
} - \frac{1} {4} n^f_j {\hat t}_{jl} where S}_j=\frac{1}{2}\sum_{ss^\prime} {\tilde {\boldsymbol \sigma}_{ss^\prime} {\tilde on $j$, $ {\bf = \frac{1}{2}\sum_{ss^\prime} (c^\dagger_{js} \sigma}_{ss^\prime} c_{l s^\prime} +c^\dagger_{ls} {\boldsymbol \sigma}_{ss^\prime} s^\prime}) $ is the conduction electron “bond" spin (${\boldsymbol \sigma}$ are the Pauli and $ {\hat t}_{jl} = \sum_{s} (c^{\dagger}_{js} c^{\;}_{ls} + c^{\dagger}_{ls} c^{\;}_{js}). $ The operators f}^{\dagger}_{j f^{\dagger}_{j (1-n^f_{j {\bar \sigma}})$ do not allow for double occupancy of the $f$-orbitals. ${\tilde f}^{\dagger}_{{\bf k} \sigma}$ an electron in the $f$-band with well defined ${\bf k}$: ${\tilde f}^{\dagger}_{{\bf \sigma} = \frac{1}{\sqrt{N}} \sum_{j} e^{i k} {\bf r}_j} f}^{\dagger}_{j The broad band dispersion, k}$, is renormalized to $ {\tilde \epsilon}_{\bf k}=\epsilon_{\bf k} + \frac{V^2_{\bf k}}{\epsilon_{\bf k}} = \epsilon_{\bf k} - \frac{V^2_0}{t} \gamma_{\bf $f$-electrons acquire effective $ \bf k}=\epsilon_f - t}_f \gamma_{\bf k}, $ where $ V_0^2/t $ is the effective hopping between nearest-neighbor (here we terms ${\
} - \frac{1} {4} n^f_j {\hat t}_{jl} ),
\label{heff2}\end{Aligned}$$ wheRe ${\bf S}_J=\frAc{1}{2}\sUm_{Ss^\prIme} {\tIlde f}^\dagger_{js}  {\bOLdsyMbol \sigma}_{ss^\prime}  {\tilde f}_{Js^\priMe}$ IS the $F$-ElEctroN spin on SItE $J$, $
{\Bf s}_{Jl} = \FrAc{1}{2}\sUm_{SS^\pRime} (c^\DagGer_{js} {\boLdsymbol \siGma}_{Ss^\Prime} c_{l s^\primE}
+C^\dAgger_{ls} {\bolDsyMbol \sigma}_{ss^\pRimE} c_{j s^\prImE})
$ is THe conDucTion eLectroN “Bond" spIn (${\boldsymBoL \Sigma}$ aRE the PauLI maTricEs), and $
{\hat t}_{jl} =  \sum_{s} (c^{\dAGgER}_{js} c^{\;}_{ls} + c^{\dagger}_{lS} c^{\;}_{js}).
$ ThE cONsTRAinEd oPerators ${\tiLdE f}^{\dagGEr}_{j \sigmA}= F^{\dAGGEr}_{j \SIgma} (1-n^f_{j {\bar \sigMa}})$ do not alloW For Double OcCupANcy of tHe $f$-orBiTAls. ${\Tilde f}^{\daggeR}_{{\bf k} \Sigma}$ creaTes an eLEctron iN The $f$-banD with wEll DefIned MOmEnTum ${\Bf K}$: ${\TilDE f}^{\DagGEr}_{{\bF k} \sigma} = \fRaC{1}{\sQrt{N}} \sUm_{j} e^{I {\BF K} \Cdot {\Bf r}_J} {\tilDe f}^{\daGger}_{j \sigma}$.

The OriGinaL BroAd banD dispErsiOn, $\EpsilOn_{\bf k}$, iS renoRmAlized to $
{\tilde \epSiloN}_{\bf k}=\epsilOn_{\bF k} + \FraC{V^2_{\Bf k}}{\epSIlon_{\bf K}} = \epSilOn_{\bf k} - 2 \frAc{V^2_0}{t} \gamMA_{\bf K}$. THE $F$-ElEctrons acquire an efFeCTIvE dispersIon $
{\tilDE\ePsILon}_{f \bf k}=\ePsIloN_f - 2 {\tiLDE t}_f \gaMma_{\bF K},
$ wHere $
{\tildE{t}_f} = - V_0^2/t
$ iS ThE eFfectivE hOpping BeTweEn nEaresT-NeigHbor $f$-oRbitals (hEre we NEglected terms ${\
} - \frac{1} {4} n^f_j {\h at t}_{jl}  ),
\ lab el{ he ff2} \end {aligned}$$ wh e re $ {\bf S}_j=\frac{1}{2}\ sum_{ ss ^ \pri m e}  {\ti lde f}^ \ da g g er_ {j s}   { \b o ld symbo l \ sigma}_ {ss^\prime } {\ tilde f}_{js ^ \p rime}$ isthe  $f$-electro n s pin on  s ite $j$,$
{ \bf s }_{jl} = \fra c{1}{2}\s um _ {ss^\p r ime} (c ^ \ da gger _{js} {\boldsymbo l  \ s igma}_{ss^\pri me} c_ {l s^ \ p rim e}+c^\dagger _{ ls} { \ boldsym b ol \ s igm a }_{ss^\prime}  c_{j s^\pr i me} )
$ is  t hec onduct ion e le c tro n “bond" sp in ( ${\boldsy mbol \ s igma}$a re thePauliMat ric es), an d$
{ \h a t t } _{ jl} = \sum_{s}  ( c^ {\dag ger} _ { j s } c^ {\; }_{l s} +c^{\dagger}_{ ls}  c^{ \ ;}_ {js}) .
$ T he c on strai ned op erato rs  ${\tilde f}^{\ dagg er}_{j \s igm a} = f ^{ \dagg e r}_{j\si gma } (1-n^ f_{j {\ b ar\s i g m a} })$ do not allow f or d ou ble occu pancyo fth e  $f$-orb it als . ${ \ t ildef}^{ \ da gger}_{{ \bf k} \s ig ma}$ cr ea tes an  e lec tro n int he $ f$-ban d with w ell d e fined momentum ${\bf k}$: ${ \ ti l d ef }^{\ dag ger}_{{\bfk} \ s igma } =\ fr ac{ 1 }{\sq rt{N} }\ su m _{j} e^{i {\bf k} \ cd ot {\b f r}_ j} {\tilde f} ^{\dagger} _ { j  \sigma} $.
T he original broad  band  dispersio n , $\epsi lon_{ \bf k}$,  is renor m a lized to  ${\t ild e \ e p si lon}_{\bf k}= \ e psil on _{\bf k } +  \frac{ V^2 _{\ bfk}} {\ epsilon_{ \bf k}}=\e ps il on_ {\bfk } - 2 \f ra c{V ^2 _0} {t} \ g amma_{ \bf k }$.Th e$ f$- electro n sa c quir ean  eff ect iv e dis pers i on$
{\til de\epsilo n}_ { f \b fk} =\epsil on_f - 2 {\ti ld e t}_f \ga mm a_{ \bf k} , $ where$
{\tilde{t}_f} = - V_0 ^ 2/t
$ i s t he ef fect ive hoppi ngbetwee n n e arest- neighb or $f $- orb i t als ( h e re  we  n eglected t e r ms${\
} -_\frac{1} {4}_n^f_j {\hat t}_{jl} ),
\label{heff2}\end{aligned}$$_where ${\bf_S}_j=\frac{1}{2}\sum_{ss^\prime}_{\tilde f}^\dagger_{js}__{\boldsymbol \sigma}_{ss^\prime} _{\tilde f}_{js^\prime}$ is_the $f$-electron spin on_site $j$, $
{\bf_s}_{jl}_= \frac{1}{2}\sum_{ss^\prime} (c^\dagger_{js} {\boldsymbol \sigma}_{ss^\prime} c_{l s^\prime}
+c^\dagger_{ls} {\boldsymbol \sigma}_{ss^\prime} c_{j s^\prime})
$ is the conduction_electron_“bond" spin_(${\boldsymbol_\sigma}$_are the Pauli Matrices), and_$
{\hat t}_{jl} =  \sum_{s}_(c^{\dagger}_{js} c^{\;}_{ls}_+ c^{\dagger}_{ls} c^{\;}_{js}).
$ The constrained operators ${\tilde f}^{\dagger}_{j_\sigma}=_f^{\dagger}_{j \sigma} (1-n^f_{j_{\bar \sigma}})$ do not allow for double occupancy of_the $f$-orbitals. ${\tilde f}^{\dagger}_{{\bf k} \sigma}$_creates an electron_in_the_$f$-band with well defined_momentum ${\bf k}$: ${\tilde f}^{\dagger}_{{\bf k}_\sigma} = \frac{1}{\sqrt{N}} \sum_{j} e^{i {\bf_k} \cdot {\bf r}_j} {\tilde f}^{\dagger}_{j \sigma}$.

The_original broad band dispersion, $\epsilon_{\bf k}$,_is renormalized to $
{\tilde \epsilon}_{\bf_k}=\epsilon_{\bf k}_+ \frac{V^2_{\bf k}}{\epsilon_{\bf k}} =_\epsilon_{\bf k} -_2 \frac{V^2_0}{t}_\gamma_{\bf k}$. The_$f$-electrons acquire an effective dispersion $
{\tilde\epsilon}_{f_\bf k}=\epsilon_f -_2 {\tilde t}_f \gamma_{\bf k},
$ where_$
{\tilde{t}_f}_= - V_0^2/t
$_is_the_effective hopping_between nearest-neighbor $f$-orbitals_(here_we neglected_terms_${\