tilde{E}(\tilde{T})-\tilde{E}_{bp}$ for the parameter values $\tilde{\omega}_i$ presented in Fig. \[fig1\], \[fig2\]. \[fig3\]](fig3_2.eps "fig:"){width="80.00000%"}

The dependencies obtained enable us to find the heat capacity of a TI bipolaron gas: $C_v(\tilde{T})=d\tilde{E}/d\tilde{T}$. With the use of (20) $C_v(\tilde{T})$ for $\tilde{T}\leq\tilde{T}_c$ is expressed as:

$$\label{22}
C_v(\tilde{T})=\frac{\tilde{T}^{3/2}}{2C_{bp}}
\left[\frac{\tilde{\omega}^2}{\tilde{T}^2}F_{1/2}\left(\frac{\tilde{\omega}}{\tilde{T}}\right)+
6\left(\frac{\tilde{\omega}}{\tilde{T}}\right)F_{3/2}\left(\frac{\tilde{\omega}}{\tilde{T}}\right)+
5F_{5/2}\left(\frac{\tilde{\omega}}{\tilde{T}}\right)\right],$$

$$F_{1/2}(\alpha)=\frac{2}{\sqrt{\pi}}\int^{\infty}_0\frac{1}{\sqrt{x}}\ \frac{dx}{e^{x+\alpha}-1}$$

Expression (22) yields a well-known exponential dependence of the heat capacity at low temperatures $C_v\backsim\text{exp}(-\omega_0/T)$, caused by the availability of the energy gap $\omega_0$.

![ Temperature dependencies of the heat capacity for various values of the parameters $\omega_i$: $\omega_0=0$; $\ \tilde{T}_{C_0}=25,2$; $\ C_v(\tilde{T}_{c1})=2$; $\omega_1=0,2$; $\ \tilde{T}_{C_1}=27,3$; $\ C_v(\tilde{T}_{c1}-0)=2,16$; $\ C_v(\tilde{T}_{c1}+0)=1,9$; $\omega_2=1$; $\ \ \tilde{T}_{C_2}=30$; $\ C_v(\tilde{T}_{c2}-0)=2,46
tilde{E}(\tilde{T})-\tilde{E}_{bp}$ for the parameter values $ \tilde{\omega}_i$ presented in Fig.   \[fig1\ ], \[fig2\ ]. \[fig3\]](fig3_2.eps " fig:"){width="80.00000% " } 

 The dependencies prevail enable us to receive the heat capability of a TI bipolaron flatulence: $ C_v(\tilde{T})=d\tilde{E}/d\tilde{T}$. With the use of (20) $ C_v(\tilde{T})$ for $ \tilde{T}\leq\tilde{T}_c$ is expressed as: 

 $ $ \label{22 } 
 C_v(\tilde{T})=\frac{\tilde{T}^{3/2}}{2C_{bp } } 
 \left[\frac{\tilde{\omega}^2}{\tilde{T}^2}F_{1/2}\left(\frac{\tilde{\omega}}{\tilde{T}}\right)+ 
 6\left(\frac{\tilde{\omega}}{\tilde{T}}\right)F_{3/2}\left(\frac{\tilde{\omega}}{\tilde{T}}\right)+ 
 5F_{5/2}\left(\frac{\tilde{\omega}}{\tilde{T}}\right)\right],$$ 

 $ $ F_{1/2}(\alpha)=\frac{2}{\sqrt{\pi}}\int^{\infty}_0\frac{1}{\sqrt{x}}\ \frac{dx}{e^{x+\alpha}-1}$$ 

 Expression (22) yield a well - know exponential dependence of the heat capacity at low temperature $ C_v\backsim\text{exp}(-\omega_0 / T)$, caused by the availability of the department of energy gap $ \omega_0$. 

! [ Temperature dependencies of the heat capacity for various values of the parameters $ \omega_i$: $ \omega_0=0 $; $ \ \tilde{T}_{C_0}=25,2 $; $ \ C_v(\tilde{T}_{c1})=2 $; $ \omega_1=0,2 $; $ \ \tilde{T}_{C_1}=27,3 $; $ \ C_v(\tilde{T}_{c1}-0)=2,16 $; $ \ C_v(\tilde{T}_{c1}+0)=1,9 $; $ \omega_2=1 $; $ \ \ \tilde{T}_{C_2}=30 $; $ \ C_v(\tilde{T}_{c2}-0)=2,46
tilfe{E}(\tilde{T})-\tilde{E}_{bp}$ for tht parameter values $\tilde{\mmega}_i$ presentdd in Fig. \[fig1\], \[fig2\]. \[fig3\]](fig3_2.eps "fmg:"){wieth="80.00000%"}

Tht dependencies obtakned enabpe us to finv the heat capacmfy of a TI billlarmi gas: $C_v(\tilde{T})=d\jilde{E}/d\tilde{D}$. With the use ow (20) $C_v(\tilde{T})$ for $\tilde{T}\leq\tilde{T}_c$ is evpressec ws:

$$\label{22}
C_v(\tilde{J})=\frac{\uilqe{T}^{3/2}}{2C_{gi}}
\ltft[\frac{\tilde{\omega}^2}{\tilde{T}^2}F_{1/2}\left(\frac{\tjlde{\omeja}}{\tilde{T}}\right)+
6\legt(\frac{\tilde{\omega}}{\tilde{T}}\righh)F_{3/2}\levt(\frac{\tilde{\omega}}{\tipde{T}}\right)+
5F_{5/2}\lgrt(\fwqc{\tilde{\omega}}{\gilde{T}}\righu)\rnght],$$

$$F_{1/2}(\alpha)=\ftac{2}{\sqrt{\pi}}\int^{\infty}_0\frac{1}{\sqrt{x}}\ \frac{db}{e^{x+\al'ha}-1}$$

Expressiin (22) ylglds a well-kiown eqponential deivndence of the heat capacity at loq temperatures $C_v\backvim\text{exp}(-\omega_0/T)$, cwused by dhz availability of the ebergy gap $\omeeq_0$.

![ Tdmptravurs depejdeicies of ths heat capaxity for various vakuqw of the parajeters $\oiega_i$: $\omega_0=0$; $\ \tilde{T}_{C_0}=25,2$; $\ C_v(\tilde{T}_{c1})=2$; $\omega_1=0,2$; $\ \uilde{F}_{C_1}=27,3$; $\ C_v(\tilde{T}_{c1}-0)=2,16$; $\ C_v(\tilde{T}_{x1}+0)=1,9$; $\omega_2=1$; $\ \ \tilde{T}_{C_2}=30$; $\ C_v(\jilde{T}_{c2}-0)=2,46
tilde{E}(\tilde{T})-\tilde{E}_{bp}$ for the parameter values $\tilde{\omega}_i$ presented \[fig1\], \[fig3\]](fig3_2.eps "fig:"){width="80.00000%"} dependencies obtained enable capacity a TI bipolaron $C_v(\tilde{T})=d\tilde{E}/d\tilde{T}$. With the of (20) $C_v(\tilde{T})$ for $\tilde{T}\leq\tilde{T}_c$ is as: $$\label{22} C_v(\tilde{T})=\frac{\tilde{T}^{3/2}}{2C_{bp}} \left[\frac{\tilde{\omega}^2}{\tilde{T}^2}F_{1/2}\left(\frac{\tilde{\omega}}{\tilde{T}}\right)+ 6\left(\frac{\tilde{\omega}}{\tilde{T}}\right)F_{3/2}\left(\frac{\tilde{\omega}}{\tilde{T}}\right)+ 5F_{5/2}\left(\frac{\tilde{\omega}}{\tilde{T}}\right)\right],$$ $$F_{1/2}(\alpha)=\frac{2}{\sqrt{\pi}}\int^{\infty}_0\frac{1}{\sqrt{x}}\ \frac{dx}{e^{x+\alpha}-1}$$ Expression (22) yields a well-known dependence of the heat capacity at low temperatures $C_v\backsim\text{exp}(-\omega_0/T)$, caused by the availability the gap ![ dependencies of the heat capacity for various values of the parameters $\omega_i$: $\omega_0=0$; $\ \tilde{T}_{C_0}=25,2$; $\ $\omega_1=0,2$; $\ \tilde{T}_{C_1}=27,3$; $\ C_v(\tilde{T}_{c1}-0)=2,16$; $\ C_v(\tilde{T}_{c1}+0)=1,9$; $\omega_2=1$; \ \tilde{T}_{C_2}=30$; $\ C_v(\tilde{T}_{c2}-0)=2,46
tilde{E}(\tilde{T})-\tilde{E}_{bp}$ for thE parameter ValueS $\tiLde{\OmEga}_i$ PresEnted in Fig. \[fig1\], \[fIG2\]. \[fig3\]](Fig3_2.eps "fig:"){width="80.00000%"}

The depenDenciEs OBtaiNEd EnablE us to fiND tHE HeaT cApAciTy OF a tI bipOlaRon gas: $C_V(\tilde{T})=d\tiLde{e}/d\Tilde{T}$. With thE UsE of (20) $C_v(\tilde{t})$ foR $\tilde{T}\leq\tiLde{t}_c$ is exPrEssED as:

$$\laBel{22}
c_v(\tilDe{T})=\fraC{\Tilde{T}^{3/2}}{2c_{bp}}
\left[\frAc{\TIlde{\omEGa}^2}{\tilde{t}^2}f_{1/2}\LeFt(\frAc{\tilde{\omega}}{\tilde{t}}\RiGHt)+
6\left(\frac{\tildE{\omega}}{\TiLDe{t}}\RIghT)F_{3/2}\lEft(\frac{\tilDe{\Omega}}{\TIlde{T}}\riGHt)+
5f_{5/2}\LEFt(\fRAc{\tilde{\omega}}{\tIlde{T}}\right)\rIGht],$$

$$f_{1/2}(\alpha)=\FrAc{2}{\sQRt{\pi}}\inT^{\inftY}_0\fRAc{1}{\sQrt{x}}\ \frac{dx}{e^{X+\alpHa}-1}$$

ExpressIon (22) yieLDs a well-KNown expOnentiAl dEpeNdenCE oF tHe hEaT CapACiTy aT Low TemperatUrEs $c_v\bacKsim\TEXT{Exp}(-\oMegA_0/T)$, caUsed bY the availabilIty Of thE EneRgy gaP $\omegA_0$.

![ TemPeRaturE depenDenciEs Of the heat capaciTy foR various vAluEs Of tHe ParamETers $\omEga_I$: $\omEga_0=0$; $\ \tildE{T}_{C_0}=25,2$; $\ C_v(\tiLDe{T}_{C1})=2$; $\oMEGA_1=0,2$; $\ \tIlde{T}_{C_1}=27,3$; $\ C_v(\tilde{T}_{c1}-0)=2,16$; $\ C_v(\tIlDE{t}_{c1}+0)=1,9$; $\Omega_2=1$; $\ \ \tilDe{T}_{C_2}=30$; $\ C_v(\TIlDe{t}_{C2}-0)=2,46
tilde{E}(\tilde{T})-\tilde {E}_{bp}$for t hepar am eter  val ues $\tilde{\o m ega} _i$ presented in Fig.\[fig 1\ ] , \[ f ig 2\].\[fig3\ ] ]( f i g3_ 2. ep s " fi g :" ){wid th= "80.000 00%"}

The  de pe ndencies obt a in ed enableusto find thehea t capa ci tyo f a T I b ipola ron ga s : $C_v (\tilde{T }) = d\tild e {E}/d\t i l de {T}$ . With the use of (2 0 ) $C_v(\tilde{ T})$ f or $\ t i lde {T} \leq\tilde {T }_c$i s expre s se d a s:$$\label{22}C_v(\tilde{ T })= \frac{ \t ild e {T}^{3 /2}}{ 2C _ {bp }}
\left[\f rac{ \tilde{\o mega}^ 2 }{\tild e {T}^2}F _{1/2} \le ft( \fra c {\ ti lde {\ o meg a }} {\t i lde {T}}\rig ht )+ 
6\le ft(\ f r a c {\ti lde {\om ega}} {\tilde{T}}\r igh t)F_ { 3/2 }\lef t(\fr ac{\ ti lde{\ omega} }{\ti ld e{T}}\right)+
5 F_{5 /2}\left( \fr ac {\t il de{\o m ega}}{ \ti lde {T}}\ri ght)\ri g ht] ,$ $ $$ F_{1/2}(\alpha)=\f ra c { 2} {\sqrt{\ pi}}\i n t^ {\ i nfty}_0\ fr ac{ 1}{\ s q rt{x} }\ \ f ra c{dx}{e^ {x+\al p ha }- 1}$$

E xp ressio n(22 ) y ields a we ll-kno wn expon entia l  dependence of the heat capa c it y at lowtem peratures $ C_v\ b acks im\t e xt {ex p }(-\o mega_ 0/ T )$ ,  caused by the avai la bility  of t he energy gap  $\omega_0 $ . 
![ Temp erat u re dependencies o f the  heat capa c ity forvario us value s of thep a rameters  $\ ome ga_ i$: $ \o mega_0=0$; $\ \ tild e{ T}_{C_0 }=2 5,2$; $ \ C _v( \ti lde {T }_{c1})=2 $; $\ome ga _1 =0 ,2 $;$\ \t i lde{T}_{ C_ 1}= 27 ,3$ ; $\C _v(\ti lde{T }_{c 1} -0 ) =2, 16$; $\ C_ v ( \til de {T }_{c 1}+ 0) =1,9$ ; $\ o meg a_2=1$;  $\ \ \ti lde { T}_{ C_ 2} =30$; $ \ C_v(\tilde{ T} _{c2}-0)=2 ,4 6
tilde{E}(\tilde{T})-\tilde{E}_{bp}$ for_the parameter_values $\tilde{\omega}_i$ presented in_Fig. \[fig1\], \[fig2\]._\[fig3\]](fig3_2.eps_"fig:"){width="80.00000%"}

The dependencies_obtained_enable us to_find the heat_capacity of a TI_bipolaron gas: $C_v(\tilde{T})=d\tilde{E}/d\tilde{T}$._With_the use of (20) $C_v(\tilde{T})$ for $\tilde{T}\leq\tilde{T}_c$ is expressed as:

$$\label{22}
C_v(\tilde{T})=\frac{\tilde{T}^{3/2}}{2C_{bp}}
\left[\frac{\tilde{\omega}^2}{\tilde{T}^2}F_{1/2}\left(\frac{\tilde{\omega}}{\tilde{T}}\right)+
6\left(\frac{\tilde{\omega}}{\tilde{T}}\right)F_{3/2}\left(\frac{\tilde{\omega}}{\tilde{T}}\right)+
5F_{5/2}\left(\frac{\tilde{\omega}}{\tilde{T}}\right)\right],$$

$$F_{1/2}(\alpha)=\frac{2}{\sqrt{\pi}}\int^{\infty}_0\frac{1}{\sqrt{x}}\ \frac{dx}{e^{x+\alpha}-1}$$

Expression (22) yields a_well-known_exponential dependence_of_the_heat capacity at low temperatures_$C_v\backsim\text{exp}(-\omega_0/T)$, caused by the availability_of the_energy gap $\omega_0$.

![ Temperature dependencies of the heat_capacity_for various values_of the parameters $\omega_i$: $\omega_0=0$; $\ \tilde{T}_{C_0}=25,2$; $\ C_v(\tilde{T}_{c1})=2$;_$\omega_1=0,2$; $\ \tilde{T}_{C_1}=27,3$; $\ C_v(\tilde{T}_{c1}-0)=2,16$; $\_C_v(\tilde{T}_{c1}+0)=1,9$; $\omega_2=1$; $\_\_\tilde{T}_{C_2}=30$;_$\ C_v(\tilde{T}_{c2}-0)=2,46