psi_{\downarrow} (-k') \psi_{\uparrow} (k'),$$

where the bare coupling constant $g$ is negative, to express the attractive (s-wave) interaction between the spin $\uparrow$ and $\downarrow$ fermionic atoms.

After the mean field (MF) approximation and the subsequent diagonalization of the expression for $H_{MF}$, we arrive at the following expression for thermodynamic potential:

$$\begin{aligned}
\Omega &=&-\frac{\Delta^2}{g}+\int_{k<k_1,~k>k_2} \frac{d^2 k}{(2\pi)^2} (\epsilon_k^b - {\cal E} _k^{\beta})
\nonumber
\\
&+& \int_{k_1}^{k_2} \frac{d^2 k}{(2\pi)^2} \epsilon_k^b
\nonumber
\\
&=&-\frac{\Delta^2}{g}+\int_{k<k_1,~k>k_2} \frac{d^2 k}{(2\pi)^2} (\epsilon_k^+ - E_k)
\nonumber
\\
&+& \int_{k_1}^{k_2} \frac{d^2 k}{(2\pi)^2} \epsilon_k^b,
\label{Omega1}\end{aligned}$$

where, for simplicity of notation we have labeled $\downarrow=a$, $\uparrow=b$. Here we have defined ${\cal{E}}_k^{a,b}= E_k \pm \epsilon_k^{-}$ are the quasiparticle excitations, with $E_k=\sqrt{ {\epsilon_k^{+}}^2+\Delta^2 }$, $\epsilon_k^{\pm} = \frac {\epsilon_k^a \pm \epsilon_k^b}{2}$ and the constant pairing gap is given by $\Delta = - g \int \frac{d^2 k}{(2\pi)^2} \langle \psi_{\downarrow}^{\dagger}(-k) \psi_{\uparrow}^{\dagger}(k) \rangle = \Delta^*$.

In Eq. (\[Omega1\])

$$\label{roots}
k_{1,2}^2 = \frac{k_a^2+k_b^2}{2} \pm \frac{1}{2}\sqrt{(k_b^2
psi_{\downarrow } (-k') \psi_{\uparrow } (k'),$$ 

 where the bare coupling constant $ g$ is negative, to carry the attractive (south - wave) interaction between the spin $ \uparrow$ and $ \downarrow$ fermionic atoms. 

 After the mean discipline (MF) approximation and the subsequent diagonalization of the expression for $ H_{MF}$, we arrive at the surveil expression for thermodynamic potential: 

 $ $ \begin{aligned } 
 \Omega & = & -\frac{\Delta^2}{g}+\int_{k < k_1,~k > k_2 } \frac{d^2 k}{(2\pi)^2 } (\epsilon_k^b - { \cal E } _ k^{\beta }) 
 \nonumber 
 \\ 
 & + & \int_{k_1}^{k_2 } \frac{d^2 k}{(2\pi)^2 } \epsilon_k^b 
 \nonumber 
 \\ 
 & = & -\frac{\Delta^2}{g}+\int_{k < k_1,~k > k_2 } \frac{d^2 k}{(2\pi)^2 } (\epsilon_k^+ - E_k) 
 \nonumber 
 \\ 
 & + & \int_{k_1}^{k_2 } \frac{d^2 k}{(2\pi)^2 } \epsilon_k^b, 
 \label{Omega1}\end{aligned}$$ 

 where, for ease of notation we have labeled $ \downarrow = a$, $ \uparrow = b$. Here we have define $ { \cal{E}}_k^{a, b}= E_k \pm \epsilon_k^{-}$ are the quasiparticle excitations, with $ E_k=\sqrt { { \epsilon_k^{+}}^2+\Delta^2 } $, $ \epsilon_k^{\pm } = \frac { \epsilon_k^a \pm \epsilon_k^b}{2}$ and the constant pair gap is given by $ \Delta = - g \int \frac{d^2 k}{(2\pi)^2 } \langle \psi_{\downarrow}^{\dagger}(-k) \psi_{\uparrow}^{\dagger}(k) \rangle = \Delta^*$. 

 In Eq.   (\[Omega1\ ]) 

 $ $ \label{roots } 
 k_{1,2}^2 = \frac{k_a^2+k_b^2}{2 } \pm \frac{1}{2}\sqrt{(k_b^2
psi_{\fownarrow} (-k') \psi_{\uparrow} (k'),$$

dhere the bare eiuplinj constznt $g$ is negative, to express the attcactuve (s-qave) interaction betwedn the sppn $\uparroq$ anv $\downarrow$ fermmknic atoms.

Aftsv the nean field (MF) spproximathon and the sutsdqbent diagonalization of the expressijn for $N_{MV}$, we arrive at the sollkwing expression for thermodynamic potentpal:

$$\begin{aligned}
\Omrga &=&-\frac{\Delta^2}{g}+\int_{k<k_1,~k>k_2} \frac{f^2 k}{(2\pl)^2} (\epsilon_k^b - {\cal E} _k^{\beta})
\nonumver
\\
&+& \ybt_{k_1}^{k_2} \frac{d^2 k}{(2\oi)^2} \epsilon_k^b
\nonumber
\\
&=&-\frzc{\Delta^2}{g}+\int_{k<k_1,~k>k_2} \frac{d^2 k}{(2\pi)^2} (\epsilun_k^+ - Z_k)
\nonumber
\\
&+& \unr_{k_1}^{k_2} \xrac{d^2 k}{(2\pi)^2} \e'silon_h^b,
\label{Omega1}\cmd{aligted}$$

wherr, for simpliciby of noration we have labelev $\downarrow=a$, $\uparrow=f$. Here we hcve defined ${\cal{E}}_k^{a,b}= E_j \pm \epvilot_k^{-}$ afw tfe suesilarticpe xxcitations, with $E_k=\sqrr{ {\epsilon_k^{+}}^2+\Delta^2 }$, $\epxijin_k^{\pm} = \frac {\elsilon_h^a \pm \epsilon_k^b}{2}$ and the constant pairing cap is given by $\Delta = - g \unt \frac{d^2 k}{(2\pi)^2} \langle \ksi_{\downarrjw}^{\dagger}(-k) \psi_{\uparrow}^{\dagger}(k) \rangle = \Delta^*$.

In Eq. (\[Omeca1\])

$$\labxl{foous}
j_{1,2}^2 = \frxx{k_w^2+k_b^2}{2} \pm \frac{1}{2}\sqrt{(k_b^2
psi_{\downarrow} (-k') \psi_{\uparrow} (k'),$$ where the bare $g$ negative, to the attractive (s-wave) and fermionic atoms. After mean field (MF) and the subsequent diagonalization of the for $H_{MF}$, we arrive at the following expression for thermodynamic potential: $$\begin{aligned} \Omega \frac{d^2 k}{(2\pi)^2} (\epsilon_k^b - {\cal E} _k^{\beta}) \nonumber \\ &+& \int_{k_1}^{k_2} \frac{d^2 k}{(2\pi)^2} \nonumber &=&-\frac{\Delta^2}{g}+\int_{k<k_1,~k>k_2} k}{(2\pi)^2} - E_k) \nonumber \\ &+& \int_{k_1}^{k_2} \frac{d^2 k}{(2\pi)^2} \epsilon_k^b, \label{Omega1}\end{aligned}$$ where, for simplicity of notation we labeled $\downarrow=a$, $\uparrow=b$. Here we have defined ${\cal{E}}_k^{a,b}= \pm \epsilon_k^{-}$ are the excitations, with $E_k=\sqrt{ {\epsilon_k^{+}}^2+\Delta^2 }$, = {\epsilon_k^a \pm and constant gap is given $\Delta = - g \int \frac{d^2 k}{(2\pi)^2} \langle \psi_{\downarrow}^{\dagger}(-k) \psi_{\uparrow}^{\dagger}(k) \rangle = \Delta^*$. In Eq. (\[Omega1\]) $$\label{roots} = \frac{k_a^2+k_b^2}{2}
psi_{\downarrow} (-k') \psi_{\uparrow} (k'),$$

wHere the barE coupLinG coNsTant $G$ is nEgative, to expreSS the Attractive (s-wave) interacTion bEtWEen tHE sPin $\upArrow$ anD $\DoWNArrOw$ FeRmiOnIC aToms.

AFteR the meaN field (MF) apProXiMation and the SUbSequent diaGonAlization of tHe eXpressIoN foR $h_{MF}$, we ArrIve at The folLOwing eXpression FoR ThermoDYnamic pOTEnTial:

$$\Begin{aligned}
\Omega &=&-\FRaC{\delta^2}{g}+\int_{k<k_1,~k>k_2} \fRac{d^2 k}{(2\pI)^2} (\ePSiLON_k^b - {\Cal e} _k^{\beta})
\nonuMbEr
\\
&+& \int_{K_1}^{K_2} \frac{d^2 k}{(2\PI)^2} \ePSILon_K^B
\nonumber
\\
&=&-\frac{\delta^2}{g}+\int_{k<k_1,~K>K_2} \frAc{d^2 k}{(2\pi)^2} (\EpSilON_k^+ - E_k)
\noNumbeR
\\
&+& \iNT_{k_1}^{k_2} \Frac{d^2 k}{(2\pi)^2} \epsIlon_K^b,
\label{OmEga1}\end{ALigned}$$

wHEre, for sImplicIty Of nOtatIOn We HavE lABelED $\dOwnARroW=a$, $\uparroW=b$. heRe we hAve dEFINEd ${\caL{E}}_k^{A,b}= E_k \Pm \epsIlon_k^{-}$ are the quAsiPartICle ExcitAtionS, witH $E_K=\sqrt{ {\EpsiloN_k^{+}}^2+\DelTa^2 }$, $\Epsilon_k^{\pm} = \frac {\ePsilOn_k^a \pm \epsIloN_k^B}{2}$ anD tHe conSTant paIriNg gAp is givEn by $\DelTA = - g \iNt \FRAC{d^2 K}{(2\pi)^2} \langle \psi_{\downarRoW}^{\DAgGer}(-k) \psi_{\uParrow}^{\DAgGeR}(K) \rangle = \DElTa^*$.

IN Eq. (\[OMEGa1\])

$$\labEl{roOTs}
K_{1,2}^2 = \frac{k_a^2+k_B^2}{2} \pm \fraC{1}{2}\SqRt{(K_b^2
psi_{\downarrow} (-k') \ps i_{\uparro w} (k '), $$
w here  the  bare coupling cons tant $g$ is negative,to ex pr e ss t h eattra ctive ( s -w a v e)in te rac ti o nbetwe enthe spi n $\uparro w$an d $\downarro w $fermionicato ms.

After t hemean f ie ld( MF) a ppr oxima tion a n d thesubsequen td iagona l ization o ftheexpression for $H _ {M F }$, we arriveat the  f o ll o w ing  ex pression f or  ther m odynami c  p o t e nti a l:

$$\begin{ aligned}
\O m ega  &=&-\ fr ac{ \ Delta^ 2}{g} +\ i nt_ {k<k_1,~k>k _2}\frac{d^2  k}{(2 \ pi)^2}( \epsilo n_k^b- { \ca l E} _k ^{ \be ta } )
\ n on umb e r
\ \
&+& \i nt _{ k_1}^ {k_2 } \ f rac{ d^2  k}{ (2\pi )^2} \epsilon _k^ b
\n o num ber
\ \
&=& -\fr ac {\Del ta^2}{ g}+\i nt _{k<k_1,~k>k_2}  \fr ac{d^2 k} {(2 \p i)^ 2}  (\ep s ilon_k ^+- E _k)
\no number\ \
& +& \ i nt _{k_1}^{k_2} \frac {d ^ 2  k }{(2\pi) ^2} \e p si lo n _k^b,
\l ab el{ Omeg a 1 }\end {ali g ne d}$$

wh ere, f o rsi mplicit yof not at ion  we  have labe led $\ downarro w=a$, $\uparrow=b$.H ere we have d e fi n e d$ {\ca l{E }}_k^{a,b}=  E_k \pm\eps i lo n_k ^ {-}$are t he qu a siparticle excitati on s, wit h $E_ k=\sqrt{ {\ep silon_k^{+ } } ^ 2+\Delta ^2 } $ ,$ \epsilon_k^{\p m} =\frac {\ep s ilon_k^a  \pm\epsilon _k^b}{2}$ a nd the c ons tan t p air i n ggap is givenb y  $\D el ta = -g \ int \fr ac{ d^2  k} {(2 \p i)^2} \la ngle \ps i_ {\ do wn arr ow}^{ \ dagger}( -k ) \ ps i_{ \upar r ow}^{\ dagge r}(k )\r a ngl e = \De l ta ^ * $.
In  E q. ( \[O me ga1\] )

$ $ \la bel{roo ts}
k_{1, 2}^ 2  = \ fr ac {k_a^2+ k_b^2}{2} \pm  \ frac{1}{2} \s qrt {(k_b^ 2 
psi_{\downarrow} (-k')_\psi_{\uparrow} (k'),$$

where_the bare coupling constant_$g$ is_negative,_to express_the_attractive (s-wave) interaction_between the spin_$\uparrow$ and $\downarrow$ fermionic_atoms.

After the mean_field_(MF) approximation and the subsequent diagonalization of the expression for $H_{MF}$, we arrive at_the_following expression_for_thermodynamic_potential:

$$\begin{aligned}
\Omega &=&-\frac{\Delta^2}{g}+\int_{k<k_1,~k>k_2} \frac{d^2 k}{(2\pi)^2} (\epsilon_k^b_- {\cal E} _k^{\beta})
\nonumber
\\
&+& \int_{k_1}^{k_2}_\frac{d^2 k}{(2\pi)^2}_\epsilon_k^b
\nonumber
\\
&=&-\frac{\Delta^2}{g}+\int_{k<k_1,~k>k_2} \frac{d^2 k}{(2\pi)^2} (\epsilon_k^+ - E_k)
\nonumber
\\
&+& \int_{k_1}^{k_2} \frac{d^2_k}{(2\pi)^2}_\epsilon_k^b,
\label{Omega1}\end{aligned}$$

where, for simplicity_of notation we have labeled $\downarrow=a$, $\uparrow=b$. Here we_have defined ${\cal{E}}_k^{a,b}= E_k \pm \epsilon_k^{-}$_are the quasiparticle_excitations,_with_$E_k=\sqrt{ {\epsilon_k^{+}}^2+\Delta^2 }$, $\epsilon_k^{\pm}_= \frac {\epsilon_k^a \pm \epsilon_k^b}{2}$ and_the constant pairing gap is given_by $\Delta = - g \int \frac{d^2_k}{(2\pi)^2} \langle \psi_{\downarrow}^{\dagger}(-k) \psi_{\uparrow}^{\dagger}(k) \rangle =_\Delta^*$.

In Eq. (\[Omega1\])

$$\label{roots}
k_{1,2}^2 = \frac{k_a^2+k_b^2}{2} \pm_\frac{1}{2}\sqrt{(k_b^2