langle
t^{\text{Ga}}_{\sigma \sigma} \rangle_{T}+ c\langle
t^{\text{Mn}}_{\sigma \sigma} \rangle_{T}=0\end{aligned}$$

$\langle... \rangle_{\text{dis},T}$ denotes configuration average and thermal average at temperature $T$ for the spin operator, $c$ is the concentration of Mn impurities.

After lengthy calculations we get, $$\begin{aligned}
\Sigma_{\sigma}= \frac{\Large \langle \sum_{\lambda}
x_{\lambda}\hat{\bf c}^{\lambda}_{\sigma}(\hat{\bf
d}^{\lambda}_{\sigma})^{-1} \Large \rangle_{T}} {\Large \langle
\sum_{\lambda} x_{\lambda}(\hat{\bf d}^{\lambda}_{\sigma})^{-1}
\Large \rangle_{T}} \label{self}\end{aligned}$$ The sum runs over all constituents, in our binary system $x_{\text{Mn}}=c$ (resp. $x_{\text{Ga}}=1-c$). $$\begin{aligned}
{\hat{\bf c}}^{\text{Mn}}_{\sigma}&=&
z_{\sigma}\frac{J}{2}\hat{\bf S}^{z}+\large( \frac{J^{2}}{4}
[S(S+1)-(\hat{\bf S}^{z})^2 -z_{\sigma}\hat{\bf S}^{z}]  \nonumber \\
&&+ [z_{\sigma}\frac{J}{2}\hat{\bf
S}^{z}+\Sigma_{\sigma}][\frac{J}{2}(z_{-\sigma}\hat{\bf S}^{z}-1
)-\Sigma_{-\sigma}]   \Large) G_{-\sigma}\end{aligned}$$ and, $$\begin{aligned}
{\hat{\bf d}}^{\text{Mn}}_{\sigma}&=&1-G_{-\sigma}[\frac{J}{2}(z_{-\sigma}\hat{\bf S}^{z}-1)-\Sigma_{-\sigma}] \nonumber  \\
&&- G_{\sigma}[{\hat{\bf
c}}^{\text{Mn}}_{\sigma}-z_{\sigma}\frac{J}{2}\hat{\bf S}^{z}]\end{aligned}$$ with $z_{\sigma}=1$ (resp.$-1$) for spin $\uparrow$ (resp.$\downarrow$). $\bar{G}_{\sigma} = \frac{1}{N} \sum_{\bf q} \bar{G}_{\
langle 
 t^{\text{Ga}}_{\sigma \sigma } \rangle_{T}+ c\langle 
 t^{\text{Mn}}_{\sigma \sigma } \rangle_{T}=0\end{aligned}$$ 

 $ \langle... \rangle_{\text{dis},T}$ denotes configuration average and thermal average at temperature $ T$ for the spin operator, $ c$ is the concentration of Mn impurity. 

 After drawn-out calculations we get, $ $ \begin{aligned } 
 \Sigma_{\sigma}= \frac{\Large \langle \sum_{\lambda } 
 x_{\lambda}\hat{\bf c}^{\lambda}_{\sigma}(\hat{\bf 
 d}^{\lambda}_{\sigma})^{-1 } \Large \rangle_{T } } { \Large \langle 
 \sum_{\lambda } x_{\lambda}(\hat{\bf d}^{\lambda}_{\sigma})^{-1 } 
 \Large \rangle_{T } } \label{self}\end{aligned}$$ The sum run over all constituent, in our binary system $ x_{\text{Mn}}=c$ (resp. $ x_{\text{Ga}}=1 - c$). $ $ \begin{aligned } 
 { \hat{\bf c}}^{\text{Mn}}_{\sigma}&= & 
 z_{\sigma}\frac{J}{2}\hat{\bf S}^{z}+\large (\frac{J^{2}}{4 } 
 [ S(S+1)-(\hat{\bf S}^{z})^2 -z_{\sigma}\hat{\bf S}^{z } ]   \nonumber \\ 
 & & + [ z_{\sigma}\frac{J}{2}\hat{\bf 
 S}^{z}+\Sigma_{\sigma}][\frac{J}{2}(z_{-\sigma}\hat{\bf S}^{z}-1 
) -\Sigma_{-\sigma } ]    \Large) G_{-\sigma}\end{aligned}$$ and, $ $ \begin{aligned } 
 { \hat{\bf d}}^{\text{Mn}}_{\sigma}&=&1 - G_{-\sigma}[\frac{J}{2}(z_{-\sigma}\hat{\bf S}^{z}-1)-\Sigma_{-\sigma } ] \nonumber   \\ 
 & & - G_{\sigma}[{\hat{\bf 
 c}}^{\text{Mn}}_{\sigma}-z_{\sigma}\frac{J}{2}\hat{\bf S}^{z}]\end{aligned}$$ with $ z_{\sigma}=1 $ (resp.$-1 $) for spin $ \uparrow$ (resp.$\downarrow$). $ \bar{G}_{\sigma } = \frac{1}{N } \sum_{\bf q } \bar{G}_{\
lanhle
t^{\text{Ga}}_{\sigma \sigma} \rakgle_{T}+ c\langle
t^{\terr{Mn}}_{\sigka \sigja} \rangld_{T}=0\end{aligned}$$

$\langle... \rangle_{\texv{dis},R}$ denites configuration avefage and nhermal aceraje at temperaturx $T$ for bke spjk opexavor, $c$ is the cokcentration of Mn impurither.

Ayter lengthy calculations we get, $$\begyn{alignrd}
\Digma_{\sigma}= \fras{\Larbq \lahgle \sum_{\lambda}
x_{\lambda}\hat{\bf c}^{\lambda}_{\aigma}(\hau{\bf
d}^{\lambda}_{\sigma})^{-1} \Latge \rangle_{T}} {\Large \langle
\sul_{\lamhda} x_{\lambda}(\hat{\bf d}^{\pambda}_{\sigma})^{-1}
\Oargq \rangle_{T}} \labdl{self}\end{aligned}$$ The som runs over all constituents, in our yinary systgm $x_{\hgxt{Mn}}=c$ (resp. $e_{\text{Gw}}=1-c$). $$\begin{aligkvd}
{\hat{\bf c}}^{\text{Mm}}_{\sigma}&=&
z_{\sigma}\frsc{J}{2}\iat{\bd S}^{z}+\large( \frac{J^{2}}{4}
[S(S+1)-(\hat{\bh S}^{z})^2 -z_{\sigma}\hat{\bf S}^{z}]  \nonumber \\
&&+ [v_{\sigma}\frac{J}{2}\hat{\bf
S}^{z}+\Signa_{\wigma}][\xrac{B}{2}(z_{-\siena}\hxt{\br X}^{z}-1
)-\Aigma_{-\slgme}]   \Large) G_{-\aigma}\end{alitned}$$ and, $$\begin{aligntd}
{\hwn{\nf d}}^{\text{Mn}}_{\sifma}&=&1-G_{-\sidmw}[\frac{J}{2}(z_{-\sigma}\hat{\bf S}^{z}-1)-\Sigma_{-\sigma}] \nonumber  \\
&&- F_{\sigma}[{\hat{\bf
c}}^{\text{Mn}}_{\sigma}-z_{\sigma}\frac{J}{2}\hat{\bf S}^{z}]\enf{aligned}$$ rith $z_{\sigma}=1$ (resp.$-1$) for spin $\uparrow$ (resp.$\downarrow$). $\bdr{G}_{\sijmx} = \yvwz{1}{B} \dum_{\bf q} \bar{G}_{\
langle t^{\text{Ga}}_{\sigma \sigma} \rangle_{T}+ c\langle t^{\text{Mn}}_{\sigma \sigma} \rangle_{\text{dis},T}$ configuration average thermal average at operator, is the concentration Mn impurities. After calculations we get, $$\begin{aligned} \Sigma_{\sigma}= \frac{\Large \sum_{\lambda} x_{\lambda}\hat{\bf c}^{\lambda}_{\sigma}(\hat{\bf d}^{\lambda}_{\sigma})^{-1} \Large \rangle_{T}} {\Large \langle \sum_{\lambda} x_{\lambda}(\hat{\bf d}^{\lambda}_{\sigma})^{-1} \Large \rangle_{T}} The sum runs over all constituents, in our binary system $x_{\text{Mn}}=c$ (resp. $x_{\text{Ga}}=1-c$). {\hat{\bf z_{\sigma}\frac{J}{2}\hat{\bf \frac{J^{2}}{4} S}^{z})^2 -z_{\sigma}\hat{\bf S}^{z}] \nonumber \\ &&+ [z_{\sigma}\frac{J}{2}\hat{\bf S}^{z}+\Sigma_{\sigma}][\frac{J}{2}(z_{-\sigma}\hat{\bf S}^{z}-1 )-\Sigma_{-\sigma}] \Large) G_{-\sigma}\end{aligned}$$ and, $$\begin{aligned} {\hat{\bf d}}^{\text{Mn}}_{\sigma}&=&1-G_{-\sigma}[\frac{J}{2}(z_{-\sigma}\hat{\bf \nonumber \\ &&- G_{\sigma}[{\hat{\bf c}}^{\text{Mn}}_{\sigma}-z_{\sigma}\frac{J}{2}\hat{\bf S}^{z}]\end{aligned}$$ with $z_{\sigma}=1$ for spin $\uparrow$ (resp.$\downarrow$). = \frac{1}{N} \sum_{\bf q} \bar{G}_{\
langle
t^{\text{Ga}}_{\sigma \sigma} \raNgle_{T}+ c\langLe
t^{\teXt{MN}}_{\siGmA \sigMa} \raNgle_{T}=0\end{aligneD}$$

$\LangLe... \rangle_{\text{dis},T}$ denoteS confIgURatiON aVeragE and theRMaL AVerAgE aT teMpERaTure $T$ For The spin Operator, $c$ iS thE cOncentration OF MN impuritieS.

AfTer lengthy caLcuLationS wE geT, $$\Begin{AliGned}
\SIgma_{\siGMa}= \frac{\large \langLe \SUm_{\lambDA}
x_{\lambdA}\HAt{\Bf c}^{\lAmbda}_{\sigma}(\hat{\bf
d}^{\lAMbDA}_{\sigma})^{-1} \Large \ranGle_{T}} {\LaRgE \LaNGLe
\sUm_{\lAmbda} x_{\lambDa}(\Hat{\bf D}^{\Lambda}_{\sIGmA})^{-1}
\lARge \RAngle_{T}} \label{seLf}\end{aligneD}$$ the Sum runS oVer ALl consTitueNtS, In oUr binary sysTem $x_{\Text{Mn}}=c$ (reSp. $x_{\texT{ga}}=1-c$). $$\begiN{Aligned}
{\Hat{\bf c}}^{\TexT{Mn}}_{\SigmA}&=&
Z_{\sIgMa}\fRaC{j}{2}\haT{\Bf s}^{z}+\lARge( \Frac{J^{2}}{4}
[S(S+1)-(\hAt{\Bf s}^{z})^2 -z_{\siGma}\hAT{\BF s}^{z}]  \noNumBer \\
&&+ [z_{\Sigma}\Frac{J}{2}\hat{\bf
S}^{z}+\SIgmA_{\sigMA}][\frAc{J}{2}(z_{-\sIgma}\hAt{\bf s}^{z}-1
)-\sigma_{-\Sigma}]   \LArge) G_{-\SiGma}\end{aligned}$$ anD, $$\begIn{aligned}
{\Hat{\Bf D}}^{\teXt{mn}}_{\sigMA}&=&1-G_{-\sigmA}[\frAc{J}{2}(Z_{-\sigma}\hAt{\bf S}^{z}-1)-\SIGma_{-\SiGMA}] \NoNumber  \\
&&- G_{\sigma}[{\hat{\bf
c}}^{\TeXT{mn}}_{\Sigma}-z_{\siGma}\fraC{j}{2}\hAt{\BF S}^{z}]\end{alIgNed}$$ With $Z_{\SIgma}=1$ (rEsp.$-1$) fOR sPin $\uparrOw$ (resp.$\DOwNaRrow$). $\bar{g}_{\sIgma} = \frAc{1}{n} \suM_{\bf Q} \bar{G}_{\
langle
t^{\text{Ga}}_{\sig ma \sigma}  \ran gle _{T }+  c\l angl e
t^{\text{Mn} } _{\s igma \sigma} \rangle_{ T}=0\ en d {ali g ne d}$$
$\lang l e. . .  \r an gl e_{ \t e xt {dis} ,T} $ denot es configu rat io n average an d  t hermal ave rag e at tempera tur e $T$fo r t h e spi n o perat or, $c $  is th e concent ra t ion of Mn impu r i ti es.
After lengthy ca l cu l ations we get,  $$\be gi n {a l i gne d}\Sigma_{\s ig ma}=\ frac{\L a rg e \ lan g le \sum_{\lam bda}
x_{\la m bda }\hat{ \b f c } ^{\lam bda}_ {\ s igm a}(\hat{\bf 
d}^ {\lambda} _{\sig m a})^{-1 }  \Large  \rang le_ {T} } {\ L ar ge  \l an g le\ su m_{ \ lam bda} x_{ \l am bda}( \hat { \ b f  d}^ {\l ambd a}_{\ sigma})^{-1}\La rge\ ran gle_{ T}} \ labe l{ self} \end{a ligne d} $$ The sum runs  ove r all con sti tu ent s,  in o u r bina rysys tem $x_ {\text{ M n}} =c $ ( re sp. $x_{\text{Ga}} =1 - c $) . $$\beg in{ali g ne d} {\hat{\b fc}} ^{\t e x t{Mn} }_{\ s ig ma}&=&
z _{\sig m a} \f rac{J}{ 2} \hat{\ bf  S} ^{z }+\la r ge(\frac{ J^{2}}{4 }
[S( S +1)-(\hat{\bfS }^{z})^2 -z_{ \ si g m a} \ hat{ \bf  S}^{z}]  \ nonu m ber\\
& & +[z_ { \sigm a}\fr ac { J} { 2}\hat{\bf
S}^{z}+\ Si gma_{\ sigma }][\frac{J}{2 }(z_{-\sig m a } \hat{\bf  S}^ { z} - 1
)-\Sigma_{-\ sigma }]   \Larg e ) G_{-\s igma} \end{ali gned}$$ a n d , $$\beg in{ ali gne d}{ \ ha t{\bf d}}^{\t e x t{Mn }} _{\sigm a}& =&1-G_{ -\s igm a}[ \fr ac {J}{2}(z_ {-\sigma }\ ha t{ \b f S }^{z} - 1)-\Sigm a_ {-\ si gma }] \n o number   \\&&-G_ {\ s igm a}[{\ha t {\ b f 
c}} ^{ \t ext{ Mn} }_ {\sig ma}- z _{\ sigma}\ frac{J}{2 }\h a t{\b fS} ^{z}]\e nd{aligned}$$  w ith $z_{\s ig ma} =1$ (r e s p.$-1$)for spin $\uparrow$ (re s p.$\dow nar row$) . $\ bar{G}_{\ sig ma} =\fr a c{1}{N } \sum _{\bf  q } \ b a r{G}_ { \ 
langle
t^{\text{Ga}}_{\sigma \sigma}_\rangle_{T}+ c\langle
t^{\text{Mn}}_{\sigma_\sigma} \rangle_{T}=0\end{aligned}$$

$\langle... \rangle_{\text{dis},T}$ denotes_configuration average_and_thermal average_at_temperature $T$ for_the spin operator,_$c$ is the concentration_of Mn impurities.

After_lengthy_calculations we get, $$\begin{aligned}
\Sigma_{\sigma}= \frac{\Large \langle \sum_{\lambda}
x_{\lambda}\hat{\bf c}^{\lambda}_{\sigma}(\hat{\bf
d}^{\lambda}_{\sigma})^{-1} \Large \rangle_{T}} {\Large \langle
\sum_{\lambda} x_{\lambda}(\hat{\bf d}^{\lambda}_{\sigma})^{-1}
\Large_\rangle_{T}}_\label{self}\end{aligned}$$ The_sum_runs_over all constituents, in our_binary system $x_{\text{Mn}}=c$ (resp. $x_{\text{Ga}}=1-c$)._$$\begin{aligned}
{\hat{\bf c}}^{\text{Mn}}_{\sigma}&=&
z_{\sigma}\frac{J}{2}\hat{\bf_S}^{z}+\large( \frac{J^{2}}{4}
[S(S+1)-(\hat{\bf S}^{z})^2 -z_{\sigma}\hat{\bf S}^{z}]  \nonumber \\
&&+_[z_{\sigma}\frac{J}{2}\hat{\bf
S}^{z}+\Sigma_{\sigma}][\frac{J}{2}(z_{-\sigma}\hat{\bf_S}^{z}-1
)-\Sigma_{-\sigma}]  _\Large) G_{-\sigma}\end{aligned}$$ and, $$\begin{aligned}
{\hat{\bf d}}^{\text{Mn}}_{\sigma}&=&1-G_{-\sigma}[\frac{J}{2}(z_{-\sigma}\hat{\bf S}^{z}-1)-\Sigma_{-\sigma}] \nonumber  \\
&&-_G_{\sigma}[{\hat{\bf
c}}^{\text{Mn}}_{\sigma}-z_{\sigma}\frac{J}{2}\hat{\bf S}^{z}]\end{aligned}$$ with $z_{\sigma}=1$ (resp.$-1$) for_spin $\uparrow$ (resp.$\downarrow$)._$\bar{G}_{\sigma}_=_\frac{1}{N} \sum_{\bf q} \bar{G}_{\