  \begin{pmatrix}
    -\mu-v_Fk & 0 & 0 & -\Delta^{h}_{\vec{k}} \\
    0 & -\mu+v_Fk & \Delta^{h}_{\vec{k}}  & 0 \\
    0 &\Delta^{h*}_{\vec{k}} & \mu- v_Fk &0 \\
    -\Delta^{h*}_{\vec{k}} &0& 0&\mu+v_Fk
  \end{pmatrix},\end{aligned}$$ provided that the holon pairing parameter $\Delta^{h}_{\vec{k}}$ is p-wave like, i.e., $\Delta^h_{-\vec{k}} = -\Delta^h_{\vec{k}}$. The spectrum of quasiparticles consists of two decoupled branches, $$\begin{aligned}
\label{eq:67}
\epsilon_{\pm,\vec{k}}^h = \sqrt{(v_Fk\pm \mu)^2+|\Delta^h_{\vec{k}}|^2}.\end{aligned}$$ The free energy at temperature $T$ then reads $$\begin{aligned}
\label{eq:68}
F=-T\sum_{i=\pm, \vec{k}} \left[
  \ln(1+e^{-\frac{\epsilon_{i \vec{k}}}{T}})+\ln(1+
  e^{\frac{\epsilon_{i \vec{k}}}{T}}) \right].\end{aligned}$$ According to Hellman-Feynman theorem, we have the gap equation for order parameter $\Delta_{\vec{q}}^h$, $$\begin{aligned}
\label{eq:69}
&&\langle \hat{b}_{-\vec{q}}\hat{a}_{\vec{q}} \rangle = \frac{\partial F}{\partial
  \Delta_{\vec{q}}^{h*}} =-\sum_{i=\pm}
\frac{\Delta_{\vec{q}}^{h}}{2 \epsilon_{i,\vec{q}}^h}\tanh \left(
  \frac{\epsilon_{i,\vec{q}}}{2T} \right) \nonumber\\
&&\Delta_{\vec{k}}^h=\sum_{\vec{q}}V_{\text{eff}}(\vec{k}-\vec{q})
\langle \hat{b}_{-\vec{q}}\hat{a}_{\vec{q}} \rangle.\
\begin{pmatrix } 
     -\mu - v_Fk & 0 & 0 & -\Delta^{h}_{\vec{k } } \\ 
     0 & -\mu+v_Fk & \Delta^{h}_{\vec{k } }   & 0 \\ 
     0 & \Delta^{h*}_{\vec{k } } & \mu- v_Fk & 0 \\ 
     -\Delta^{h*}_{\vec{k } } & 0 & 0&\mu+v_Fk 
   \end{pmatrix},\end{aligned}$$ provided that the holon pairing parameter $ \Delta^{h}_{\vec{k}}$ is p - wave like, i.e., $ \Delta^h_{-\vec{k } } = -\Delta^h_{\vec{k}}$. The spectrum of quasiparticle consist of two decoupled branches, $ $ \begin{aligned } 
 \label{eq:67 } 
 \epsilon_{\pm,\vec{k}}^h = \sqrt{(v_Fk\pm \mu)^2+|\Delta^h_{\vec{k}}|^2}.\end{aligned}$$ The spare department of energy at temperature $ T$ then reads $ $ \begin{aligned } 
 \label{eq:68 } 
 F=-T\sum_{i=\pm, \vec{k } } \left [ 
   \ln(1+e^{-\frac{\epsilon_{i \vec{k}}}{T}})+\ln(1 + 
   e^{\frac{\epsilon_{i \vec{k}}}{T } }) \right].\end{aligned}$$ According to Hellman - Feynman theorem, we have the opening equation for order argument $ \Delta_{\vec{q}}^h$, $ $ \begin{aligned } 
 \label{eq:69 } 
 & & \langle \hat{b}_{-\vec{q}}\hat{a}_{\vec{q } } \rangle = \frac{\partial F}{\partial 
   \Delta_{\vec{q}}^{h * } } = -\sum_{i=\pm } 
 \frac{\Delta_{\vec{q}}^{h}}{2 \epsilon_{i,\vec{q}}^h}\tanh \left (
   \frac{\epsilon_{i,\vec{q}}}{2 metric ton } \right) \nonumber\\ 
 & & \Delta_{\vec{k}}^h=\sum_{\vec{q}}V_{\text{eff}}(\vec{k}-\vec{q }) 
 \langle \hat{b}_{-\vec{q}}\hat{a}_{\vec{q } } \rangle.\
  \bfgin{pmatrix}
    -\mu-v_Fk & 0 & 0 & -\Delta^{h}_{\vec{k}} \\
    0 & -\mu+v_Hk & \Delfa^{h}_{\vec{k}}  & 0 \\
    0 &\Delta^{h*}_{\vec{k}} & \mu- v_Fk &0 \\
    -\Deluc^{h*}_{\vec{k}} &0& 0&\mu+v_Fk
  \end{pmxtrix},\end{apigned}$$ peovived that the holon pairiky parzletex $\Velta^{h}_{\vec{k}}$ is p-eave like, h.e., $\Delta^h_{-\vec{k}} = -\Ddlca^h_{\vec{k}}$. The spectrum of quasiparticlqs consoshs of two decoopled frandhes, $$\begin{aligned}
\label{eq:67}
\epsilon_{\pm,\vsc{k}}^h = \sxrt{(v_Fk\pm \mu)^2+|\Dekta^h_{\vec{k}}|^2}.\end{aligned}$$ The fref enfrgy at temperaturf $T$ then reqds $$\fwgin{aligned}
\lxbel{eq:68}
F=-T\sum_{i=\pm, \vec{k}} \lert[
  \ln(1+e^{-\frac{\epsilon_{i \vec{k}}}{T}})+\ln(1+
  e^{\ffac{\epxilon_{i \vec{j}}}{T}}) \rihvt].\end{alignev}$$ Accogding to Hellman-Feynmdn theotem, we have thc gap eqyation for order paraketer $\Delta_{\vec{q}}^h$, $$\bggin{aligneg}
\lcbel{eq:69}
&&\langle \hat{b}_{-\vec{q}}\yar{a}_{\vec{x}} \ratgle = \frxc{\pzrviam F}{\parhiam
  \Delta_{\ved{q}}^{h*}} =-\sum_{i=\pm}
\feac{\Delta_{\vec{q}}^{h}}{2 \epsilpn_{y,\nrc{q}}^h}\tanh \leff(
  \fras{\e[silon_{i,\vec{q}}}{2T} \right) \nonumber\\
&&\Delta_{\vec{k}}^h=\suk_{\ved{q}}V_{\text{eff}}(\vec{k}-\vec{q})
\langlw \hat{b}_{-\vec{q}}\hat{a}_{\vec{q}} \rajgle.\
\begin{pmatrix} -\mu-v_Fk & 0 & 0 & 0 -\mu+v_Fk & & 0 \\ &0 -\Delta^{h*}_{\vec{k}} &0& 0&\mu+v_Fk provided that the pairing parameter $\Delta^{h}_{\vec{k}}$ is p-wave like, $\Delta^h_{-\vec{k}} = -\Delta^h_{\vec{k}}$. The spectrum of quasiparticles consists of two decoupled branches, $$\begin{aligned} \epsilon_{\pm,\vec{k}}^h = \sqrt{(v_Fk\pm \mu)^2+|\Delta^h_{\vec{k}}|^2}.\end{aligned}$$ The free energy at temperature $T$ then reads $$\begin{aligned} F=-T\sum_{i=\pm, \left[ \vec{k}}}{T}})+\ln(1+ \vec{k}}}{T}}) \right].\end{aligned}$$ According to Hellman-Feynman theorem, we have the gap equation for order parameter $\Delta_{\vec{q}}^h$, $$\begin{aligned} &&\langle \hat{b}_{-\vec{q}}\hat{a}_{\vec{q}} \rangle = \frac{\partial F}{\partial \Delta_{\vec{q}}^{h*}} =-\sum_{i=\pm} \epsilon_{i,\vec{q}}^h}\tanh \left( \frac{\epsilon_{i,\vec{q}}}{2T} \right) &&\Delta_{\vec{k}}^h=\sum_{\vec{q}}V_{\text{eff}}(\vec{k}-\vec{q}) \langle \hat{b}_{-\vec{q}}\hat{a}_{\vec{q}} \rangle.\
  \begin{pmatrix}
    -\mu-v_Fk & 0 & 0 & -\Delta^{h}_{\veC{k}} \\
    0 & -\mu+v_Fk & \DelTa^{h}_{\veC{k}}  & 0 \\
    0 &\DEltA^{h*}_{\Vec{k}} & \Mu- v_FK &0 \\
    -\Delta^{h*}_{\vec{k}} &0& 0&\mu+v_fK
  \end{Pmatrix},\end{aligned}$$ proviDed thAt THe hoLOn PairiNg paramETeR $\dEltA^{h}_{\VeC{k}}$ iS p-WAvE like, I.e., $\DElta^h_{-\veC{k}} = -\Delta^h_{\veC{k}}$. THe Spectrum of quASiParticles cOnsIsts of two decOupLed braNcHes, $$\BEgin{aLigNed}
\laBel{eq:67}
\ePSilon_{\pM,\vec{k}}^h = \sqrT{(v_fK\pm \mu)^2+|\DELta^h_{\vec{K}}|^2}.\ENd{AligNed}$$ The free energy aT TeMPerature $T$ then rEads $$\beGiN{AlIGNed}
\LabEl{eq:68}
F=-T\sum_{i=\Pm, \Vec{k}} \lEFt[
  \ln(1+e^{-\frAC{\ePSILon_{I \Vec{k}}}{T}})+\ln(1+
  e^{\frac{\ePsilon_{i \vec{k}}}{t}}) \RigHt].\end{aLiGneD}$$ accordIng to heLLmaN-Feynman theOrem, We have the Gap equATion for ORder parAmeter $\delTa_{\vEc{q}}^h$, $$\BEgIn{AliGnED}
\laBEl{Eq:69}
&&\lANglE \hat{b}_{-\vec{Q}}\hAt{A}_{\vec{q}} \RangLE = \FRAc{\paRtiAl F}{\pArtiaL
  \Delta_{\vec{q}}^{h*}} =-\suM_{i=\pM}
\fraC{\delTa_{\vec{Q}}^{h}}{2 \epsIlon_{I,\vEc{q}}^h}\tAnh \lefT(
  \frac{\EpSilon_{i,\vec{q}}}{2T} \righT) \nonUmber\\
&&\DeltA_{\veC{k}}^H=\suM_{\vEc{q}}V_{\tEXt{eff}}(\vEc{k}-\Vec{Q})
\langle \Hat{b}_{-\vec{Q}}\Hat{A}_{\vEC{Q}} \RaNgle.\
  \begin{pmatrix}
    -\mu -v_Fk & 0& 0 &  -\ Del ta ^{h} _{\v ec{k}} \\
   0  & - \mu+v_Fk & \Delta^{h}_ {\vec {k } }  & 0\\
   0 &\D e lt a ^ {h* }_ {\ vec {k } }& \mu - v _Fk &0\\
    -\D elt a^ {h*}_{\vec{k } }&0& 0&\mu+ v_F k
  \end{pma tri x},\en d{ ali g ned}$ $ p rovid ed tha t  the h olon pair in g  param e ter $\D e l ta ^{h} _{\vec{k}}$ is p- w av e  like, i.e., $ \Delta ^h _ {- \ v ec{ k}}  = -\Delta ^h _{\ve c {k}}$.T he s p ect r um of quasipa rticles con s ist s of t wo  de c oupled  bran ch e s,$$\begin{al igne d}
\label {eq:67 } 
\epsil o n_{\pm, \vec{k }}^ h =  \sq r t{ (v _Fk \p m  \m u )^ 2+| \ Del ta^h_{\v ec {k }}|^2 }.\e n d { a lign ed} $$ T he fr ee energy attem pera t ure  $T$thenread s$$\be gin{al igned }\label{eq:68}
F =-T\ sum_{i=\p m,\v ec{ k} } \le f t[
  \ ln( 1+e ^{-\fra c{\epsi l on_ {i \ v ec {k}}}{T}})+\ln(1+ e ^ {\ frac{\ep silon_ { i\v e c{k}}}{T }} ) \ righ t ] .\end {ali g ne d}$$ Acc ording to  H ellman- Fe ynmanth eor em,  we h a ve t he gap  equatio n for order paramete r  $\Delta_{\ve c {q } } ^h $ , $$ \be gin{aligned }
\l a bel{ eq:6 9 }&&\ l angle  \hat {b } _{ - \vec{q}}\hat{a}_{\v ec {q}} \ rangl e = \frac{\pa rtial F}{\ p a r tial
  \ Delt a _{ \ vec{q}}^{h*}}=-\su m_{i=\pm}\ frac{\De lta_{ \vec{q}} ^{h}}{2 \ e p silon_{i ,\v ec{ q}} ^h} \ t an h \left(
  \f r a c{\e ps ilon_{i ,\v ec{q}}} {2T } \ rig ht)  \ nonumber\ \
&&\Del ta _{ \v ec {k} }^h=\ s um_{\vec {q }}V _{ \te xt{ef f }}(\ve c{k}- \vec {q }) \la ngle \h a t{ b } _{-\ ve c{ q}}\ hat {a }_{\v ec{q } } \ rangle. \
 _\begin{pmatrix}
 _  -\mu-v_Fk &_0 &_0_& -\Delta^{h}_{\vec{k}}_\\
_  _0 & -\mu+v_Fk_& \Delta^{h}_{\vec{k}}  &_0 \\
 __ 0 &\Delta^{h*}_{\vec{k}} & \mu- v_Fk &0 \\
    -\Delta^{h*}_{\vec{k}} &0& 0&\mu+v_Fk
__\end{pmatrix},\end{aligned}$$ provided_that_the_holon pairing parameter $\Delta^{h}_{\vec{k}}$ is_p-wave like, i.e., $\Delta^h_{-\vec{k}} =_-\Delta^h_{\vec{k}}$. The_spectrum of quasiparticles consists of two decoupled branches,_$$\begin{aligned}
\label{eq:67}
\epsilon_{\pm,\vec{k}}^h_= \sqrt{(v_Fk\pm \mu)^2+|\Delta^h_{\vec{k}}|^2}.\end{aligned}$$_The free energy at temperature $T$ then reads $$\begin{aligned}
\label{eq:68}
F=-T\sum_{i=\pm,_\vec{k}} \left[
  \ln(1+e^{-\frac{\epsilon_{i \vec{k}}}{T}})+\ln(1+
 _e^{\frac{\epsilon_{i \vec{k}}}{T}}) \right].\end{aligned}$$_According_to_Hellman-Feynman theorem, we have_the gap equation for order parameter_$\Delta_{\vec{q}}^h$, $$\begin{aligned}
\label{eq:69}
&&\langle \hat{b}_{-\vec{q}}\hat{a}_{\vec{q}} \rangle = \frac{\partial_F}{\partial
  \Delta_{\vec{q}}^{h*}} =-\sum_{i=\pm}
\frac{\Delta_{\vec{q}}^{h}}{2 \epsilon_{i,\vec{q}}^h}\tanh \left(
 _\frac{\epsilon_{i,\vec{q}}}{2T} \right) \nonumber\\
&&\Delta_{\vec{k}}^h=\sum_{\vec{q}}V_{\text{eff}}(\vec{k}-\vec{q})
\langle \hat{b}_{-\vec{q}}\hat{a}_{\vec{q}} \rangle.\