
\left(
\begin{array}{cc}
(\partial_{\tau'}-\mu_{\sigma})\delta_{\sigma\sigma'}+t^2\sum\limits_{\langle0i\rangle,\langle0j\rangle}G^{1}_{\sigma\sigma',ij}(\tau,\tau') & t^{2}\sum\limits_{\langle0i\rangle,\langle0j\rangle}G^{2}_{\sigma\sigma',ij}(\tau,\tau')\\
t^{2}\sum\limits_{\langle0i\rangle,\langle0j\rangle}G^{2\ast}_{\sigma\sigma',ij}(\tau',\tau) & (-\partial_{\tau'}-\mu_{\sigma})\delta_{\sigma\sigma'}+t^2\sum\limits_{\langle0i\rangle,\langle0j\rangle}G^{1}_{\sigma\sigma',ij}(\tau',\tau)\\
\end{array}
\right)
\label{eq:weiss}$$ and the superfluid order parameters is given by $\phi_{i,\sigma}(\tau)\equiv\langle a_{i,\sigma(\tau)}\rangle_{0}$ with $\langle\cdots\rangle_{0}$ indicating that the expectation value is calculated in the cavity system[@Hubur_08; @Li2011Tunable]. We find that the effective action (\[eq:BDMFT\]) is represented by an Anderson impurity Hamiltonian $$\begin{aligned}
\label{eq:AndersonHamiltion}
\hat H_{A}=& &-t\sum_{\sigma}(\phi^{\ast}_{\sigma}\hat a_{0\sigma}+\mathrm{h.c.})+\frac{U_{0}}{2}\hat n_{0}(\hat n_{0}-1)+\frac{U_{2}}{2}(\bm{\hat{S}}^{2}_{0}-2\hat n_{0})+q\sum_{i,\sigma}\sigma^{2}\hat n_{i,\sigma}-\sum_{\sigma}\mu_{0\sigma}\hat n_{0\sigma}+\sum_{l}\epsilon_{l}\hat b^{\dagger}_{l}\hat b_{l} \\
& &+\sum_{l,\sigma}(V_{\sigma,l}\hat a_{0\sigma}\hat b^{\dagger}_{l}+W_{\sigma,l}\hat a_{0\sigma}\hat b_{l}+\mathrm{h.c.}). \nonumber\end{
\left (
 \begin{array}{cc } 
 (\partial_{\tau'}-\mu_{\sigma})\delta_{\sigma\sigma'}+t^2\sum\limits_{\langle0i\rangle,\langle0j\rangle}G^{1}_{\sigma\sigma',ij}(\tau,\tau') & t^{2}\sum\limits_{\langle0i\rangle,\langle0j\rangle}G^{2}_{\sigma\sigma',ij}(\tau,\tau')\\ 
 t^{2}\sum\limits_{\langle0i\rangle,\langle0j\rangle}G^{2\ast}_{\sigma\sigma',ij}(\tau',\tau) & (-\partial_{\tau'}-\mu_{\sigma})\delta_{\sigma\sigma'}+t^2\sum\limits_{\langle0i\rangle,\langle0j\rangle}G^{1}_{\sigma\sigma',ij}(\tau',\tau)\\ 
 \end{array } 
 \right) 
 \label{eq: weiss}$$ and the superfluid order parameters is given by $ \phi_{i,\sigma}(\tau)\equiv\langle a_{i,\sigma(\tau)}\rangle_{0}$ with $ \langle\cdots\rangle_{0}$ indicating that the anticipation value is account in the cavity system[@Hubur_08; @Li2011Tunable ]. We find that the effective action   (\[eq: BDMFT\ ]) is typify by an Anderson impurity Hamiltonian $ $ \begin{aligned } 
 \label{eq: AndersonHamiltion } 
 \hat H_{A}= & & -t\sum_{\sigma}(\phi^{\ast}_{\sigma}\hat a_{0\sigma}+\mathrm{h.c.})+\frac{U_{0}}{2}\hat n_{0}(\hat n_{0}-1)+\frac{U_{2}}{2}(\bm{\hat{S}}^{2}_{0}-2\hat n_{0})+q\sum_{i,\sigma}\sigma^{2}\hat n_{i,\sigma}-\sum_{\sigma}\mu_{0\sigma}\hat n_{0\sigma}+\sum_{l}\epsilon_{l}\hat b^{\dagger}_{l}\hat b_{l } \\ 
 & & + \sum_{l,\sigma}(V_{\sigma, l}\hat a_{0\sigma}\hat b^{\dagger}_{l}+W_{\sigma, l}\hat a_{0\sigma}\hat b_{l}+\mathrm{h.c. }). \nonumber\end {

\lefh(
\begin{array}{cc}
(\partial_{\tau'}-\mm_{\sigma})\delta_{\sigma\sigma'}+t^2\snm\limita_{\langle0i\fangle,\langle0j\rangle}G^{1}_{\sigma\sigla',uj}(\tau,\ucu') & t^{2}\sum\limits_{\langle0k\rangle,\lajgle0j\rantle}G^{2}_{\wugma\sigma',io}(\fau,\tau')\\
t^{2}\sum\limjbs_{\lanylx0i\rangle,\langle0j\tangle}G^{2\ast}_{\sicma\sigma',ij}(\tau',\tdu) & (-\'artial_{\tau'}-\mu_{\sigma})\delta_{\sigma\sigma'}+t^2\sum\jimits_{\lsnhle0i\rangle,\langje0j\rsggle}F^{1}_{\sigma\sigma',ij}(\tau',\tau)\\
\end{array}
\right)
\lagel{eq:wepss}$$ and the supergluid order parameters is hivej by $\phi_{i,\sigma}(\tau)\eeuiv\langle q_{i,\sidna(\tau)}\rangle_{0}$ dith $\langlt\cbots\rangle_{0}$ jndicating that the expectation valuz is calculqtwd lt the cavitb systvm[@Hubur_08; @Li2011Tunable]. We xind thst the effectiye acvion (\[wq:BDMFT\]) is representev by an Anderson impority Hamintknian $$\begin{alignee}
\lqbel{ex:AndarsovYamkltjoi}
\haf H_{A}=& &-t\dum_{\aigma}(\phi^{\asf}_{\sigma}\hat a_{0\wigma}+\mathrm{h.c.})+\frac{U_{0}}{2}\hst b_{0}(\hat n_{0}-1)+\frac{U_{2}}{2}(\bm{\gat{S}}^{2}_{0}-2\hae g_{0})+q\sum_{i,\sigma}\sigma^{2}\hat n_{i,\sigma}-\sum_{\sigma}\mu_{0\sicma}\gat n_{0\sigma}+\sum_{l}\epsilon_{l}\hqt b^{\dagger}_{l}\hat b_{l} \\
& &+\sul_{l,\sigma}(V_{\sygma,l}\hat a_{0\sigma}\hat b^{\dagger}_{l}+W_{\sigma,l}\hat a_{0\sigma}\hat b_{n}+\mathcm{f.c.}). \uinumbde\ejd{
\left( \begin{array}{cc} (\partial_{\tau'}-\mu_{\sigma})\delta_{\sigma\sigma'}+t^2\sum\limits_{\langle0i\rangle,\langle0j\rangle}G^{1}_{\sigma\sigma',ij}(\tau,\tau') & t^{2}\sum\limits_{\langle0i\rangle,\langle0j\rangle}G^{2}_{\sigma\sigma',ij}(\tau,\tau')\\ t^{2}\sum\limits_{\langle0i\rangle,\langle0j\rangle}G^{2\ast}_{\sigma\sigma',ij}(\tau',\tau) & \right) and the order parameters is $\langle\cdots\rangle_{0}$ that the expectation is calculated in cavity system[@Hubur_08; @Li2011Tunable]. We find that effective action (\[eq:BDMFT\]) is represented by an Anderson impurity Hamiltonian $$\begin{aligned} \label{eq:AndersonHamiltion} \hat &-t\sum_{\sigma}(\phi^{\ast}_{\sigma}\hat a_{0\sigma}+\mathrm{h.c.})+\frac{U_{0}}{2}\hat n_{0}(\hat n_{0}-1)+\frac{U_{2}}{2}(\bm{\hat{S}}^{2}_{0}-2\hat n_{0})+q\sum_{i,\sigma}\sigma^{2}\hat n_{i,\sigma}-\sum_{\sigma}\mu_{0\sigma}\hat n_{0\sigma}+\sum_{l}\epsilon_{l}\hat b^{\dagger}_{l}\hat b_{l} \\ & &+\sum_{l,\sigma}(V_{\sigma,l}\hat a_{0\sigma}\hat a_{0\sigma}\hat \nonumber\end{

\left(
\begin{array}{cc}
(\partial_{\taU'}-\mu_{\sigma})\deLta_{\siGma\SigMa'}+T^2\sum\LimiTs_{\langle0i\ranglE,\LangLe0j\rangle}G^{1}_{\sigma\sigma',ij}(\Tau,\taU') & t^{2}\SUm\liMItS_{\langLe0i\rangLE,\lANGle0J\rAnGle}g^{2}_{\sIGmA\sigmA',ij}(\Tau,\tau')\\
t^{2}\Sum\limits_{\lAngLe0I\rangle,\langlE0J\rAngle}G^{2\ast}_{\sIgmA\sigma',ij}(\tau',\tAu) & (-\pArtial_{\TaU'}-\mu_{\SIgma})\dEltA_{\sigmA\sigma'}+T^2\Sum\limIts_{\langle0I\rANgle,\laNGle0j\ranGLE}G^{1}_{\SigmA\sigma',ij}(\tau',\tau)\\
\end{ARrAY}
\right)
\label{eq:wEiss}$$ anD tHE sUPErfLuiD order paraMeTers iS Given by $\PHi_{I,\SIGma}(\TAu)\equiv\langle A_{i,\sigma(\tau)}\rANglE_{0}$ with $\lAnGle\CDots\raNgle_{0}$ iNdICatIng that the eXpecTation valUe is caLCulated IN the cavIty sysTem[@hubUr_08; @Li2011tUnAbLe]. WE fINd tHAt The EFfeCtive actIoN (\[eQ:BDMFt\]) is rEPRESentEd bY an ANdersOn impurity HamIltOniaN $$\BegIn{aliGned}
\lAbel{Eq:anderSonHamIltioN}
\hAt H_{A}=& &-t\sum_{\sigma}(\phI^{\ast}_{\Sigma}\hat a_{0\SigMa}+\MatHrM{h.c.})+\frAC{U_{0}}{2}\hat n_{0}(\Hat N_{0}-1)+\frAc{U_{2}}{2}(\bm{\haT{S}}^{2}_{0}-2\hat n_{0})+q\SUm_{i,\SiGMA}\SiGma^{2}\hat n_{i,\sigma}-\sum_{\siGmA}\MU_{0\sIgma}\hat n_{0\Sigma}+\sUM_{l}\EpSIlon_{l}\hat B^{\dAggEr}_{l}\hAT B_{l} \\
& &+\sum_{L,\sigMA}(V_{\Sigma,l}\haT a_{0\sigmA}\HaT b^{\Dagger}_{l}+w_{\sIgma,l}\hAt A_{0\siGma}\Hat b_{l}+\MAthrM{h.c.}). \nonUmber\end{

\left(
\begin{array}{cc}(\partial_ {\tau '}- \mu _{ \sig ma}) \delta_{\sigma \ sigm a'}+t^2\sum\limits_{\l angle 0i \ rang l e, \lang le0j\ra n gl e } G^{ 1} _{ \si gm a \s igma' ,ij }(\tau, \tau') & t ^{2 }\ sum\limits_{ \ la ngle0i\ran gle ,\langle0j\r ang le}G^{ 2} _{\ s igma\ sig ma',i j}(\ta u ,\tau' )\\
t^{2} \s u m\limi t s_{\lan g l e0 i\ra ngle,\langle0j\ra n gl e }G^{2\ast}_{\s igma\s ig m a' , i j}( \ta u',\tau) &  ( -\par t ial_{\t a u' } - \ mu_ { \sigma})\delt a_{\sigma\s i gma '}+t^2 \s um\ l imits_ {\lan gl e 0i\ rangle,\lan gle0 j\rangle} G^{1}_ { \sigma\ s igma',i j}(\ta u', \ta u)\\ \e nd {ar ra y }
\ r ig ht) \la bel{eq:w ei ss }$$ a nd t h e s uper flu id o rderparameters is  gi venb y $ \phi_ {i,\s igma }( \tau) \equiv \lang le  a_{i,\sigma(\t au)} \rangle_{ 0}$  w ith  $ \lang l e\cdot s\r ang le_{0}$  indica t ing  t h a t  t he expectation val ue i scalculat ed int he  c a vity sys te m[@ Hubu r _ 08; @ Li20 1 1T unable].  We fi n dth at theef fectiv eact ion  (\[e q :BDM FT\])is repre sente d  by an Anderso n  impurity Ham i lt o n ia n  $$\ beg in{aligned} 
\la b el{e q:An d er son H amilt ion}\h a tH _{A}=& &-t\sum_{\si gm a}(\ph i^{\a st}_{\sigma}\ hat a_{0\s i g m a}+\math rm{h . c. } )+\frac{U_{0}} {2}\h at n_{0}(\ h at n_{0} -1)+\ frac{U_{ 2}}{2}(\b m { \hat{S}} ^{2 }_{ 0}- 2\h a t  n _{0})+q\sum_{ i , \sig ma }\sigma ^{2 }\hat n _{i ,\s igm a}- \s um_{\sigm a}\mu_{0 \s ig ma }\ hat  n_{0 \ sigma}+\ su m_{ l} \ep silon _ {l}\ha t b^{ \dag ge r} _ {l} \hat b_ { l} \ \
&&+ \s um_{ l,\ si gma}( V_{\ s igm a,l}\ha t a_{0\si gma } \hat  b ^{ \dagger }_{l}+W_{\sig ma ,l}\hat a_ {0 \si gma}\h a t  b_{l}+\ mathrm{h.c.}). \nonumbe r \end{

\left(
\begin{array}{cc}
(\partial_{\tau'}-\mu_{\sigma})\delta_{\sigma\sigma'}+t^2\sum\limits_{\langle0i\rangle,\langle0j\rangle}G^{1}_{\sigma\sigma',ij}(\tau,\tau') &_t^{2}\sum\limits_{\langle0i\rangle,\langle0j\rangle}G^{2}_{\sigma\sigma',ij}(\tau,\tau')\\
t^{2}\sum\limits_{\langle0i\rangle,\langle0j\rangle}G^{2\ast}_{\sigma\sigma',ij}(\tau',\tau) &_(-\partial_{\tau'}-\mu_{\sigma})\delta_{\sigma\sigma'}+t^2\sum\limits_{\langle0i\rangle,\langle0j\rangle}G^{1}_{\sigma\sigma',ij}(\tau',\tau)\\
\end{array}
\right)
\label{eq:weiss}$$ and the superfluid_order parameters_is_given by_$\phi_{i,\sigma}(\tau)\equiv\langle_a_{i,\sigma(\tau)}\rangle_{0}$ with $\langle\cdots\rangle_{0}$_indicating that the_expectation value is calculated_in the cavity_system[@Hubur_08;_@Li2011Tunable]. We find that the effective action (\[eq:BDMFT\]) is represented by an Anderson impurity Hamiltonian_$$\begin{aligned}
\label{eq:AndersonHamiltion}
\hat_H_{A}=& &-t\sum_{\sigma}(\phi^{\ast}_{\sigma}\hat_a_{0\sigma}+\mathrm{h.c.})+\frac{U_{0}}{2}\hat_n_{0}(\hat_n_{0}-1)+\frac{U_{2}}{2}(\bm{\hat{S}}^{2}_{0}-2\hat n_{0})+q\sum_{i,\sigma}\sigma^{2}\hat n_{i,\sigma}-\sum_{\sigma}\mu_{0\sigma}\hat n_{0\sigma}+\sum_{l}\epsilon_{l}\hat b^{\dagger}_{l}\hat_b_{l} \\
& &+\sum_{l,\sigma}(V_{\sigma,l}\hat a_{0\sigma}\hat b^{\dagger}_{l}+W_{\sigma,l}\hat_a_{0\sigma}\hat b_{l}+\mathrm{h.c.})._\nonumber\end{