bf q} = V/N$ where $V$ is the strength of the impurity potential. In what follows we use standard finite temperature Green’s function formalism [@DS88; @BF04]. Because of the existence of two sub-lattices, the Green’s function can be written as a $2\times 2$ matrix: $$\begin{aligned}
\bm G_{\sigma}(\bm k,{\bf p},\tau) = \left(\begin{array}{cc}
G_{AA,\sigma} (\bm k,{\bf p}, \tau) \hspace{0.5cm} & G_{AB,\sigma}(\bm k,{\bf p},\tau) \\
G_{BA,\sigma} (\bm k,{\bf p},\tau) \hspace{0.5cm} & G_{BB,\sigma}(\bm k,{\bf p},\tau)  
\end{array}\right) \,,\end{aligned}$$ where $$\begin{aligned}
G_{AA,\sigma}(\bm k,{\bf p}, \tau) &=& - \langle {\cal T} a_{\bm k,\sigma}(\tau) a_{\bm p,\sigma}^{\dag}(0)
\rangle \,,
\nonumber
\\
G_{AB,\sigma}(\bm k,{\bf p}, \tau) &=& - \langle {\cal T} a_{\bm k,\sigma}(\tau) b_{\bm p,\sigma}^{\dag}(0)
\rangle \,,
\nonumber
\\
G_{BA,\sigma}(\bm k,{\bf p},\tau) &=& - \langle {\cal T} b_{\bm k,\sigma}(\tau) a_{\bm p,\sigma}^{\dag}(0)
\rangle \,,
\nonumber
\\
G_{BB,\sigma}(\bm k,{\bf p},\tau) &=& - \langle {\cal T} b_{\bm k,\sigma}(\tau) b_{\bm p,\sigma}^{\dag}(0)
\rangle \,,
\label{defg}\end{aligned}$$ where $\tau$ is the “imaginary” time, and ${\cal T}$ is the time ordering operator.

For a single impurity the Green’s function can be written as $\bm G(\bm k,{\bf p},\tau) = \delta_{\bm k,\bm p} \bm G(\bm k,\tau)$, where [@DS88]: $$\bm G(\bm k,\omega_n)=\bm G^0(\bm
bf q } = V / N$ where $ V$ is the strength of the impurity potential. In what follows we practice standard finite temperature Green ’s affair formalism [ @DS88; @BF04 ]. Because of the existence of two sub - lattice, the Green ’s routine can be written as a $ 2\times 2 $ matrix: $ $ \begin{aligned } 
 \bm G_{\sigma}(\bm k,{\bf p},\tau) = \left(\begin{array}{cc } 
 G_{AA,\sigma } (\bm k,{\bf p }, \tau) \hspace{0.5 cm } & G_{AB,\sigma}(\bm k,{\bf p},\tau) \\ 
 G_{BA,\sigma } (\bm k,{\bf p},\tau) \hspace{0.5 curium } & G_{BB,\sigma}(\bm k,{\bf p},\tau)  
 \end{array}\right) \,,\end{aligned}$$ where $ $ \begin{aligned } 
 G_{AA,\sigma}(\bm k,{\bf p }, \tau) & = & - \langle { \cal T } a_{\bm k,\sigma}(\tau) a_{\bm p,\sigma}^{\dag}(0) 
 \rangle \, , 
 \nonumber 
 \\ 
 G_{AB,\sigma}(\bm k,{\bf p }, \tau) & = & - \langle { \cal T } a_{\bm k,\sigma}(\tau) b_{\bm p,\sigma}^{\dag}(0) 
 \rangle \, , 
 \nonumber 
 \\ 
 G_{BA,\sigma}(\bm k,{\bf p},\tau) & = & - \langle { \cal T } b_{\bm k,\sigma}(\tau) a_{\bm p,\sigma}^{\dag}(0) 
 \rangle \, , 
 \nonumber 
 \\ 
 G_{BB,\sigma}(\bm k,{\bf p},\tau) & = & - \langle { \cal T } b_{\bm k,\sigma}(\tau) b_{\bm p,\sigma}^{\dag}(0) 
 \rangle \, , 
 \label{defg}\end{aligned}$$ where $ \tau$ is the “ fanciful ” time, and $ { \cal T}$ is the time ordering operator. 

 For a single impurity the Green ’s routine can be written as $ \bm G(\bm k,{\bf p},\tau) = \delta_{\bm k,\bm p } \bm G(\bm k,\tau)$, where [ @DS88 ]: $ $ \bm G(\bm k,\omega_n)=\bm G^0(\bm
bf e} = V/N$ where $V$ is the strtngth of the imputiry potxntial. Jn what wollows we use standard finive twmperqture Green’s function wormalism [@DS88; @BF04]. Bwcauww of the eejstence of twk sub-nettices, the Gregn’s function can be writtet xs a $2\times 2$ matrix: $$\begin{aligned}
\bm G_{\sidma}(\bm k,{\nf p},\tau) = \left(\begyn{argar}{cc}
G_{ZA,\sigma} (\bm k,{\bf p}, \tau) \hspace{0.5cm} & G_{AB,\aigma}(\bm k,{\bf p},\tau) \\
G_{BA,\xigma} (\bm k,{\bf p},\tau) \hspace{0.5cm} & G_{BH,\sigma}(\bm k,{\bf p},\tau)  
\end{array}\ritht) \,,\qbd{aligned}$$ whdre $$\begin{aligned}
G_{AA,\sigja}(\bm k,{\bf p}, \tau) &=& - \langle {\cal T} a_{\cm k,\sngma}(\tau) a_{\bm p,\wigld}^{\dag}(0)
\rangle \,,
\ionumbvr
\\
G_{AB,\sigma}(\bm k,{\bf p}, \tag) &=& - \lanble {\cal T} a_{\bm l,\sijma}(\tqu) b_{\bm p,\sigma}^{\dag}(0)
\ranglx \,,
\nonumber
\\
G_{BA,\sigma}(\bm k,{\bf p},\tau) &=& - \langle {\cal T} b_{\bm k,\wigma}(\jau) a_{\tm p,\rugmx}^{\daf}(0)
\rengme \,,
\nonkmbxr
\\
G_{BB,\sigma}(\bj k,{\bf p},\tau) &=& - \langle {\cal T} b_{\bm l,\sytma}(\tau) b_{\bm p,\sjgma}^{\dad}(0)
\rwngle \,,
\label{defg}\end{aligned}$$ where $\tau$ is uhe “ijaginary” time, and ${\cal T}$ is the time ordering operator.

Sor a single impurity the Green’s function can be frittxn as $\bi T(\bl k,{\bf p},\tau) = \delta_{\bm k,\bm p} \bm G(\bm k,\tau)$, where [@DF88]: $$\bk N(\bm k,\omega_n)=\bm G^0(\bi
bf q} = V/N$ where $V$ is of impurity potential. what follows we function [@DS88; @BF04]. Because the existence of sub-lattices, the Green’s function can be as a $2\times 2$ matrix: $$\begin{aligned} \bm G_{\sigma}(\bm k,{\bf p},\tau) = \left(\begin{array}{cc} G_{AA,\sigma} k,{\bf p}, \tau) \hspace{0.5cm} & G_{AB,\sigma}(\bm k,{\bf p},\tau) \\ G_{BA,\sigma} (\bm k,{\bf p},\tau) & k,{\bf \end{array}\right) where $$\begin{aligned} G_{AA,\sigma}(\bm k,{\bf p}, \tau) &=& - \langle {\cal T} a_{\bm k,\sigma}(\tau) a_{\bm p,\sigma}^{\dag}(0) \rangle \nonumber \\ G_{AB,\sigma}(\bm k,{\bf p}, \tau) &=& - {\cal T} a_{\bm k,\sigma}(\tau) p,\sigma}^{\dag}(0) \rangle \,, \nonumber \\ k,{\bf &=& - {\cal b_{\bm a_{\bm p,\sigma}^{\dag}(0) \rangle \nonumber \\ G_{BB,\sigma}(\bm k,{\bf p},\tau) &=& - \langle {\cal T} b_{\bm k,\sigma}(\tau) b_{\bm p,\sigma}^{\dag}(0) \rangle \,, \label{defg}\end{aligned}$$ $\tau$ is time, and T}$ the ordering operator. For impurity the Green’s function can be G(\bm k,{\bf p},\tau) = \delta_{\bm k,\bm p} \bm k,\tau)$, where $$\bm G(\bm k,\omega_n)=\bm G^0(\bm
bf q} = V/N$ where $V$ is the strength oF the impuriTy potEntIal. in What FollOws we use standaRD finIte temperature Green’s fuNctioN fORmalISm [@dS88; @BF04]. BEcause oF ThE EXisTeNcE of TwO SuB-lattIceS, the GreEn’s functioN caN bE written as a $2\tIMeS 2$ matrix: $$\begIn{aLigned}
\bm G_{\sigMa}(\bM k,{\bf p},\tAu) = \LefT(\Begin{ArrAy}{cc}
G_{aA,\sigmA} (\Bm k,{\bf p}, \Tau) \hspace{0.5Cm} & g_{aB,\sigmA}(\Bm k,{\bf p},\tAU) \\
g_{Ba,\sigMa} (\bm k,{\bf p},\tau) \hspace{0.5CM} & G_{bb,\sigma}(\bm k,{\bf p},\taU)  
\end{arRaY}\RiGHT) \,,\enD{alIgned}$$ where $$\BeGin{alIGned}
G_{AA,\SIgMA}(\BM k,{\bF P}, \tau) &=& - \langle {\cal t} a_{\bm k,\sigma}(\tAU) a_{\bM p,\sigmA}^{\dAg}(0)
\rANgle \,,
\noNumbeR
\\
G_{ab,\siGma}(\bm k,{\bf p}, \taU) &=& - \lanGle {\cal T} a_{\bM k,\sigmA}(\Tau) b_{\bm p,\SIgma}^{\dag}(0)
\Rangle \,,
\NonUmbEr
\\
G_{Ba,\SiGmA}(\bm K,{\bF P},\taU) &=& - \LaNglE {\Cal t} b_{\bm k,\sigMa}(\TaU) a_{\bm p,\SigmA}^{\DAG}(0)
\RangLe \,,
\nOnumBer
\\
G_{Bb,\sigma}(\bm k,{\bf p},\tAu) &=& - \lAnglE {\Cal t} b_{\bm k,\Sigma}(\Tau) b_{\Bm P,\sigmA}^{\dag}(0)
\raNgle \,,
\lAbEl{defg}\end{aligneD}$$ wheRe $\tau$ is thE “imAgInaRy” Time, aND ${\cal T}$ iS thE tiMe orderIng operATor.

foR A SInGle impurity the GreeN’s FUNcTion can bE writtEN aS $\bM g(\bm k,{\bf p},\tAu) = \DelTa_{\bm K,\BM p} \bm G(\Bm k,\tAU)$, wHere [@DS88]: $$\bm g(\bm k,\omEGa_N)=\bM G^0(\bm
bf q} = V/N$ where $V$ isthe streng th of  th e i mp urit y po tential. In wh a t fo llows we use standardfinit et empe r at ure G reen’sf un c t ion  f or mal is m  [ @DS88 ; @ BF04].Because of  th eexistence of tw o sub-latt ice s, the Green ’sfuncti on  ca n  be w rit ten a s a $2 \ times2$ matrix :$ $\begi n {aligne d } 
\ bm G _{\sigma}(\bm k,{ \ bf p},\tau) = \le ft(\be gi n {a r r ay} {cc }
G_{AA,\s ig ma} ( \ bm k,{\ b fp } ,  \t a u) \hspace{0. 5cm} & G_{A B ,\s igma}( \b m k , {\bf p },\ta u) \\G_{BA,\sigm a} ( \bm k,{\b f p},\ t au) \hs p ace{0.5 cm} &G_{ BB, \sig m a} (\ bmk, { \bf p} ,\t a u) 
\end{a rr ay }\rig ht)\ , , \ end{ ali gned }$$ w here $$\begin {al igne d }
G _{AA, \sigm a}(\ bm  k,{\ bf p},  \tau )&=& - \langle { \cal  T} a_{\b m k ,\ sig ma }(\ta u ) a_{\ bmp,\ sigma}^ {\dag}( 0 )
\ ra n g l e\,,
\nonumber
\\
G _{ A B ,\ sigma}(\ bm k,{ \ bf  p } , \tau)&= & -  \la n g le {\ calT }a_{\bm k ,\sigm a }( \t au) b_{ \b m p,\s ig ma} ^{\ dag}( 0 )
\r angle\,,
\non umber \\
G_{BA,\sigm a }(\bm k,{\bfp }, \ t au )  &=&  -\langle {\c al T }  b_{ \bmk ,\ sig m a}(\t au) a _{ \ bm p,\sigma}^{\dag}(0) 
\ rangle  \,,\nonumber
\\G_{BB,\sig m a } (\bm k,{ \bfp }, \ tau) &=& - \la ngle{\cal T} b _ {\bm k,\ sigma }(\tau)b_{\bm p, \ s igma}^{\ dag }(0 )
\ ran g l e\,,
\label{de f g }\en d{ aligned }$$  where$\t au$  is  th e“imaginar y” time,  a nd  $ {\ cal  T}$i s the ti me  or de rin g ope r ator.
Fora si ng le imp urity t h eG r een’ sfu ncti onca n bewrit t enas $\bm  G(\bm k, {\b f  p}, \t au ) = \de lta_{\bm k,\b mp} \bm G(\ bm  k, \tau)$ , where [@ DS88]: $$\bm G(\bm k,\o m ega_n)= \bm  G^0( \bm
bf q}_= V/N$_where $V$ is the_strength of_the_impurity potential._In_what follows we_use standard finite_temperature Green’s function formalism_[@DS88; @BF04]. Because_of_the existence of two sub-lattices, the Green’s function can be written as a $2\times_2$_matrix: $$\begin{aligned}
\bm_G_{\sigma}(\bm_k,{\bf_p},\tau) = \left(\begin{array}{cc}
G_{AA,\sigma} (\bm k,{\bf_p}, \tau) \hspace{0.5cm} & G_{AB,\sigma}(\bm_k,{\bf p},\tau)_\\
G_{BA,\sigma} (\bm k,{\bf p},\tau) \hspace{0.5cm} & G_{BB,\sigma}(\bm k,{\bf_p},\tau)_ 
\end{array}\right) \,,\end{aligned}$$_where $$\begin{aligned}
G_{AA,\sigma}(\bm k,{\bf p}, \tau) &=& - \langle {\cal_T} a_{\bm k,\sigma}(\tau) a_{\bm p,\sigma}^{\dag}(0)
\rangle \,,
\nonumber
\\
G_{AB,\sigma}(\bm_k,{\bf p}, \tau)_&=&_-_\langle {\cal T} a_{\bm_k,\sigma}(\tau) b_{\bm p,\sigma}^{\dag}(0)
\rangle \,,
\nonumber
\\
G_{BA,\sigma}(\bm k,{\bf p},\tau)_&=& - \langle {\cal T} b_{\bm_k,\sigma}(\tau) a_{\bm p,\sigma}^{\dag}(0)
\rangle \,,
\nonumber
\\
G_{BB,\sigma}(\bm k,{\bf p},\tau) &=&_- \langle {\cal T} b_{\bm k,\sigma}(\tau)_b_{\bm p,\sigma}^{\dag}(0)
\rangle \,,
\label{defg}\end{aligned}$$ where $\tau$_is the_“imaginary” time, and ${\cal T}$_is the time_ordering operator.

For_a single impurity_the Green’s function can be written_as $\bm G(\bm_k,{\bf p},\tau) = \delta_{\bm k,\bm p}_\bm_G(\bm k,\tau)$, where_[@DS88]:_$$\bm_G(\bm k,\omega_n)=\bm_G^0(\bm