#### Abstract

We examine the thermodynamic features of an exotic fluid known as modified cosmic Chaplygin gas in the context of homogeneous isotropic universe model. For this purpose, the behavior of physical parameters is discussed that help to analyze nature of the universe. Using specific heat formalism, the validity of third law of thermodynamics is checked. Furthermore, with the help of thermodynamic entities, the thermal equation of state is also discussed. The thermodynamic stability is explored by means of adiabatic, specific heat and isothermal conditions from classical thermodynamics. It is concluded that the considered fluid configuration is thermodynamically stable and expands adiabatically for an appropriate choice of parameters.

#### 1. Introduction

The discovery of accelerated expansion of the universe has unambiguously been proved by a diverse set of high-precision observational data accumulated from various astronomical sources [1–3]. Dark energy (DE) is considered as the root cause behind this tremendous change in cosmic history. It possesses negatively large pressure which violates the strong energy condition (, where and are pressure and energy density, respectively) but its complete characteristics are still not known. Planck’s observational data reveals that about % of our universe is filled with this mysterious form of energy while the other cosmic budget includes % ordinary matter and % dark matter (DM) [4, 5]. To explore the perplexing nature of DE, there began a search for different candidates that can play their role as an alternative for DE. The simplest candidate is the cosmological constant while other favorable approaches include quintessence, k-essence, and Chaplygin gas (CG) known as DE matter models [6–8].

Chaplygin gas is an intriguing model presented by a Russian physicist Chaplygin as a convenient soluble model to study the lifting force on the wing of an aeroplane in aerodynamics. It efficiently describes the cosmic expansion and elegantly discusses DM and DE in a unified form. The distinct feature of this model is its positive and bounded squared speed of sound leading to stable results as compared to other fluids with negative pressure. Chaplygin gas acts as an alternative for dust dominated era with small value of the scale factor and tends to cosmic expansion for its large value while primordial universe cannot be discussed in this scenario [9–12]. Despite the fact that it does not meet the strong energy condition, it successfully shows consistency with the observational results accumulated from various cosmic probes such as Hubble space telescope, Wilkinson microwave anisotropy probe, and cosmic background explorer.

To discuss the cosmic history as well as get more accuracy with observational data, several modifications of CG have been presented which are obtained by introducing new parameters in its equation of state (EoS). Bento et al. [13] established the generalization of this model named generalized Chaplygin gas (GCG) which interprets the same evolutionary phases of the universe as CG. Benaoum [14] introduced the modified Chaplygin gas (MCG) which illustrates the radiation dominated era. González-Díaz [15] proposed the generalization of cosmic CG models known as generalized cosmic Chaplygin gas (GCCG) in such a way that it avoids big-rip (singularity at a finite time) which was previously presented in the DE models representing phantom era. This generalization provides stable and physical behavior models even when the vacuum matter configuration fulfills the phantom energy condition (, , , where is the EoS parameter). The other proposed CG models include variable Chaplygin gas (VCG), variable modified Chaplygin gas (VMCG), new variable MCG, extended CG, and modified cosmic Chaplygin gas (MCCG) [16–19].

Chaplygin gas models have stimulated many researchers to investigate their thermal stability. Santos et al. [20, 21] explored thermal stability of GCG as well as MCG and deduced that these fluid models verify the third law of thermodynamics along with the adiabatic expansion. Myung [22] proved the third law of thermodynamics for CG model and illustrated that it can represent a unified picture of DM and DE without any phase transition. Kahya and Pourhassan [23] analyzed extended CG model cosmologically as well as thermodynamically and found stable results without any phase transition against density perturbations. They also concluded that all laws of thermodynamics are satisfied for this exotic fluid which came out to be thermodynamically stable at all times. Panigrahi [24] studied thermodynamic behavior of VCG and observed that it is thermodynamically stable throughout the evolution. He also discussed thermal EoS which is an explicit function of temperature only and checked the validity of the third law. Panigrahi and Chatterjee [25] found that current accelerated expansion of the universe can be explained using VMCG model. Sharif and Sarwar [26] explored how GCCG can explain accelerated expansion of the universe by interpreting different physical parameters and showed that the fluid is adiabatically stable.

Here, we investigate thermodynamic stability of MCCG model in the background of isotropic and homogeneous universe model. In Section 2, we discuss the behavior of physical parameters such as pressure and EoS as well as deceleration parameters and analyze the stability using speed of sound. Section 3 deals with the thermodynamic stability of MCCG. The results are summarized in the last section.

#### 2. Physical Parameters for MCCG

In this section, we discuss the behavior of MCCG in the background of FRW universe model for different physical parameters and examine its stability through squared speed of sound. The line element for FRW universe model is given bywhere is the scale factor. The EoS for MCCG is defined aswhere ( is an arbitrary constant) and is a positive constant. This EoS reduces to GCCG as [26] while GCG is recovered in the limit along with [13]. The energy density of the fluid configuration is given bywhere and represent the internal energy and volume, respectively. Classical thermodynamics provides a useful relationship among the quantities , , and in the formUsing (2) in the above expression, we obtainwhich is a nonlinear ordinary differential equation. Its solution is given bywhere we have used the binomial expansion up to first order and is an integration constant which is either universal constant or a function of entropy . The above equation can also be written aswhere , , and . It is clearly observed that internal energy of the fluid can only be discussed when is a whole number between the above-mentioned range for positive values of . Using (3) and (7), the energy density of MCCG becomes

In the following, we use this equation to discuss different physical parameters.

##### 2.1. Pressure

The pressure of MCCG in terms of can be obtained using (2) and (8) asThe graphical analysis of this equation is shown in Figure 1 for different values of with . The positive and negative behavior of pressure correspond to decelerated and accelerated eras of the universe, respectively. For , it is observed that the accelerating universe at small volume tends to dust dominated universe at large volume. We note that the decelerating universe tends to accelerate as volume increases for , whereas accelerated phase is obtained only for . The same behavior of pressure is observed for different values of the parameter .

##### 2.2. EoS Parameter

Here we discuss the effective EoS parameter of MCCG. Using (8) and (9), we haveWe study the following two extremal cases for volume to analyze the behavior of above equation. (i)For small volume , the above equation reduces towhich is a barotropic EoS. In this case, will depend entirely on the value of .(ii)For large volume , (10) takes the formFor , let volume be denoted by , which is given by

The EoS parameter discusses both accelerated and decelerated phases of the universe and successfully describes the phase transitions (dubbed as flip) at a critical value between these cosmic phases. The proper flip occurs when while the inequality leads to real flip.

Figure 2 shows the behavior of for different values of . It is found that the value of decreases as becomes negatively large. The negative values of show the decelerated cosmic phase at small volume undergoing acceleration at large volume while demonstrates only the decelerating phase. We observe that, at large volume, the considered values of , and correspond to stiff matter, dust dominated, CDM, and phantom eras, respectively. Thus, the EoS parameter can interpret different evolutionary phases of the universe.

##### 2.3. Deceleration Parameter

The deceleration parameter is given byUsing (10), this parameter for MCCG takes the formFor small volume, it reduces towhich implies that the universe undergoes deceleration at its early stage since while for large volume, (15) becomesIn this case, the flip occurs when deceleration parameter vanishes and the corresponding flip volume is given byprovided that and the inequality leads to proper flip. Figure 3 shows the evolution of deceleration parameter against volume for different values of . At small volume, the universe undergoes deceleration while accelerating behavior is observed at large volume for considered negative values of . The flip occurs at and 3 for and , respectively. For , the deceleration parameter switches from acceleration to deceleration at while no flip is observed for .

##### 2.4. Speed of Sound

Here we analyze the stability of MCCG using speed of sound aswhose feasible range is . This equation reduces to at early universe while, for , we haveFigure 4 shows the behavior of squared speed of sound against the positive parameter for different values of . It is observed that the viable ranges for are , , and corresponding to , 4, and 7, respectively. Thus, the stable results are found for the considered values of in a particular range of .

#### 3. Thermodynamic Stability

In this section, we discuss thermodynamic stability of MCCG during its evolution. The stability conditions are given by the following [27]:(i)The pressure reduces for both adiabatic as well as isothermal expansions aswhere represents temperature.(ii)Specific heat at constant volume is positive.

Differentiation of (9) with respect to volume yieldsWhen volume is very small, the above equation reduces to zero while, for large volume, we have the following expression:Figure 5 shows that the adiabatic condition is fulfilled for all the considered values of .

To investigate the positivity of specific heat at constant volume, we consider specific heat in terms of temperature and entropy aswhere the temperature of MCCG is obtained from the following relation:Differentiating (6) with respect to , we haveSubstituting this relation in (25), the expression of becomesWhen is a universal constant (), the temperature vanishes while it varies for CG expansion, so we consider . Here, we assume the case to have a positive temperature which is cooled down through adiabatic expansion. Using the concept of dimensional analysis, (6) givesUsing the relation , the above equation becomesTaking , it follows thatwhere and are constants having the dimensions of temperature and volume, respectively. Differentiating (30) with respect to , we obtainSubstituting this value in (27), the temperature of MCCG takes the formwhere . Using (6) and (30) in the above equation, we have

For a positive definite entropy, we assume and . It is worth mentioning here that when , the entropy vanishes which indicates that the considered fluid obeys third law of thermodynamics (if the temperature of a physical system approaches zero then entropy becomes zero). Differentiating (33) with respect to , we obtainwhere . Inserting (33) and (34) in (24), it follows thatFor to be real, and to be positive, the inequality must hold. The graphical analysis of (35) is shown in Figure 6 for different values of with . We observe that the positivity of specific heat is obtained for in the range while MCCG is thermally stable for both values of and throughout the evolution. It is also noted that when temperature is zero, thermal capacity vanishes which also assures the validity of third law of thermodynamics.

Finally, we analyze the behavior of considered model through isothermal condition. For this purpose, we assume and by solving (9) and (30), we haveThe corresponding EoS parameter takes the formTo check the isothermal condition, we should have and . In our case, it is difficult to have a thermal EoS for MCCG as a function of temperature only since (33) is a complicated equation such that the explicit expression for in terms of cannot be extracted. For this reason, we are unable to analyze the isothermal condition in this scenario.

#### 4. Conclusions

In this paper, we have analyzed thermodynamic stability of MCCG within the framework of FRW universe model. We have examined the expanding evolution of the universe through different physical parameters like pressure, effective EoS, and deceleration parameters as well as speed of sound. The results of these parameters can be summarized as follows:(i)The consistent behavior of pressure with the evolutionary picture of the universe is obtained for the considered values of whereas inconsistent evolution is observed for its positive values (Figure 1).(ii)The EoS parameter for MCCG depicts that decelerated and accelerated phases of our universe can be discussed for different values of parameter (Figure 2). We have also calculated the critical value at and found that its value increases as increases from its negative values to zero.(iii)The evolution of deceleration parameter against volume gives the decelerated universe when for negative values of while accelerating behavior is observed when . For positive values of , we have found only deceleration while at , acceleration occurs before the flip (Figure 3).(iv)We have analyzed the stability of MCCG through speed of sound and obtained stable regions at large volume for the considered values of (Figure 4). For VMCG model, the squared speed of sound could be positive or negative [25] whereas for GCCG model, the stable regions do not exist in late universe [26].

Finally, we have investigated thermodynamic stability of considered fluid configuration using adiabatic, isothermal, and specific heat conditions. We have found the validity of adiabatic as well as positivity of specific heat for the considered values of (Figures 5 and 6). It is worth mentioning here that third law of thermodynamics is obeyed for MCCG. We conclude that MCCG expands adiabatically and the expansion is thermodynamically stable for a suitable choice of the parameters.

#### Data Availability

No data were used to support this study.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.